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The high voltage x-radiation from the betatron has been employed to produce (y, ”) reactions 
in several elements of atomic number up to 47, the reaction observed by detection of induced 
radioactivity. The peak x-ray energy was controlled by integrating the voltage on the main 
coils with an RC circuit, which actuated the orbit expander at a predetermined electron energy. 
Upon irradiating samples at sequences of energies, measuring their activity with beta-counters, 
and plotting activities against peak energy, smooth curves were obtained from which thresholds 
were estimated. The accuracy of the method is limited by low activity in almost all cases and by 
uncertainties in energy output. The following thresholds in Mev have been measured: C", 18.7 
to 19.4; 11.1+0.5; 16.3+0.4; Fe®, 14.2+0.4; Cu®, 10.9+0.3; Zn™, 114+0.4; Se” o 
9.8+0.5 for the lower, short period isomer; Mo” %, 13.5+0.4; Ag!*, 9.3+0.5; 9.5. The 
first three are in fair agreement with other data, but because of exceptionally weak activities are 
not sufficiently reliable to afford a test of the method. Rough measurements on the uranium 
photo-fission threshold are also reported. Since the radiation from the betatron is continuous and 
its spectral distribution is unknown, the excitation curves are not susceptible to interpretation, 
and cross sections for the above reactions cannot be estimated at present. 


fying the particular isotopes of the target element 
involved in the reaction. 

To produce the photo-effect it is necessary that 
the energy of the incident quantum exceed the 
binding energy of the particle to be ejected. The 
following reaction is typical: 


INTRODUCTION 


ITH the high energy x-rays now available 

from the betatron, it is energetically pos- 

sible to photo-disintegrate nearly all nuclei. In all 
cases so far observed, photo-disintegration of a 
target nucleus of charge Z and mass number A 
leads to an isotope of the target with mass 
number A—1, upon ejection of a neutron from 
the target nucleus. The process may be detected 


by direct observation of the ejected neutron; 
however, frequently the product nucleus is 
radioactive and the reaction is then most con- 
veniently detected by observation of the induced 
activity in an irradiated sample—this method 
having the added advantage of definitely identi- 


*Now at General Electric Research Laboratory, 
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Conservation of energy requires that 
hy> M(6C*) + M(on!) 


where M denotes the atomic mass. Substitution 
of numerical values’? for the masses of N“ and 


1R. T. Birge, Rev. Mod. Phys. 13, 233 (1941). 
2 i J. Livingood and G. T. Seaborg, Rev. Mod. Phys. 12, 
30 (1940). 
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Fic. 1. Circuit scheme for control of peak x-ray energy. 


C* and for Q leads to 
hy=10.65+0.15 Mev 


for the threshold of the photo-effect in N'*. Simi- 
lar calculations give thresholds of 15.6++0.6 for 
O' and 18.7+0.3 for C. Calculation of the 
threshold is possible only in several cases, how- 
ever. Since a knowledge of the binding energies of 
particles in nuclei is of considerable theoretical 
interest, we have undertaken measurements of 
photo-effect thresholds in a sequence of nuclei 
from Z=6 to Z=43. This work has been reported 
elsewhere.’ 

Previous work on the photo-effect, except in 
the cases of beryllium and deuterium, has been 
limited by the unavailability of x-rays of suffi- 
ciently high energy. For heavier nuclei at least 8 
Mev are required; in all cases reported herein, 
considerably higher energies are necessary. 
Gamma-rays from proton bombardment of lith- 
ium and boron have been the chief photo-dis- 
integrative agents employed in past work. Bothe 
and Gentner* have measured the cross section for 


the reaction Cu®(y, 2)Cu® with lithium and 


*G. C. Baldwin and H. W. Koch, Phys. Rev. 63, 59A 


and 426A (1943). 
‘W. Bothe and W. Gentner, Zeits. f. Physik 112, 45 


(1939). 
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. toring of x-ray intensity, and (3) detection of the 


5D. W. Kerst, Phys. Rev. 60, 47 (1941); Rev. Sci. Inst. 


boron gamma-rays, and have compared the yields 
induced by both radiations in other elements, 
Inhomogeneity of the boron radiation, together 
with low yields of activity, makes interpretation 
of their results uncertain. Thresholds could not be 
estimated in this sort of measurement. 


METHOD 


The threshold measurements reported in this . 
paper were made by irradiating samples at a 
sequence of energies, the energy being raised in 
small steps until the process was observed to © 
occur, as evidenced by radioactivity induced in 
the sample. The yield of activity was plotted asa 
function of peak x-ray energy, giving a curve 
from which the threshold was obtained by in- 
spection. This method involved the problems of 
(1) accurate control of peak energy, (2) moni- 


reaction. 

For energy control, use was made of the fact 
that the energy acquired by electrons during ac- 
celeration in the betatron is proportional to the 
magnetic flux linking their orbit, which is in turn 
proportional to the flux linking the main coils.* A 
resistor and capacitor were placed in series across 
one of the main coils, the time constant large 
compared with the period of the alternating flux. 
The condenser voltage, proportional to the inte- 
gral of the applied voltage, was proportional to 
the flux linkage and hence to the energy. This 
voltage was applied in series with an adjustable 
bias to a trigger circuit to actuate the orbit ex- 
pander at a chosen critical input voltage read on 
the bias voltmeter. The block diagram of Fig. 1 
illustrates the circuit scheme. With the constants 
employed (R= 10° ohms, C= 10-7 farad) a voltage 
of 88.5 v on the integrating condenser corre- 
sponded to an electron energy of 20 Mev. The 
trigger point, reproducible to within 0.2 v, 
showed negligible drift after adequate warm-up. 
All voltages were electronically stabilized. 

This circuit was calibrated by a search coil in 
the orbit plane. It was known that the magnet 
flux is an accurate sine wave, the third harmonic 
amounting to only 0.07 percent at the highest 
amplitude. The flux within the search coil could 
thus be readily determined by a voltage measure- 


13, 387 (1942). 
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ment and simultaneously compared with the 
voltage integrator response and with the magnet 
voltage read on a 2-turn coil about the magnet 
yoke. The latter reading was made a permanent 
substandard to which each setting was referred 
by determining the magnet voltage at which the 
integrator circuit would cause orbit expansion at 
peak flux. From these calibration data the Hr 
value at which the electrons struck the target 
could be computed. The probable error of this 
calibration was 1.5 percent; in actual energy 
settings further error arose. The total error in 
most cases did not exceed 2 percent. 

Samples were irradiated at a standard position 
in the x-ray beam, 46 cm from the target, ordi- 
narily to half-saturation. Irregularities in the re- 
sulting activity could be largely accounted for by 


- variations in x-ray intensity from run to run. The 


intensity, monitored by an ionization chamber, 
was recorded in roentgens per minute at 1 meter 
from the target. The method of correction to 
standard intensity will be evident from a con- 
sideration of the process of radioactive growth 
under decay. During an infinitesimal time inter- 
val dt, the number of active nuclei N will change 
by 

dN=kIdt—Ndt/r (1) 


where k is a constant, J the x-ray intensity, and r 
the mean life of the expected activity. This gives 
for the number of active nuclei produced between 
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Fic. 2. Excitation curve for production of C" by x-rays. 
The cross-hatched area represents the statistical spread in 
natural background of the Geiger-Miiller counter. 
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Fic. 3. Excitation curve for N®: (A) Corrected for intensity 
variations; (B) Uncorrected points near threshold. 


t=0 and t=T 
f I(t) exp (t/r)dt=kR. (2) 
0 


In practice the integral was replaced by a 
weighted sum ; the monitor was read at regular in- 
tervals during irradiation, each reading multiplied 
by the weighting factor, exp (¢/r)-exp (—T7/r), 
and the results added. The resulting sum R is 
called the “‘irradiation.”’ 

To reduce as far as possible the effect of in- 
tensity fluctuations all points were then corrected 
to the irradiation which the sample ordinarily re- 
ceived near threshold. This did not change points 
near threshold, where the effects were small or 
doubtful, but smoothed the curves somewhat and 
thereby facilitated threshold estimation. 

The resulting smoothed curves _are not sus- 
ceptible to ready interpretation. It is not yet 
possible to determine from them the trend of the 
cross section for the process. That portion of the 
x-ray spectrum which is responsible for the reac- 
tion is responsible for only a small portion of the 
monitor response. If the cross section of the 
process be given by o(£) and the unknown spec- 
tral distribution by N(E), the yield at some peak 
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replotted to larger scale in the inset. The first two points 130 40 ao aaa cou 
are below background because of the shielding effect of the MEV 
water sample in which the counter was immersed. Foc. 5. Excitation curve for Fe*, Th 
wit 
energy Emax is proportional to and beta-decay data. These reactions would pro- cov 
i vide a check on the accuracy of the method, ex- pla 
N(E)o(E)dE, cept that the resulting activities are the weakest Na 
Eo of all those studied. All other elements give I 
higher yields and sharper thresholds, but reliable tot 
where E> is the threshold for the reaction. The mass and beta-decay data are not available for rect 
monitor response is given by a similar expression, comparison. plo 
but the lower limit and the energy-dependence of diff 
| the integrand are entirely different. There is thus CARBON but 
| no simple correlation between activities observed This reaction should appear above 18.75+0.3 whi 
f at two different values of Emax, even though both Mev,® hence, it has not been produced with We 
i were produced by the same total irradiation. natural sources of gamma-radiation. reac 
i Activity, usually rather weak, was detected Cylinders of graphite 7.5 cm long, 2.5 cm in stat 
i with §-ray counters. The counter walls were diameter, and 2 mm thick were irradiated at 21 tain 
HI silvered glass 7 cm long, 16 mm in diameter, and Mev for 20 minutes and gave an initial activity nati 
| about 0.2 mm thick. They were filled with 9:1 of 200 counts per minute with a half-life of 2043 Thi: 
argon-alcohol to 9-cm pressure. The counter was_ minutes. At lower energies the yields became so dict 
coupled to a Neher-Harper circuit followed by a’ small that it was uncertain whether the activity 
scale-of-16 counter circuit. Background of about was actually 21-minute C"—almost any other 
35 min.—! was reduced to 15 min.~' by enclosure substance, such as adsorbed air, becomes active >) 
in a lead shield 10 cm thick. Still further reduc- at considerably lower energy. The total counts have 
tion of background is desirable. Samples were obtained at various energies, reduced to irradia- weal 
usually thick cylinders which fitted snugly over tions of 50R, are shown in Fig. 2. The cross- = 
the counters. hatched area is the statistical spread of counter (193) 
80 
‘F. W. Aston, Mass Spectra and Isotopes (Lon s, 
H arkas, an reutz, ys. Kev. ; A. 
For the first three of the following elements it Allison’ Phys. Rev. SS, 624 (1939); "L. P. Smith, Phys. i 
was possible to predict the threshold from mass __ Rev. 56, 548 (1939). 962 ( 
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background. The threshold of such a curve is 
taken as the energy below which points con- 
sistently fall within the background spread. The 
points below 19.4 Mev may or may not be 
genuine carbon activity. The threshold is proba- 
bly below 19.4 and above 18.7 Mev, in fair 
agreement with the predicted value—a high 
result naturally being observed with such a weak 
activity. More accurate work with carbon must 
await the production of higher intensity x-radia- 
tion than has yet been obtained with the betatron. 


NITROGEN 


As with carbon, yields from nitrogen” ® are ex- 
tremely weak, but there is less danger of con- 
tamination of samples, as the threshold is low. 
At 16-Mev peak energy the observed activity was 
only 600 counts per minute. Samples of NaNs, 
irradiated for 10-minute periods, were held for 
counting between concentric aluminum cylinders. 
The inner cylinder, 0.8 mm thick, was perforated 
with 600 evenly-spaced 2-mm holes, which were 
covered by 0.001” aluminum foil. This was 
placed over the counter. The annular section of 
NaN; thus employed was 4.5 mm thick. 

Figure 3, curve A, shows 30-minute counting 


totals corrected to an irradiation of 25R. Uncor- 


rected counts in the vicinity of threshold are 
plotted to larger scale in the inset, curve B. The 
difficulty of estimating a threshold is apparent, 
but no activity is observable below 11.1 Mev, 
while the sample is definitely active at 11.5 Mev. 
We have therefore placed the threshold for the 
reaction N'“(y,”)N™ at 11.1+0.5 Mev. The 
stated uncertainty includes estimated uncer- 
tainties both in energy settings and in determi- 
nation of the intercept of the excitation curve. 
This result is to be compared with the value pre- 
dicted in the introduction, above. 


OXYGEN 
The 2.1-minute positron activity of O'* should 


have a threshold of 15.6+0.6 Mev." It has been 


weakly produced with lithium gamma-rays.‘ *!° 


assy) Bothe and W. Gentner, Zeits. f. Physik 106, 236 

§O. Huber, P. Scherrer, O Lienhard, and H. Waffler, 
Helv. Phys. Acta XV, 312 (1942). 

*W. A. Fowler, L. A. Delsasso, and C. C. Lauritsen, 
Phys. Rev. 49, 561 (1936). , 

1°'W. Chang, M. Goldhaber, and R. Sagane, Nature 139, 
962 (1937). 


Preliminary measurements’ indicated a higher 
threshold, between 16.3 and 16.7 Mev. These 
measurements were repeated with improved 
technique, using water samples and an immersion 
counter of the same dimensions as the other 
counters used in this work. The volume of each 
water sample was 130 cm’, which surrounded the 
counter to a thickness of 1 cm. 

Samples were irradiated for 4 minutes to a 
total irradiation of 50R; counting was begun 30 
seconds after irradiation and continued for three 
half-lives to give optimum measure of activity. 
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Fic. 6. Excitation curve for Cu®. The energy values 
should be increased by about 0.1 Mev, as explained in the 
text. 


Figure 4 shows the resulting excitation curve up 
to 19.5 Mev, with points near threshold replotted 
to larger scale. The threshold is 16.3+0.4 Mev. 


IRON 


Iron gave a moderately strong 9-minute ac- 
tivity, which has been assigned to Fe™*.2* Bothe 
and Gentner*’ did not report this activity; 
Carlson and Henderson" failed to obtain it with 
higher intensity lithium gamma-rays. Huber 
et al.® recently obtained it with 17-Mev gamma- 


"uP, R. Carlson and J. E. Henderson, Phys. Rev. 58, 
193A (1940). 
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rays and reported a yield of 38 percent that ob- 
served with copper, based on the original isotope 
Fe*, It is surprising that it was completely 
missed in previous work—the yields obtained 
with the betatron were quite strong. Flat iron 
plates, 10X8X1 cm, were irradiated for 10 
minutes to 29R. For detection they were placed 
5 mm from an unshielded counter. The initial 
counting rate at 1 Mev above threshold was three 
times the background of 26 min.—! From Fig. 5 
the threshold is 14.2+0.4 Mev. 


COPPER 


The 10-minute Cu® is the strongest activity 
observed in a (y, ) reaction,” * and its threshold 
is quite sharp. Bothe and Gentner* measured the 
cross section for lithium gamma-rays as 5 X 10~** 


cm’, and with boron gamma-rays to be 2.3 times 


smaller. 


Copper tubes 7.5 cm long, 1 mm thick, and 2.5 - 


cm in diameter were used as samples. Irradiation 
was for 10-minute periods. The data are shown in 
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Fic. 7. Excitation curve for Zn®. 


Fig. 6, reduced to an irradiation of 15R. There is 
a slight systematic error in energy-values here 
because of phase-shift in the integrating unit, 
which is appreciable only when expansion occurs 
well below peak flux. The correction to be applied 


measurements of the (m,2m) threshold have 


in this case was found to be not more than 0,1 
Mev. The threshold has accordingly been taken 
as 10.9+0.3 Mev. This is considerably lower than 


indicated.” 
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Fic. 8. Excitation curve for total activity in selenium. The 
observed count is owing to two isomeric activities. 


ZINC 


Cylinders of zinc 6 cm long, 2.5 cm in diameter, 
of '3-mm wall-thickness were irradiated for 25 
minutes. Strong yields of Zn®™ were obtained.** 
The counts for irradiations of 50R are shown in 
Fig. 7. The threshold is 11.6+-0.4 Mev. 


SELENIUM 


X-rays induce two activities in this element, 
one of about 17-minute half-life (also reported as 
19-minute),? the other about an hour. Bothe and 
Gentner‘ obtained both by 1-hour irradiation of 
selenium with 17-Mev gamma-rays and also pro- 
duced them with slow neutrons. Langsdorf and 
Segré showed that they are isomers of either 
Se’® or Se*. These isomers were separated by 
chemical methods, and it was established that the 
57-minute selenium decays into the 17-minute 


#R, ne, Phys. Rev. 53, 492 (1938). 
13 A, pone and E. Segré, Phys. Rev. 57, 105 (1938). 
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isomer, emitting a 98-kev gamma-ray; the lower 
isomer decays by emission of a 1.5-Mev beta- 

icle. There is no evidence that the upper 
state decays directly by beta-emission. 

Plates of pure selenium, 7.5X5X1 cm, were 
irradiated for 15-minute periods, then placed 1 
mm from the counter surface. Total activity in 
35-minute counting periods, reduced to irradia- 
tion of 15R (computed on the basis of a 17- 
minute half-life), is plotted in Fig. 8. No activity 
results below 9.80.5 Mev (there was an uncer- 
tainty of as high as 3 percent in the energy 
settings in this experiment). If the short period 
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HOURS 
Fic. 9. Decay of selenium activity, Se7® * ™. 


is the lower isomer, this threshold should be that 
of the 17-minute period. Analysis of decay curves 
indicates that this is probably true. 

Designating the upper state by the subscript y 
and the lower by 8, the differential equations for 
production are 


—N,=—-—k,]I, (3) 
Ty 

Ne N. 

(4) 
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the lower isomer being produced by decay of the 
upper state as well as by the (vy, ) process. These 
equations could be solved for the initial popula- 
tions of the two states and for the ratio of initial 
activities which should be observed if the con- 
stants kg and k, were given. The beta-decay 
exhibits both half-lives. Since the beta-ray 
counter employed in the experiment had very 
low efficiency for gamma-rays, we need not con- 
sider the gamma-activity. With the simplifying 
assumption that the x-ray intensity is constant, 
it can be shown that 


ky 


where 


initial activity decaying with lifetime r¢ 


: initial activity decaying with lifetime 7, 
Note that in the absence of direct photoelectric 
production of the lower isomer (i.e., for kg=0), r 
would be negative, indicating that growth, rather 
than decay, should be initially observed for the 
activity with lifetime 

From decay curves (Fig. 9) it is evident that r 
is always positive and tends to increase at lower 
energies. Near threshold the curves become diffi- 
cult to analyze because of large relative statistical 
fluctuations, but it is still possible to estimate © 
upper and lower limits for the 57-minute activity. 
These estimates indicate that k,/kg decreases 
with decreasing energy, hence that the 57-minute 
threshold is probably higher than the 17-minute. 
This trend is shown in Fig. 10; as the energy is 
increased, k,/kg appears to attain a limiting value 
of about 0.6. In making any deductions from this 
trend, it should be borne in mind that the activity 
observed at a given peak energy is the integrated 
effect of all x-ray quanta above threshold energy. 
This treatment must be regarded as essentially 
qualitative as it was not possible to make a 
reliable analysis of decay curves near threshold, 
but it clearly indicates the trend. A chemical 
separation would be out of the question with the 
weak activities obtained with this reaction. 

MOLYBDENUM 


The 17-minute isotope of this element is the 
only one of convenient lifetime for the present 
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experiments.” Activities were moderately strong. 
A shorter period of about 1-minute half-life was 
observed above 14 Mev; it was allowed to decay 
for 5 minutes before counting was begun on the 
17-minute activity, and did not interfere with 
determination of the latter threshold. Analysis of 
the molybdenum samples revealed no impurity 
other than a slight trace of tungsten. No attempt 
has yet been made to identify this period. 
Samples of rolled sheet molybdenum, approxi- 
mately 0.8 mm thick, measuring 6 cm X7 cm, 
were irradiated for 10-minute periods. The counts 
plotted in Fig. 11 have been corrected to an 


‘irradiation of 30R. The threshold is 13.5+0.4 


Mev. 


SILVER 


Two activities are induced in silver by photo- 
disintegration.” 7 Ag! is positron-active with a 
24.5-minute half-life (the isomeric 8.2-day period 
is too long to be observed in the present study) ; 
Ag’* emits electrons with a half-life of 2.3 
minutes.’ To determine a threshold for each, it is 
necessary to separate their effects. This was done 
by analyzing the decay curves by somewhat the 
same procedure as that used for selenium. 


Fic. 10. Ratio of rate of production for upper isomer to that for lower isomer 
of Se7? oF 81, 


MEV 


A silver plate 7 X6.5 cm, 1 mm thick, was rolled 
into a cylinder which fit snugly over the counter. 
This was allowed to decay completely between 
successive runs. A series of 25-minute runs gave 
the decay curves shown in Fig. 12. From this and 
from Fig. 13A it is evident that the short period 
predominates at lower energies. Below 10 Mev it 
is difficult to estimate the longer activity, but an 
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Fic. 11. Excitation curve for 17-min. Mo activity. 
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upper limit for it can be set as in the case of 
selenium. Part B of Fig. 13 gives totals in a series 
of 3-minute runs to irradiations of 3.2R, indi- 
cating a threshold of 9.30.5 Mev for the 2.3- 
minute period. 

Since each isotope is produced independently 
of the other, growth under irradiation follows 
Eq. (1). The ratio of the constants kio¢ and kos is 


ki06 

where the R’s are the weighted irradiations for 
each isotope and 


initial 24.5-minute activity 


r= . 
initial 2.3-minute activity 


Were each isotope irradiated to saturation at the 


TaBLE I. Threshold energy values for nuclear photo-effects. 


Target Active Half-life Threshold 
element isotope (minutes) (Mev) 
6 Carbon cu 21 18.7 to 19.4 
7 Nitrogen N® 10 11.1+0.5 
8 Oxygen or 2.1 16.3+0.4 
26 Iron 9 14.2+0.4 
29 Copper Cu® 10.3 10.9+0.3 
30 Zinc Zn® 39 11.6+0.4 
34. Selenium Se? or 17 9.8+0.5 
Se79 or 8 57 higher 
42 Molybdenum Moors 17 13.5+0.4 
47 Silver Agi 2.3 9.3+0.5 
Agi 24.5 above 9.5 


same x-ray intensity, k106/k10s would be the ratio 
of activities observed. Values of this ratio com- 
puted from decay-curve data are plotted in 
Fig. 14. At and below 9.5 Mev the 24.5-minute 
activity is not produced. It is therefore concluded 
that the 2.3-minute threshold lies below the 
24.5-minute threshold by at least 0.25 Mev. 


CONCLUSIONS 


Thresholds which. have been measured so far 
are tabulated in Table I. Probable errors listed 
include uncertainties both in energy settings and 
in determination of threshold from inspection of 
the excitation curve. The chief difficulty in the 
latter arises from statistical fluctuation in counter 
background. A threshold determination is limited 
by the sensitivity of the means of detecting the 
effect. Since the effects are comparable with 
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Fic. 12. Decay of silver activity at various 
irradiation energies. 
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Fic. 13. (A) Initial activities of Ag! and Ag’ os. 
irradiation energy. (B) Excitation curve for silver with 
short irradiation period. 
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Fic. 14. Ratio of production rates for Ag'** and Ag’. 


background fluctuations in several cases, one 
cannot be completely satisfied that the true 
threshold may not lie considerably below the 
point where the effect first becomes detectable. 
The first three elements, yielding weak activities, 
gave thresholds consistently higher than pre- 
dicted values. This may indicate that the true 
thresholds are below those read from the curves. 
It should be pointed out, however, that emission 
of gamma-rays after the beta-decay would make 
the calculated threshold seem too low. This might 
be the case for oxygen. 

In other elements, the relatively better data 
warrant confidence that the true thresholds do 
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Fic. 15. Decay curves for photo- 
fission of uranium, observed by 
collection of fission recoils. 
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NET COUNTS 


not differ from those listed by more than the 
errors given. 

No interpretation of the shapes of the excita- 
tion curves has been attempted. The spectra] 


* distribution of the radiation employed, which 


arises from a target of intermediate, variable 
thickness, must be known before actual cross 
sections and true excitation curves can be found. 
There is a general similarity of shape in all the 
curves. 

All the measured thresholds are appreciably 
above the mean binding energy per particle, 
which is under 8 Mev in all cases studied. The 
trend is toward lower values in the heavier 
elements, but the list is not yet large enough to 
warrant generalizations from the data. An ex- 
tension of the list to heavier elements is con- 
templated, together with a thorough rechecking 
of the present list. Present conditions make it 
unlikely that this work can be completed in the 
near future, however. 


PHOTO-FISSION 


Photo-fission of uranium was also investigated 
by collecting fission recoils on a paper catcher and 
measuring the collected activity. Uranium oxide 
was coated on a cylindrical form. A cylinder of 
paper was slipped over this sample and held by 


MEV 


Fic. 16. Excitation curve for pho- 
to-fission of uranium from measured 
activity of collected recoils. 
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spacers 1 mm from the uranium. This was 
irradiated for 10-minute intervals at various 
energies, and the activity was measured with a 
counter. Blank tests on ordinary bond paper 
showed no activity below the oxygen threshold 
at 16.3 Mev. 

Decay of the collected fission fragments was 
followed for 45 minutes after runs at 13 and 15 
Mev. The points shown in Fig. 15 were fitted 
by a double exponential of half-lives 2.2 and 18 
minutes, though the full decay scheme is much 
more complex. The activity is not strong enough 
to permit resolving the various activities present. 

Net 30-minute counting totals (background of 
500 subtracted out), corrected to an x-ray in- 


tensity of 1R per minute at one meter (or an 
irradiation of 10R, if no allowance be made for 
decay), are plotted in Fig. 16. Below 8 Mev the 
collected activity becomes too weak to follow. 
The threshold is apparently below 7 Mev, 
however. 

Attempts to observe photo-fission in lead by 
this method were unsuccessful at energies up to 
16 Mev, where the collector became active. 
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The non-relativistic wave function (in spherical polar coordinates) for a charged particle in a 


Coulomb field is expressed in a form suitable for problems in which the particle has a small 
positive energy. This formulation amounts to expanding the radial part of the wave function in 
powers of the energy E and is achieved by simple algebraic manipulation of power series and 
recurrence formulas. The coefficients of the expansion are functions of the radial coordinate and 


FUNDAMENTAL EQUATIONS 
PARTICLE of mass m and charge! Z’e is 
moving in the Coulomb field of a nucleus 

with charge Ze. If the energy of the particle is 
E(>0) and if it is described in terms of spherical 
polar coordinates, 7, 6, @ (origin at Ze), then a 
suitable wave function for the particle is? 


Viim(, 6, ¢) = (cos 6) Lx, i(r), 
r (1+ 1 +i) 


X (2kr)'e** 21+2, -2itr) 


B=ZZ'me*/h*?; k?=2mE/h* 


* Publication niet the Ernest Kempton Adams 
Fund for Physical R of Columbia University. 

1¢@ is positive throughout, so that, e.g., Z’=—1 for an 
electron or +2 for an alpha-particle. 

2N. F. Mott and H. S. W. Massey, 7 of Atomic 


Liil(r) 


/v(2t+-2)] 
1 


Collisions (Oxford University Press, 1933), p. 39. 


are identified with Bessel functions of integral order. 


in which the constant factors in L;, ,(r) have been 
chosen so that asymptotically 


sin 
arg In | 


The confluent hypergeometric function ,F; is 
defined? as 


a@ a(a+1) 2? 
1F (a, 6, 2) 


To simplify the analysis consider, instead of 
Ly, :(r), the function [essentially :(r)] 


(2) 


+E. T. Whittaker and G. N. Watson, Modern Analysis 
(Cambridge University Press, 1927), p. 338. 


which is real and satisfies the differential equation 

ru’ +-(21+-2)u’ + —28)u=0. The power series 
solution (regular at r=0) of this differential 
equation is: 


Ux, = b,(k?, 


bo=1 


The condition b)=1 follows from the definitions 
(2) and (3). 6, has been written as a function of k? 
since, as is apparent from the recurrence formula, 
only the square of & is involved. 

In calculating numerically the radial part of 
v(r, 0, @) one can use Eq. (4) with a reasonable 
amount of labor (provided r is not too large). 
When &? is small, however, it is possible to 
simplify the calculation and to derive an expres- 
sion for u of the form: 


p=0 


The purpose of this paper is to develop formulas 


for the functions u»,,(r), using the properties of 
the expansion coefficients },(k?, 2). 


COEFFICIENTS b, EXPANDED IN POWERS OF k? 
From Eq. (4) it follows that 


+h [(1—n)/(n +1420) 


If b,-1 is replaced by the value of },_1 obtained 
from the same recurrence formula, that is, 


=(28/(n +21) Jona +21) ]on-s, 
and terms in k? are collected, the result is: 


28 (2— 
+ 2141) (n+21) 


Now repeat this procedure, substituting for the 


first b,_2 on the right the expression D,_2= - 
+k?---b,_4. Again collect terms in k? and obtain 
nally repeat this »—3 times until the first term 
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has the factor bo (which = 1). The general theorem 
thus derived is: 


1 ED! 
(m+ 


From this formula it is possible to develop an 
expansion of b, in powers of k?. The coefficient of 


the zeroth power of &? is just the first term on the 
right side of the equation. The coefficient of & is 


1). (6) 


* obtained by substituting the zeroth power coeffi- 


cient into the finite sum. Coefficients of higher 
powers are obtained in a similar manner. In each 
case the only difficulty is the evaluation of the 
finite sum over v. In terms of the finite sums the 
result of this process is an expansion of the }, in 
powers of k?: 


b,(k?, 1) = > b,*(0, 2) 


b,*(0, L) = 


7 
X 


. 


Now substitute this expansion into Eq. (4) and 
compare the result with Eq. (5). It is immediately 


apparent that the functions Uo, i(r) are just 
uo, 1(r) = > b,*(0, 1)r"/n! 
=(=1)9(28)-* [(21+1)!/ 
n=0 


Sa” 
(n+21+1) (8) 


The power series of (8) is to be compared to the 
expansion of the Bessel function A,[2(—26r)*] 
defined‘ as 


26r)-?!2J 
[p!/(n+P) (9) 
The comparison ‘shows that if the S,“ are not 


*E. Jahnke and F. Emde, Tables of Functions (B. G. 
Teubner, 1938), p. 128. 
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too complicated, it should be possible to reduce 
the uo, :(r) to a finite sum of A,'s. 

Inspection of the definition of S,*, Eq. (7), 
shows that the S,“ are polynomials in m of degree 
3u. In evaluating S,* it is convenient to express 
these polynomials as a sum of terms each with the 


general form 
X 


where A, \’>0 (and A+A’+2<3y since the 
‘highest power of m must be n*). This form is 
chosen so that when it is substituted in place of 
S,* in Eq. (8), the resulting series in m can be 
compared to Eq. (9) and immediately identified 
as 


EVALUATION OF THE FINITE SUMS 


Instead of the definition of S,*, Eq. (7), it is 
more convenient to consider the difference equa- 


tion satisfied by S,*: 


By solving this equation [or Eq. (7)] for the 
special cases x= 1 and u=2, one is led to surmise 
that the general solution which has the desired 
form is: 


It is apparent that this is of the appropriate 
degree (3u) in m, but that it has the necessary 
number of arbitrary constants [to be determined 
by Eq. (10)] is not immediately apparent. 

The validity of (11) is readily established by 
substituting it into Eq. (10), expanding both 
sides in powers of m, and comparing coefficients. 
The result is a series of u+1 equations defining 
each a,, in terms of a,,-1 and a,-1,, for A<». 
Starting with y=0 and then considering »=1, 2, 


+++w, each successive equation for the a,, is a 
linear equation in one unknown (assuming the 
a@,-1,, to be previously determined) : 


Ay, 0 = Ayp-1,0= 1, 
1 
(3u—») 


ay =0 for 
with 


(n+21+1)! 


wan, 


The @’s and y's can be expressed explicitly in 
terms of yu, A, J, and substituted in the a,, 
recurrence formula above. Then taking y=1, the 
sum over \ reduces to one term, \=0, and since 
@,-1,0= 1, a simple relation between a,; and 
a@,~1,1 is obtained. This can be solved in terms of 
and J, yielding In princi- 
ple it is possible to continue this process finding 
Qy2, etc., as explicit functions of and /. 
However, the algebra becomes tedious mainly 
through the occurrence of polynomials in /, so 
that in practice it is easier to calculate an ex- 
tensive table of a,, when some particular nu- 
merical value of / is used. , 
The values of a,, for 1.=0, 1, 2, 3 are: 


w=0: an=1 

a@21 = (2/5)(SJ—1), 
ae2=(1/5)(SI—1) (12b) 

a3 = (3/5)(Sl—2), 

a32 = (3/35)(35/?—211+-4), 


= (1/35)(351?—211+-4). 


ayo=1, 


ao=1, 


azo=1, 


it of 
1 the 
k? is 
effi- 
gher 
the 
b, in 
= > Yur . 
itely 
(8) 
v=0 
> the 
pr)*} 
(9) 
» not 


- EXPANSION OF ¢ IN POWERS OF k? 

Substitution of (11) into Eq. (8) yields a gen- 

eral formula for the functions wo,:in terms of sums® 
of u Bessel functions of integral order: 


(—1)4(21+1)! 


(21-+1+p+»)! 
X (13) 
with 


26r)*] = 
X (—2Br)— 12 2(— 267) 4]. 


For repulsive fields 8 is positive so the Bessel 
functions have an imaginary argument; in (13) A, 
is to be taken as p!(28r)—?/*4-?J,[ 27(26r)*]. B is 
negative for attractive fields so the Bessel func- 
tions have a real argument. (For 6>0 the Bessel 
functions are monotonically increasing functions 
and for 8<0 they are oscillatory functions with 
decreasing amplitude.) 

Thus the wave function [asymptotically 
Ne™¢P;'™! (cos 6) sin (kr---)/(kr)] for a particle 
with charge Z’e in the Coulomb field of Ze is: 


Veim(r, 0, 6) = Ne™*P,|™! (cos (2) 


[(—1)*/(3"u!) 


x|r (14144 


xr 26r)*] 
mmo 
As long as & is finite the radial part of Yrim is 
asymptotically (for 1r>>1/k)~sin (kr---)/(kr). 
For k=0 the radial function behaves at large 
distances like A214:[2(—26r)#]. 
In (14) the factor (—1)*(kr)** appears. This 
means that the expansion is an alternating series 


(14) 


5 The number of terms in each of the sums ma may De 


reduced to two terms by repeated application 
recurrence formula for A, or J»: 


Ap+i(x) = [4p(p+1)/x*][Ap(x) —Ap-1(x) 
= (2p/x) —Jp—-1(x). 


14 J. G. BECKERLEY 


in powers of (kr)?; consequently the convergence 
will be rapid only when kr<1. For kr>1 the 
asymptotic sin (kr---)/(kr) is valid, while for 
kr~1 (and 8/k>1) it is necessary to derive an 
asymptotic formula for the 1F; of Eq. (3) valid 
when the first parameter is large and imaginary, 


CONCLUSION 


As formula (14) is most useful when kr is smal] 
enough to make only the first few terms im- 
portant, these terms are explicitly presented 
here: 


Viim(7, 0, 6) = (cos (2) 


(kr)? 


(kr)* 2 gy 


Asi+s 
(i+4)(2145)) | (15) 


where A,=A,[2(—26r)*]. 


The above expansion® of 8, 6) in powers 
of (tkr)* is particularly useful in calculating 
matrix integrals involving slow positive energy 
electrons in a Coulomb field provided the inter- 
action is such that the most important region of 
integration is near r=0. In this case the con- 
vergence of the series, Eq. (14), is so rapid that 
only the first few terms are needed in the 
calculation. 

A special case of the theorem, Eq. (15), has 
been used by the author in his Ph.D. dissertation 
(unpublished) submitted to Stanford University, 
California. 

* This same expansion can no doubt be derived by 


transformation of the integral ntation for the func- 
tion u:i(r) instead of by the algebraic methods of the 


present derivation. 
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Effect of Nuclear Electric Quadrupole Moment on the Energy Levels of a 
Diatomic Molecule in a Magnetic Field*t 


Part I. Heteronuclear Molecules 


BERNARD T. FELD AND Wix.is E. Lams, Jr. 
Pupin Physics Laboratories, Columbia University, New York, New York 
(Received September 28, 1944) 


The energy levels of a diatomic molecule—in which one nucleus has electric quadrupole 
moment and the second zero spin—are derived in ‘‘weak’’ and “strong” magnetic fields, and for 
intermediate fields in the special case of nuclear spin 1. The calculations are made for large 
molecular rotational quantum number. The transitions which would be induced in a molecular 
beam resonance experiment are discussed: Formulae are given for the energy differences in- © 
volved in these transitions, and the shape of the molecular beam resonance spectrum, corre- 
sponding to the various transitions, is considered in detail. It is shown that the presence of 
nuclear electric quadrupole moment introduces subsidiary minima into the spectrum; from the 
energy differences between these minima, one can derive the value of the constant e*gQ (e is the 
electronic charge, Q is the nuclear electric quadrupole moment, and q is a constant, depending on 
the charge distribution of the rest of the molecule, which must be independently evaluated in 
each case). The possibility of experimental application of these results is discussed. From _ 
heretofore unpublished data, we obtain a value for egQ of Li? in LiF of 1.00.5 X 10 erg. This 
may be compared with 0.3 10~*! erg for D in HD and Dz. 


I. INTRODUCTION 


HE existence of a nuclear electric quadrupole 
moment was first inferred from spectro- 
scopic measurements of hyperfine structure. 
Thus, the. departure of the hyperfine energy 
differences of europium! ? from the interval rule 
was explained on this basis. The molecular beam 
“zero-moment”’ technique was first applied to the 
measurement of the nuclear quadrupole moments 
of In"™5* and Ga®*”4 at the suggestion of I. I. 
Rabi. One of the early applications of the ‘‘mo- 
lecular beam resonance method” was the determi- 
nation by Kellogg, Rabi, Ramsey, and Zacharias® 
of the quadrupole moment of the deuteron from 
measurements of the hyperfine structure of HD 
and Dz in a magnetic field. The molecular beam 
resonance method has been discussed in a number 


* Submitted by Bernard T. Feld in partial fulfillment of 
the requirements for the degree of Doctor of Philosophy in 
the Faculty of Pure Science, Columbia University, New 
York, New York. 

t Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 
(sas) Schiiler and T. Schmidt, Zeits. f. Physik 94, 457 

2H. B. G. Casimir, Physica II, 719 (1935). 

*D. R. Hamilton, Phys. Rev. 56, 30 (1939). 

4N. A. Renzetti, iy Rev. 57, 753 (1940). 

a M. P. Kellogg, I. I. Rabi, N. F. Ramsey, Jr., and 
J. R. Zacharias, Phys. Rev. 57, 677 (1940). 
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of papers.*7 However, a short review of the 
highlights will probably not be amiss. 

A beam of well-collimated molécules, at low 
pressure, is passed successively through three 
magnetic fields. The first and third are inhomo- 
geneous fields, equal in magnitude and direction 
but with gradients oppositely directed. Between 
these, there is a homogeneous magnetic field. 

Because of its magnetic moment, the molecule 
is deflected in the inhomogeneous field. If no 
change in the magnetic moment of the molecule 
takes place, while it goes through the homogene- 
ous field, the deflection in the two inhomogeneous 
fields will be equal and opposite; as a result, the 
fields will have no effect on the intensity of the 
beam, as measured by a suitable detector placed 
at the end of its path. However, should the 
orientation or magnitude of the magnetic mo- 
ment change between fields, the deflections will 
not compensate ; there will then be a decrease in 
beam intensity. 

Changes in ‘orientation are brought about by 
the application of a small, oscillating magnetic 
field perpendicular to and in the region of the 


* 1. I. Rabi, S. Millman, P. Kusch, and J. R. Zacharias, 


Phys. Rev. 55, 526 (1939). ° 
aK S. Millman, and I. I. Rabi, Phys. Rev. 57, 765 
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homogeneous field. In the homogeneous field, the 
quantum mechanical system (molecule) will tend 
to orient itself in such a way that its spin satisfies 
the stationary state quantum conditions for a 
magnetic moment in a magnetic field; the com- 
ponent of the angular momentum in the direction 
of the field must be an integral or half-integral 
multiple of /27. 

In general, corresponding to different orienta- 
tions of the magnetic moment in the magnetic 
field, the system has different energies. When h 
times the frequency »v of the oscillating field is 
very close or (better) equal to the difference in 
energy W between two states of the system and 
when the quantum mechanical conditions allow 
a transition between these states, there will occur 
a reorientation within the system, leading to a 
change of the resultant magnetic moment of the 
molecule. Thus, for values of the frequency of 


Fic. 1. Li? resonance in LiF. The ordinates are arbitrary, 
while the abscissa numbers are galvanometer readings, 
a ge ogee to the strength of the homogeneous magnetic 
field. The perturbing oscillating field ern sed was held 
constant. This is one of a number of heretofore unpublished 
curves, for which we are indebted to Drs. Kusch and 


oscillating field such that 
hv=W, I (1) 


there will be a decrease of beam intensity at the 
detector. Measurement of the resonance fre- 


B. T. FELD AND W. E. LAMB, JR. 


quencies and corresponding values of the homo. 
geneous field at the resonance minima may lead 
to valuable information concerning such prop. 
erties as gyromagnetic ratio, hyperfine’ splitting, 
electric quadrupole moment, etc. 

An important application of this method has 
been the measurement of nuclear magnetic mo. 
ments in molecules having zero electronic angular 
momentum.*® In one of these investigations,* 
it was found that the resonance minima for Na® 
and K*® (in the molecules Naz and Ke) are ac. 
companied by subsidiary minima. Furthermore, 
these subsidiary minima (for a particular nucleus) 
appear always to be separated by a constant 
energy difference, independent of the value of 
frequency and magnetic field at which the reso- 
nance occurs. This seems to imply the existence 
of a constant interaction between the nucleus and 
the rest of the molecule. In addition, the investi. 
gators® have shown that the effect is at least of an 
order of magnitude too great to be explained as 
magnetic dipole interaction between the two 
nuclei or between the nuclei and the rotational 
magnetic moment of the molecule. To quote the 
authors: “It appears likely to us that these large 
effects may be due, in part, to the presence of a 
quadrupole electrical moment in the nucleus 
which interacts electrostatically with the rest of 
the molecule.” 

Very similar effects have been observed by 
Kusch and Millman" in a number of heteronuclear 
molecules. In Fig. 1 is plotted one of their 
upublished curves for a resonance involving a 
reorientation of the nucleus of Li’ in LiF. The 
subsidiary unresolved minima have been ob- 
served always to appear in the curve under a 
variety of conditions. The figure is an example of 
resonance in a molecule of the type to which the 
calculations of this paper pertain. 

We have, in fact, calculated the energy levels 
and spectrum for heteronuclear molecules in a 
magnetic field when one member has electric 
quadrupole moment. The calculation has been 


-§P. Kusch, S. Millman, and I. I. Rabi, Phys. Rev. $5, 


1176 (1939). 
®S. Millman, P. Kusch, and I. I. Rabi, Phys. Rev. 56, 


165 (1939). 
10 J. M. B. Kellogg, I. I. Rabi, N. F. Ramsey, Jr., and 


J. R. Zacharias, Phys. Rev. 56, 728 (1939). 

1 We are greatly indebted to Drs. Millman and Kusch 
for having a omen | us access to a large amount of unpub 
lished data. 
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ENERGY LEVELS OF A 


made for fields ranging from ‘‘weak”’ to ‘‘strong,” 
- for nuclear spins from 1 to 9/2, and large rota- 
tional quantum number J. A ‘‘weak”’ field is one 
in which the energy of interaction between the 
field and the magnetic-moment of the molecule is 


Fic. 2. Diatomic molecule with rotational spin J, and one 
nucleus having spin I and electric quadrupole moment. The 
second nucleus has zero spin. The molecule is in a magnetic 
field H. 


small compared to the electric quadrupole inter- 
action energy of the molecule; thus, the magnetic 
interaction terms can be treated as perturbations 
on the energy of the system. In “‘strong’’ fields, 
the electric quadrupole interaction is small com- 
pared to the magnetic interaction. A calculation 
of this type has already been made for the mole- 
cule HD,‘ and differs from our calculation in that 
the rotational quantum number J is quite small 
in HD. Although the condition of large J intro- 
duces many quantum states into the problem, it 
results in a simplification insofar as J can now be 
treated as a classical angular momentum vector, 
and the magnitude of the nuclear spin J can be 
neglected as compared to that of J. 

For the purpose of simplification, we have 
made the assumption that the spin of the second 
nucleus comprising the molecule is zero. As long 
as the electric quadrupole interaction is large 
compared to the magnetic interaction between 
nuclei, and between nucleus and the field caused 
by molecular rotation,” this assumption may be 

2 If the nuclear magnetic moment is wu» and the distance 
between nuclei R, then the magnitudes of both the nuclear 
ppt and spin-orbit energies are of order u,?/R*. For 
a nuclear magnetic moment of one nuclear magneton and 
R of 10-* cm, this energy turns out to be of order 10-* AB, 


The quadrupole interaction has order of magnitude 1 
where ¢ is the electronic charge, g is a factor depending on 


DIATOMIC MOLECULE 17 
made without detracting from the generality of 
our treatment. For it is evident that, if the 
afore-mentioned considerations hold, the energy 
of the molecule (because of quadrupole inter- 
action) depends only on the orientation of the 
quadrupole nucleus with respect to the rotational 
angular momentum. The energy of the molecule 


’ inamagnetic field consists then of two completely 


separate parts, pertaining to the two nuclei; so 
that the effects of the quadrupole moment of one 
nucleus are independent of the state of the second, 
and energy changes caused by the reorientation 
of one nucleus are not affected by the state of the 
other. 

A diatomic molecule of this type may be 
pictured diagrammatically as in Fig. 2. The 
nuclei are rotating about the direction J; de; 
represents a small element of electronic charge. 
The rest of the charge, which is distributed in a . 
suitable fashion throughout the molecule, is not 
shown. The externally applied magnetic field is 
represented by the vector H. ; 

The part of the Hamiltonian for this molecule 
in an external magnetic field H, which is of 
interest in our problem, may be written as 


H = -I+pogsH-J+35(Q), (2) 


where yo is the nuclear magneton, g; and gy are 
the Landé g factors for the nuclear and rotational 
spin, respectively, and 3C(Q) is the Hamiltonian 
for the interaction of the electric quadrupole 
moment of the nucleus with the rest of the charges 
in the molecule. 


Il. STRONG MAGNETIC FIELDS; 
PASCHEN-BACK EFFECT 


(1) Semiclassical Treatment 


In sufficiently strong magnetic field, nuclear 
spin and molecular rotation may be considered 
as uncoupled. One can adopt a semiclassical 
approach. Allowing the subscript 7 to represent 
the nucleus with spin and quadrupole moment, 
and j the remainder of the molecule, we write 
as the quadrupole term in the energy expansion 
the electronic charge distribution, of order 1/R*, and Q is 
the value of the electric quadrupole moment. For values of 
Q of order 10-** cm? (magnitude of the Q’s observed spectro- 
scopically) e*gQ ~10-* erg. In the case of HD, the quad- 
rupole moment of the deuteron having been determined to 


be = 10-*? ~ 10-* and the magnetic effects could 
not be neglected. 
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for our system 
eQN de;(3 cos? 1) 
E(Q)= 
J II (1) 


In this expression, the quadrupole moment Q is 
defined by the integral over the nucleus 


eQ= f deg 7(3 cos* 6r;—1), II (2) 
i 


and N is a normalization factor, to be discussed. 
The angle 6;; is measured between the direction 
of nuclear spin and an element of nuclear 
charge de;. Likewise, 6;; is the angle between 
the nuclear spin direction and the element of 
charge de; (Fig. 2). 

Introducing the direction of molecular rotation 
J, we apply the addition theorem of spherical 
harmonics." Thus, we may write 


(3 cos? 1) = }(3 cos? 67;—1)(3 cos? 6r7—1) 
+terms with exp (ig,;) exp (igrz). II (3) 


Performing the integration II (1) over the mole- 
cule, we find the terms exp (ig,;) average out 


to zero, giving 
e*gQN 
E(Q) cos? 6rz3—1), II (4) 


where 


de;(3 cos? 1) 


Taking advantage of the fact that I and J, 
for strong fields, may be considered to precess 
independently about the direction of the field H 
(Fig. 2), we write 


(3 cos* 67 1) = $(3 cos? 6;—1)(3 cos? 67 —1) 
+terms with exp exp (ig). II (6) 
Averaging over all directions about H, we obtain 


cos? cos? @;—1). II (7) 


Since J is large, the molecular rotation may be 
considered from a classical point of view. Hence, 
cos 0;=m,/J ranges continuously from —1 
to +1. 


18 E. U. Condon and G. H. Shortley, Theory of Atomic 


re (Cambridge University Press, 1935), p. 53 (hence- 
referred to as Condon and Shortley). 
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The nuclear spin J, being small, must be 
treated quantum mechanically. Thus, only cer. 
tain orientations of I with respect to H, charac. 
terized by integral or half-integral m; (cos 6, 
=m,/|I|) are permitted. In addition, there ig 
the relationship P= J(J+1). 

In the definition of quadrupole moment by 
Eqs. II (1) and II (2), there is introduced a 
normalization factor NV. This factor is evaluated 
by incorporating into the definition of quadrupole 
moment the condition that when m;=TI, the 
factor eQN}(3 cos? 6;—1) becomes equal to 
From this and the relationships introduced jn 
the last paragraph, we have 


21I+1) 


Thus, we arrive at an expression for the quadru- 
pole term in the electrostatic energy of the 


molecule: 


—I(I+1)] 
1(2I—1) 


The quadrupole energy is seen to be proportional 
to the product gQ. Thus, knowledge of the 
quadrupole energy alone will not suffice to 
determine the quadrupole moment Q. 

The quantity g depends on the electronic con- 
figuration of the molecule and can be evaluated 
if this distribution is known.'® Although no 
attempt shall be made in this work to calculate g, 
it should be stressed that before the results of 
our calculations can be utilized for the deter- 
mination of nuclear electric quadrupole moments, 
an independent evaluation of g is required. 


(2) Quantum Mechanical Treatment 
(a) Energy Values 


In strong magnetic field, the Hamiltonian for 
our system may be written as the sum of two 


terms. The first 


Ho= IT (10) 


4H. B.G. Casimir, “On the interaction between atomic 
nuclei and electrons,” prize essay published by Teyler'’s 
Tweede Genootschap (1936). Also in the 
Archives du Musee Teyler [111], VIII, 201 (1936) (hence 


forth referred to as Casimir). 
16 A. Nordsieck, Phys. Rev. 58, 310 (1940); E. Teller, 


Phys. Rev. 57, 556A (1940). 
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in which the symbols are used as defined in the 
preceding section, expresses the interaction be- 
tween the magnetic moments associated with I 
and J and the external magnetic field. This 
term is diagonal in a representation characterized 
by the quantum numbers m;, my. 

The perturbation term of the Hamiltonian 
results from the electric quadrupole interaction 
between nucleus and the rest of the molecule, 


(ms, mr| F| mz, mr) 
(ms, F|ms—1, mr+1)= 


DIATOMIC MOLECULE 


and is given by the expression“ 
e*gQ 
D 
—I(I+1)J(J+1)}. (11) 


Referring to the operator in brackets as F, 
we have for its matrix elements in the mr, | my 
representation :5 


—I(I+1) 


(mys, mz| F|ms+1, mr—1) = }(2my+1)(2m1—1) — mz) (J+ms+1) (I —mr+1)}', 


(mz, mr| F|ms—2, 


—ms+2) 


(II 12) 


(mys, mr| F|ms+2, m;—2)= 


$L(J 


X }}. 


In the case <b, first-order perturbation 
theory gives for the quadrupole energy, 


e*qQ 
(0) = 
4J(2J—1)I(2IT—1) 


X J (J +1) 
which, for large J, reduces to 
—I(I+1)] 

()(0) = 
1(2I-1) 


This is, of course, identical with the result ob- 
tained in the previous section. 

From the value of the Hamiltonian and the 
above matrix elements, it is possible to calculate, 
by the use of perturbation theory,’* the energy - 
level system of the molecule for different values 
of I. The case J=} is of no interest to us, be- 
cause the charge density of the nucleus in this 
state is spherically symmetrical and, hence, 
Q=0. 

For the systems under consideration, there are 
(2J+1)(2I+1) different states. However, since 
J is large, we shall calculate only the 27+1 
states of different my, retaining m,/J (which can 
be considered as varying continuously) as a 
parameter. 

The energy values, including second-order 


T(I+1)], II (13) 


II (14) 


perturbation terms, for the case J=1, are: 


6 Reference 13, p. 30. 


(The expressions for J2 $ have the same form, 
although differing in detail.) 


e*qQ 
E(m;=1) = + — 1) 


(322-+1)(1—z") 
— wogsH 


e*gQ 


E(m;= —1) = —wogrH +mogsHmys 


1) 


(322+ 1)(1—2*) 
— wogsH 


* In the above expressions, we have substituted the 


continuously varying parameter z for m;/J. 
(6) The Spectrum 


Transitions between states of our system are 
caused by a small perturbation in the form of a 
weak, oscillating magnetic field H’, whose con- © 
tribution to the Hamiltonian (too small to 
effect the energies of the system) may be ex- 

17 = = & 


all further computations write gy=eg;, where « is a 
number. 
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field equal to the energy difference divided by h. 
The probability of a transition is, in general, 
a function of the matrix elements between states, 
as well as of the velocity distribution in the 
molecular beam and the length of the oscillating 
field. Assuming a monoenergetic beam subject 
to a small perturbation (such that only a small 
fraction of the molecules capable of undergoing 
a transition are effected) we may take the 
transition probability P to be proportional to 
the square of the corresponding matrix element.”° 


#8 Reference 13, pp. 47-49. 
19The resonance process has been discussed in the 
following references: I. 1. Rabi, Phys. Rev. 51, 652 (1937) ; 
i: Schwinger, Phys. Rev. 51, 648 (1937); F. Block and A. 
iegert, Phys. Rev. 57, 522 (1940); A. F. Stevenson, Phys. 
Rev. <8, 1061 (1940); H. C. Torrey, Phys. Rev. 59, 293 
(1941). 
20 The transition probability between two states p and q, 
if we assume no interference from other transitions, is 


beat 


where byg is the matrix element of the perturbation 
(o15'|g) divided by h, ¢ is the time spent in the perturbing 
field, and yp, is the frequency at which the resonance 
takes place (hvpg=Wp5,). At resonance, and for small 


Py 


2 
byg, Ppq= 2°l*bpq. The natural width of a line is most simply 
obtained from the uncertainty principle AEAt~h. For 
molecules (assumed molecular wt.~50 and temp.>1000° 
Abs.) in an oscillating field of length 5 cm. AE~5X 107% 
erg. Since this is considerably less than the energies in- 
volved in the Dye og interaction (see footnote 12) we 
shall henceforth neglect natural widths in our discussion 


of line shapes 


20 B. T. FELD AND W. E. LAMB, JR. 
pressed as 
In general, 
130F (o,-1) 
H’ = (H,/i+H,/j+H,'k) cos (wt-+8). II (17) 
In the approximation for which: m; and my are 
the appropriate quantum numbers, I and J may 
be pictured as precessing independently about ae P 
the direction of the homogeneous magnetic (0,-1) x#.0) 
field. The matrix elements for such angular 
momenta are given in Condon and Shortley :!8 (0) Ee 
90) 
(m|I|m)=mk, 
2 
In effect, we are dealing with a resonance a ** x= ant 
process'® for which transitions between states 
occur at values of the frequency of oscillating. : : 
° 22 2.4 x 28 t'0 


Fic. 3. Energy of transition for J=1 in strong magnetic 
fields as a function of z=m,/J, for different values of the 
perturbation x. The two sets of curves are for the transition 
m,=0—>m;= — 1 and m;= 1m; =0, respectively. (Am; =0 
in both cases.) 


In considering a transition, then, two factors 
must be taken into account: (1) the energy 
differences W are in general functions of the 
parameter z=m,/J[W=f(z) ]. This parameter 
varies continuously from —1 to +1, each value 
having equal probability of occurrence. For 
an interval AW there corresponds an interval 
Az=AW-dz/dW. Since the number of values of z 
in the interval AW, and hence the number of 
transitions within this interval, is proportional 
to Az, we may define a “density of transitions” D, 
given by 

D(z) =Az/AW=1/f"(z). II (19) 


The reciprocal of the slope of the W vs. z curve 
(Fig. 3) will thus be proportional to the number 
of states involved in transitions having a given 
energy difference; hence, the depth of the reso- 
nance minimum will vary with this quantity. 
(2) The product of D(z) by the transition proba- 


’ bility P(z) is a measure of the relative strength 


ofthe transition, as a function of the parameter z. 
D(z)P(z) will henceforth be called the “shape 
function”’ for a transition. 


moe 


| | 


oscillating frequency is varied through 
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The selection rules corresponding to a pertur- 


bation of the type II (13) are'* Am;=0, +1, 
Am;=0, +1. Since transitions of the type 
Am,;=-1 result in very small energy changes, 
and therefore correspond to very low oscillating 
field frequencies, they are of comparatively 
little interest. We, therefore, confine our atten- 
tion to the resonance spectrum corresponding to 
the transitions Am;= +1. 

Returning, now, to the case J=1, for which 
we have already calculated the energy values, 
we write for the energy differences corresponding 
to the selection rules Am; = +1, Am;=0. 


W(1, 0) = wogrH + $e*gQ(3z*—1) 
(327+ 1)(1—2?) 
wogrH(1—e) 
W(0, —1) =ogrH — $e?qQ(32?—1) 
(32+ 1)(1—z*) 
wogrH(1—«) 


In Fig. 3 is plotted W/yogrH as a function of z 
for different values of the (small) variable 
x= $e?qQ/uogrH and for e=0. 

Calculation of the shape function for the 
spectrum requires the transition probabilities, 


II (20) 


+3(Se*qQ)* 


Fic. 4. Shape of the resonance line for J = 1, corresponding 
to the transitions Am;= +1, Am,;=0, in strong magnetic 
fields, as a function of the transition energy. The ordinate 
unit is arbitrary, as is the height at which the line is cut off. 
The calculation of shape includes only first order-pertur- 
bation terms. The figure corresponds to a molecular beam 
experiment in which the field is held constant while the 
the resonance. 


which we obtain directly from the matrix ele- 
ments II (18), and the condition H’ = H,’i. 


P(1, 0) =(1, my\|35’|0, my)? 

P(0, —1)=(0, —1, mz)? 
= 

The fact that P is independent of z is caused, 
of course, only by the fact that we are considering 
transitions for which Am;=0. Thus, the shape 
function for this transition depends on z only 
through the factor D(z). 

Considering the value of W(z) up to first-order 
perturbation terms, we have f(z)=(3z2*—1). 
Hence, D(z)=k/6z. In Fig. 4 we have plotted 
the strength of the line (in arbitrary units) as a 
function of the transition energy in units of 
W/ogrH. To do this, we note that z is propor- 
tional (up to terms of first order in x) to the 
square root of the energy, measured from the 


value for z=0. This enables us to write for the 
shape function of m;=1—+m,;=0. 


II (21) 


k 1 
D=- II (22) 
7 
+1] 
$e*gQ 


in which the value of W varies from H — 
to pogril+2X %e'qQ. The value of D for the 
transition m; =0—m, = — 1 differs from the above 
only by the sign of the term 1 in the denominator. 
The final line shape is determined by super- 
posing the shape functions of both transitions. 
The idealized line shape shown in Fig. 4 has 
certain striking features. For Values of W corre- 
sponding to z=0, the denominator of the shape 
function vanishes and the line strength becomes 
infinite. It is also to be noted that the outsides 
of the lines are very sharp, while, on the inside, 
the lines merge into one another. Observed 
molecular beam spectra may be expected to 
differ from this idealized shape mainly because, 
as discussed in footnote 20, a transition has 
natural breadth, the magnitude of which in- 
creases with increased oscillating field. The 
effect of the natural breadth will be to broaden 
the lines and to “‘wash out”’ the sharpness of the 
line shape. The height at which we have cut off 
the line strength is completely arbitrary; in an 
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actual case, this height will depend on the 
strength of the perturbation, on the value of J, 
as well as on other experimental factors. In 
addition, the effect of second-order perturbation 
terms—which were disregarded in our discussion 
of line shapes—will be to distort the shape of 
the line and to shift the positions of the minima— 
without, however, changing the distance be- 
tween them—towards higher energies by an 
amount equal to $x”. 

The case J=3 differs from J=1 only in that 
there is a third transition of energy difference 
intermediate between the other two, and whose 
energy is—to first-order perturbation theory— 
independent of z. The energy differences, in- 
cluding terms of second-order perturbation, are: 

Wi, 1) 


e*gQ\? 
6 ’ 
( 4 ) pogrH (1 —e) 


— 3) 


(1 —27)(92?—1) 

II (23 
i( 4 ) (1 —e) 


Fic. 5. Shape of the resonance line for J = 3, corresponding 
to the transition Am;=+1, Am;=0, in stron netic 
fields, as a function of the transition energy. The ordinate 
unit is arbitrary, as is the height at which the wings are cut 
off. The central minimum rises to a height approximately 
four times that of the wings. The shape function includes 


only first-order rbation terms. The figure corresponds 


to a molecular beam e iment in which the field is held 
constant while the oscillating frequency is varied through 
the resonance. 


e? 


) 


Corresponding to these differences, the transi- 
tion probabilities are : 


4) = 
P(3, —3) II (24) 
P(—3, 


The shape functions for the transitions (3, 4) 
and (—3, —#) have the same form as those 
calculated for the J=1 transitions—provided 
only the terms of first-order perturbation are 
taken into account. In addition, in this approxi- 
mation, the central minimum is independent of z 
and hence the breadth of the central minimum 
will be determined only by the natural width 
and can be made quite small compared to the 
breadth of the wings. The total probability of 
the central line (and, hence, the area under- 
neath) is greater by a factor of 4/3 than for 
each of the other two lines. 

In Fig. 5 the shape functions of the various 
transitions are superposed to show the shape of 
the spectrum for J=%. The width of the central 


- minimum and the height at which the others cut 


off are arbitrary. As was the case for J=1, 
finite transition, widths will cause a general 
broadening of the lines as well as a loss of sharp- 
ness of the contours. 

In this case, second-order perturbation, be- 
sides broadening the lines, will cause the central 
minimum to shift towards higher energies by 
2x7; the wings, on the other hand, suffer only a 
distortion of shape unaccompanied by a shifting 
of the position of the minimum. 

The idealized spectral line shapes for the cases 
I=2 through 9/2 are shown in Fig. 6. As J in- 
creases, the structure of the wings becomes 
increasingly complicated. However, in an experi- 
ment of finite resolving power, it is not to be 
expected that this structure will be resolved. 
Instead, most structure will probably be run 
together, so that the patterns will be similar to 
those shown for J=1 or I=}. 
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Fic. 6. Idealized shape of the resonance lines for J=2 through 9/2, corresponding to the transitions Am;= +1, Am; =0, 
in strong magnetic fields, as a function of the transition energy. The ordinate unit is completely arbitrary, as is the height 
at which the central lines and wings are cut off. The shape function includes only first-order 
figures correspond to molecular beam experiments in which the field is held constant while 
frequency is varied through the resonance. 
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Ill. WEAK MAGNETIC FIELDS; ZEEMAN EFFECT 
A. Zero Field Energy Levels 
(1) Semiclassical Derivation 
In zero field, the coupling between nuclear spin 
and molecular rotation is strong; we may adopt 
the Russell-Saunders approximation. In this 


scheme, we define a total angular momentum 
F=I+J, where F may assume the values 


Pu S43; --+, | 
F=J+i. | III (1) 


For large J, since J and F so nearly coincide, we 
may consider I to be precessing about J, with 
angle between I and J (Fig. 2) given by 


cos Oy=i/|I|. (2) 


Adopting the same method, as in Section II (1), 
for the calculation of the quadrupole term in the 
energy expansion, one arrives again at Eq. II (4) 


E(Q) = 4e*gQN(§ cos’ III (3) 


As a result of the quantum mechanical nature of 
the vector I, it is again necessary to introduce 
the normalization factor 
Hence 


(2) Quantum Mechanical Solution 


The quantum mechanical expression for the 
quadrupole energy is given by Casimir" 


(5) 


E(Q)= 


where 
K=F(F+1)—I(I[+1) —J(J+1). III (6) 
Using the relationships F= J+7; J>zi, I, we have 
e*gQ [31?—I(I+1) ] 


a result identical with that deduced by semi- 
classical methods. 
Since E(Q) only involves 7 as 7, there are in 


zero field, for given J, 
I+1 levels—for I an even multiple of 3, 
I+4 levels—for J an odd multiple of 3, 


independent of the value of J (provided J is 
large). For transitions with the selection rule 
Ai= +1,” the energy differences are 


W(t, III (8 
Thus, for J=1 
3e*gQ 
W(+1, 0)= III (9) 


and for 
e*gQ 
W(, 3)= W(-3, III (10) 
etc. 


From the above considerations it may be seen 
that a measurement of the zero-field molecular 
beam spectrum of a diatomic molecule in which 
one nucleus has electric quadrupole moment, can 
yield considerable information. The value of the 
quadrupole moment—or, at least, of g?@—can be 
determined directly, by a measurement of the 
frequency of the lowest (or any other) transition. 
In addition, the nuclear spin can be obtained by 
counting the number of transition lines appearing 
in the spectrum. A mere counting is not, however, 
unambiguous, since even nuclear spin and odd 
nuclear spin } integer greater yield the same 
number of transitions. If in addition to counting 
transitions, one measures the ratios of the energy 
differences [see III (8) ] involved, one then has 
an unambiguous determination of nuclear spin. 

Molecular beam resonance transitions, how- 
ever, cannot be observed in zero magnetic field. 
The necessity for a field arises from the fact, as 
explained in Section I, that in order to detect 
changes in magnetic moment, the molecules must 
pass through inhomogeneous fields before and 
after their passage through the homogeneous 
field. If the strength of the central, homogene- 
ous field is very much different from that of the 
neighboring inhomogeneous fields, there will be 

21—n the section to follow, it will be shown that the 
selection rule for Russell-Saunders coupling is AF=0, +1. 
For given J, in the zero-field approximation, this is equiva- 
lent to Ai=0, +1. However, At=0 does not lead to ap- 
preciable energy differences and does not interest us. Thus, 


we — the selection rules for transitions in zero field 
Ai=+1,. 
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induced “‘non-adiabatic transitions’ by the sud- 
den field changes. To avoid this, the magnitudes 
of all fields must be made approximately equal— ; 
thus, the homogeneous field cannot easily be 
reduced to zero. 

The introduction of a small magnetic field into 
the problem results in a perturbation and splitting 
of the energy levels. This problem (the Zeeman 
effect) is treated in the following sections. 


(F, m|I|F+1, m+1) 
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B. Weak Field Energy Levels and Transitions 


In weak fields, the Hamiltonian for our system 
contains, in addition to the field-independent 
electric quadrupole term already discussed, a per- 
turbation of the form, pogr/7I.+uogsHJ,.. In the 
F, mr™ representation, the quadrupole term is 
diagonal with elements given by Eq. III (4). The 
matrix elements for I and J in this representation™ 
have the following form: 


2 


| 
4(F+1)*(2F+1)(2F+3) 


(F, m|1| F+1, m)=| 


(F, m|1I| F+1, m—1) 


4(F+1)*(2F+1)(2F+3) 


2 
I(I+1) -—J(J+1)+F(F+1) 
2F(F+1) 
2F(F+1) 
I(I+1) —J(J+1)+ F(F+1) 
2F(F+1) 


(F, m|I| F, m—1)= 


(F, m|1|F—1, m+1) 


+ 
4(F+1)*(2F+1)(2F+3) 


(i—7j), 


mk, (11) 


2 


(F, F—1, m)= 


| 
4F*(2F—1)(2F+1) 


4F?(2F—1)(2F+1) 


(F, m|1| m—1) 


2 


The matrix elements for J are obtained as 
follows: for AF=0; interchange J and J in the 
matrix elements for I. For AF= +1; the matrix 
elements for J= —the corresponding matrix ele- 
ments for I. 

The matrix elements, in the F, m represen- 
tation, for the small oscillating perturbation 


| 
4F*(2F—1)(2F+1) 


ia. 


are obtained from the 
above. The selection rules for this perturbation 
are seen to be AF=0, +1, Am= +1. However, we 
shall not consider the transitions AF=0. For 
these only involve changes in m; and since the 


2 Henceforth mp shall always be referred to as m. 
* Referencei3, pp. 63, 64, 67. 
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Fic. 7. Energy of transition for J=1 in weak magnetic 
fields as a function of z=m/J, for different values of the 
perturbation x. The two sets of curves correspond to the 
transitions F= J+1—F=J and F=J—1—F=J. The four 
transitions of Eqs. III (14) reduce to two because of the 
condition e=0. 


main (electric quadrupole) term in the energy is 
independent of m, these transitions will not in- 
volve the quadrupole moment. In addition, since 
AF=0 transitions are concerned with differences 
between small perturbation terms, the energy 
differences will be very small compared to transi- 
tions for which AF= +1. . 

In a fashion completely similar to that de- 
scribed in the preceding sections, we may calcu- 
late energies, energy differences, transition proba- 
bilities, and line shapes for the system in weak 
fields. 


(1) 


(a) Using perturbation theory including terms 
of second order, we find the energy values are: 


e*gQ 
E(J+1, m) 
+ (uogrH — 


H- 2° 
4 (1—s9), 


—e 


E(J, m)= + pogsHm 
(wog rT — wogsH)?* 


(1—z2?), III (12) 


E(J—1, m) 


— (wogrH — wogsH)z 
— wogsH)? 


The continuously varying parameter z is used for 
the quantity m/J. The variable expressing the 
strength of the perturbation we denote by 


x= pogrH = wogrH /Wo. III (13) 


We have seen that!’ We write gy = 
Although we retain terms of order « in all of the 
following expressions, most of the calculations 
shall be made for e=0. 

(b) From the expressions for the energy values 
and the definitions given above, the energy 
differences for the selection rules AF=+1, 


(1—2?). 


919? gy (€*0) 


2g,(€ = Y2) 


WA, 


Fic. 8. Shape of the resonance lines for J=1, corre- 
sponding to the transitions AF=+1, Am=+1, in weak 
magnetic fields, as a function of the transition energy. The 
shape function is calculated to terms including only first- 
order perturbation. 8 (a) is the shape function for e=0, and 
8 (b), for e=4. The figures correspond to molecular beam 
experiments in which the magnetic field is held constant 
while the oscillating field frequency is varied through the 
resonance. 
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Am= +1 follow: e 
W(I+1, m; J, m+1) 
=1+[(1—e)z— 

(12) Wo § 

0 

= 

3(1—z*)(1—¢) 

W(J—1, m; J, m+1) 5 

for W(J—1, m; J, m—1) 

the 

|* 3 
13) In Fig. 7 are plotted the energy differences as a { = 2 8 5 
egr. function of z, for e=0. The energy differences are % 2 ° 
the shown for different values of the perturbation * es 
oa x—which is varied from 0.01 to 0.25. * sé 
(c) Corresponding to the energy differences, 
lues are the transition probabilities : 
- 
P(I+1, m; J, m+1)=(J+1, m|H|J,m+1)? 
P(J+1, m; J, m—1) 2} 

= 

P(J—1, m; J,m-—1) 

The shape of the spectrum is obtained by tak- © : 2 

ing the products of P by D=1/(dW/dz). If we 7, 

consider first-order perturbation termsonly,wesee 

that W depends linearly upon z; the factor Dis, S&X a 

then, independent of z. Hence, the shape function © 3 a § 

is simply proportional to P. In Figs. 8(a) and 8(b)  — " 

are plotted the shape of the spectral lines for the a £ 

transition AF= +1, Am= +1. Figure 8(a) repre- 2 

sents the case 8(b) hasgr=2gy.Theshape 

The of the spectrum for a value of ¢ between 0 and} ™“ jw § 

a will be intermediate between these extremes. ° o 

In an experimental measurement of the spec- 

trum, the observed shapes will differ considerably ° 


from the idealizations shown. Deviations will 


s 
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result from failure to satisfy the assumptions 
upon which the calculations are based; the per- 
turbation caused by the homogeneous magnetic 
field may, for instance, be so large that it is no 
longer sufficient to consider only the first-order 
term in determining line shapes. In this case, 
introduction of second and higher order pertur- 
bation terms will cause the line shape to be a 
more complicated function of z, and of energy; in 
general, lines will be somewhat broadened and 
distorted by this effect. Increase in the magnitude 
of the oscillating field will also introduce compli- 
cations in the values for the transition proba- 
bilities P,?° and similarly, result in a distortion 
and broadening of the spectral lines. In addition, 
the velocity distribution in the molecular beam 
will add to the broadening and diffusion of the 
shape of the line, since a spread in the times spent 
in the oscillating magnetic field introduces a 
corresponding spread in the transition proba- 
bilities for molecules in the beam. As a result of 
these factors, the observed spectral line will 


—2(Ws44|5C| 


probably be rather flat and broad, with width 
about 2yogrH. 


(2) I=} 


(a) Energy values. Calculation of the energy 
values for J = % by perturbation theory is compli- 
cated by the fact that the unperturbed (zero field) 
energy levels are degenerate [see Eq. III (4) ]. 
More important, there exist non-vanishing ma- 
trix elements between the degenerate states 
F=J+3 and F=J—}. Hence, it is necessary to 
apply degenerate perturbation theory to the 
problem.* To do this, we define two new states 
characterized by the wave functions VY, and 
W_,; these are linear combinations of the original 
4, and WV,_,, chosen to satisfy the requirement 
that the matrix element (¥,/3C|_)?* vanish. 
The combinations 


V,=sin 


WV_= —cos OV 


III (16) 


where 


III (17) 


tan 20= 


satisfy this requirement, as well as the normalization and orthoganality conditions. We may apply 
second-order perturbation theory to these new states, together with the original J+ 3 and J—} 


states, and obtain the following energies: 


2 1—¢)*(1—z? 
4 
2 272 1-— 2 2 
— 6) *(1 III (18) 
— €)*(1—2?) 


e*gQ 


(b) Energy differences: Adopting the selection rule AF=+1, Am= +1, there are eight possible 
transitions (involving the quadrupole term), with energy differences : 


m; +,m+1) 
Wo 

W(J+4%, m; +, m—1) 
Wo 


=1+[(1—) — ee +3 (1 — 6) °(1 


*4 Reference 13, pp. 30-24. 
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Ve 1= 3% 
8 


Fic. 10. Shape of the resonance lines for J] =} corresponding to transitions AF= +1, Am= +1, in weak magnetic 
fields, as a function of the transition energy. The shape functions include only first-order perturbation terms. 10 (a) is for 
e=0, while 10 (b) is for «=}. 


W(J+4, m; —,m+1) 


Wo 
W(J+%, m; —,m—1) 
W —2*)x’, 
WI-3 m +1) 
—3,m; 
W(J- ’ ;+, —1 
0 
W(J—3, m; —,m+1) 
=1—[(1—¢) {§2—(1 — $2") tee +9 (1 —6)*(1 x", 
0 
W(J—i, m—1) 
where 
— » €=g7/g1, s=m/J. TIT (20) 


These energy differences are plotted in Fig. 9 as a function of z, for e=0, and for different values of 


the perturbation parameter x. 
(c) The transition probabilities corresponding to the above energy differences are: 


m; +, m1) = 1+ 


2(1 


PU+4, m; —, m1) = 


— 
III (21) 


mi +, mse) = 


P(J—4, m; —, m+1) = 1+ 


2(1—32%)4] 
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Jal) WJel,J) 


J,J+1) B(J) 


0 


© WI-1,J) B(J-2) 


The product DXP (considering only first- 
order perturbation terms) turns out to be par- 
ticularly simple in this case. In fact, the shape 
function may be obtained from the transition 
probability by dropping the last factor. Thus, 
DXP(J+3, m; +, m+1) =k(1—2)?, etc. 

In Figs. 10(a) and 10(b) are plotted the line 
shapes for e=0 and «= }, respectively. Again, we 
have assumed monoenergetic molecules and small 
perturbation. The same factors as have been 
discussed for the case J=1 will again come into 
play in this case to distort the line shape. 

The shapes of the spectra, for the cases J=2 
through J=9/2, have been investigated. Corre- 
sponding to each line in the zero-field spectrum, 
there is a group of transitions. For integral values 
of J, the shapes of the spectral lines are identical 
(except for second-order effects) with those shown 
in Fig. 8. Half-integral values of I involve two 
types of spectral lines. For transitions not in- 


Fic. 11. Complete energy matrix for J=1 in the F, m representation. 


B( Jel) W(J41,J) 
B(J) W{J,J-1) 


© MJ-1,J) B(J-1) 


volving the degenerate states (F~J+}4), the 
shape of the spectrum will again be as shown in 
Fig. 8. For transitions involving the degenerate 
states, the shape will—in general—be consider- 
ably more complicated (even than Fig. 10), 
although the general nature of the spectral group 
is like that shown in Fig. 10. However, in an ex- 


‘ perimental observation of the spectrum the more 


complicated structure will probably remain 
unresolved. Hence, we shall omit a more detailed 
discussion of these transitions. 


IV. INTERMEDIATE FIELDS; /=1, 
EXACT SOLUTION 


A. Energy Values 


In the representation of our system for which 
F and m are the appropriate quantum numbers, 
the quadrupole term in the Hamiltonian is 
diagonal. The only off-diagonal terms in the 


E 
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energy matrix arise, for constant m, from the 
matrix elements of J, and J, in this representa- 
tion. The complete energy matrix takes the form 
shown in Fig. 11 where the diagonal terms are 
represented by £ and the off-diagonal terms 
by 

Thus, the matrix divides up into a series of 
three-row, three-column submatrices (except the 
ends). Since there are no connecting elements 
between different submatrices, each one of these 
may be diagonalized separately. Hence, the 
secular equation for our problem takes the 
general form: 


B+b-—E 0 
c 
c B-—b-E 


=0, IV (1) 


p= + 


—e 
a= 


+uogsHm, 


b=ogrH(1—e)z, 
c= wogrH (1 —€) {}(1—2?)} 4. 


This determinant reduces to a cubic equation 
having, in general, three different, real roots. 


e 


qQ 


—e? 


These may be found by a straightforward alge- 
braic process. The general results are rather 
complicated, and will be omitted here. 

In special cases, the solution takes on a rather 
simplified form. ; 


(1) m=J 


For this condition, c=0; upon expansion of the 
general expressions for the energies, we obtain for 
the roots: 


Er;=a, 
En 
=8+0. 


Since, for c=0, the energy matrix is diagonal for 
all values of the field [Eq. JV (1)], we are 
enabled to identify the roots with the F,m 
quantum numbers. This simple comparison of 
roots gives 


IV (3) 


E,;=E(J,m), 
En=E(J— m), 
Emr = E(J+1, m). 


IV (4) 


(2) Weak Fields 


For small field perturbation, [b, c<(8—a) ], 
the solutions reduce to: 


+4(1/2—62+5z') 


(1—z*) 


—e)* 
(3e*gQ)? 


E(J, m) = E;= + — 


e*qQ 


(1 

H*(1 — 
(3e%gQ)? 


the third-order term in the weak field expansion 
diverges. Furthermore, the third-order term in 
the perturbation theory expression for the energy 
is also found to exibit the same property. The 


These results are seen to agree with the energies 
obtained from second-order perturbation theory, 
in Section III. 

It is interesting to note that for m=0 (z=0), 


| 
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reason for this behavior becomes clear when we 
consider that we are dealing with a degenerate 
system. Even though there are no matrix ele- 
ments connecting states of equal unperturbed 
energy, it is not permissible to apply perturbation 
theory, unless the degeneracy is removed by the first 
or second-order perturbation.*® That this condition 


—e? 


E(J+1, 0)= 


E(J—1, 0)= 
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(4) Strong Fields 


does not hold for m=0 may be observed, if we 
consider the energy matrix for b=0 [IV (1)]. 
However, for the case )=0, it is easy to solve the 
secular equation exactly. 

(3) m=0 


Solution of the equation gives 


IV (6) 


3 { — )?(1—2*)} 4. 


For this approximation, we may consider the quadrupole interaction small compared to the 
magnetic interaction. Thus, }, c>(@—a). Expanding the exact solution, we arrive at the following 


energy values: 


E(J+1, m) = —wogrH + wogsH 1) — — 32?) — 


($e*gQ)? 
(1 — e) 
—4z4(1 —2z?) 


-(1—27)(1+ 32?) 
(2e*gQ)? 


IV (7) 


E(J, m) = wogsHm-+ — 32”) + 824(1 —2?) 


E(J—1, m) = wogrE + woes H(m —1) —}e%qQ(1 — 32%) +4 


The representation in which it is convenient 
to discuss the Paschen-Back effect has as 
quantum numbers m; and my. To achieve the 
transition between representations, we recall 


that our matrix elements have been calculated . 


for constant m. Thus 
m=m,-+-m,=constant. IV (8) 


Hence, if mr=1, m;=m—1; likewise for m;=0, 
myz=m; and for —1, my=m+1. These re- 
lationships provide the connection between repre- 
sentations, and we get 


E(J+1, m)=E(—1,m,), 
E(J,m)=E(0, ms), 
E(J—1, m)=E(1, my). 
% Reference 13, first paragraph, p. 34. 


IV (9) 


($e*qQ)? 


(1—e2 2 
2”)(1+32?) 


—€)? 


—424(1—2?) 


With these identifications, we see that the energy 
values calculated above agree exactly with those 
obtained in Section II, by perturbation theory. 


B. Eigenfunctions 


The wave functions V, of our system in uni- 
form magnetic field may be expressed in terms 
of the set of wave functions ®r,, for which the 
quadrupole term is diagonal. 


Vi Ps, TV (10) 
These obey the wave equation 
i=1,2,3 IV (11) 


and also form a normal, orthogonal set. 
Using the values of the energy Ej, as calcu- 
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lated in the preceding parts of Section IV, 
together with the normalization condition, we 
may use the wave equation to derive the values 
of the coefficients. 


C. Transitions 


From the calculated wave functions and 
matrix elements in the F, m representation, we 
may compute the matrix elements corresponding 
to transitions in our system. To do this, we recall 
that the matrix elements of the oscillating 
perturbation (in the F, m representation) allow 
transitions with the selection rule AF=+1, 
Am=-+1. If we assume small perturbation, the 
transition probability between states will be 
proportional to the square of the matrix element. 
These are found to be 


P(i, m; j, M41) = (Wi, m| | Vj, m+1)? 
= — €)?{ 
+ 
4, j=1,2, 3. 


IV (12) 
where 


From these and from the expressions for W, the 
shape function DXP can be calculated in any 
specific case. 


V. Discussion 


As a result of the calculations contained in 
this paper, it appears possible that a molecular 
beam experiment can yield valuable information 
concerning nuclear electric quadrupole moments. 
We have shown in Section II that a determina- 
tion of the (constant) energy difference between 
subsidiary minima, in a strong magnetic field, 
yields the quantity e*gQ. The fact that such 
minima have been observed for the resonance of 
Li? in LiF (Fig. 1) leads to the expectation that 
the experimental conditions for their observation 
may be favorable. 

In Fig. 1, the position of the subsidiary minima 
is difficult to determine. However, for a rough 
estimate, we shall consider these to be at 18.70 
and 19.05, on the scale of the abscissa. This 
separation corresponds to‘a field difference®’ of 

~The ex iment was performed by keeping the fre- 


uency of t e perturbing oscillator constant and varying 
the iemageneess magnetic field through the resonance. 
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47 gauss—about 10 times the natural width of 
the line. Hence, with improved resolving power, 
it may be expected that the subsidiary minima 
will be well resolved. 

For the temperature at which the LiF béam 
was formed (1100°K), the most probable value 
of J is about 16. This is sufficiently greater than 
3, the spin of Li’, so that we may consider the 
approximation developed in this paper to apply. 
According to this approximation (see Fig. 5) 
the ratio of the field difference between minima 
to the field value at the central resonance is 


From this, the value of e*gQ for Li’ is calculated 
to be equal to 1X10-* erg. Although the value 
of the quadrupole moment Q cannot be deter- 
mined until g is calculated, the above figure may 
be compared with that found for deuterium in 
HD and D:.® This was about 0.3 10-* erg and 
corresponded, in the case of D, to a quadrupole 
moment of 2.7 cm?. 

The zero-field spectrum, if it could be observed, 
would give a direct and simple measurement 
of gQ. However, the experimentally necessary 
introduction of a small magnetic field leads to 
complications. For magnitudes of the field much 
less than the quadrupole interaction (Section 
III), the lines observed would not be very broad ; 
a determination of their centers of gravity would 
be sufficient to yield the information we require. 
However, when the field is of the sare order as 
the quadrupole interaction (about 50 gauss in 
LiF and Nag, 100 gauss in K,)*® the spectrum 
becomes complicated with flat, broad, uninter- 
pretable minima. It is hoped that there may soon 
be presented an opportunity to apply these 
calculations to the experimental determination 
of nuclear electric quadrupole moments. 

The authors wish to acknowledge their in- 
debtedness to Professor I. I. Rabi, who first 
suggested this problem, to Professor Edward 
Teller, for his continued and inspiring interest 
and discussions, and to the members of the 
Molecular Beam Laboratory for their help in 
the preparation and criticism of this manuscript. 
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In bombarding thin beryllium targets with protons of energy 0.30 Mev to 3.0 Mev, gamma- 


ray resonances were observed at proton energies 0.975, 1.06, 1.13, 1.36, and 2.52 Mev and a 
neutron resonance at 2.52 Mev. For 1.04-Mev protons the value obtained for the gamma-ray 
energy was 7.5 Mev which indicates a radiative capture process. At the 2.52-Mev resonance the 
measured value of the energy of the gamma-rays was about 3.0 Mev. The most probable expla- 
nation of the observed resonance at 2.52 Mev appears to be that it is a resonance for the capture 
of a proton to form B'* which later disintegrates into (B®+-n') and also into (Li**+He*) where 


Li** emits 3-Mev gamma-radiation in returning to the ground state. 


INTRODUCTION 


AMMA-RAYS from beryllium caused by 
proton bombardment have been observed 
by several investigators. Crane, Delsasso, Fowler, 
and Lauritsen! observed gamma-rays from a thick 
beryllium target for proton energies 0.45 Mev to 
0.80 Mev. They detected no resonances but found 
the energies of four gamma-ray lines to be 2.2, 
3.7, 4.8, and 6.0 Mev. Hafstad and Tuve? found 
no resonances for proton energies 0.40 to 0.90 
Mev. Herb, Kerst, and McKibben?’ extended the 
work to 1.6 Mev and obtained a broad resonance 
at approximately 1.0 Mev. Curran, Dee, and 
Petrzilka* examined the excitation curve with 
proton energies of 0.20 Mev to 1.0 Mev for thin 
beryllium targets of various thicknesses. They 
reported resonance peaks at 0.35 Mev and at 0.67 
Mev. Measurements of absorption of the gamma- 
radiation by the same investigators showed that 
for proton energies of 0.40 Mev to 0.80 Mev the 
gamma-radiation was caused chiefly by the 
radiative capture process. 

In the present investigation. the excitation 
curve was studied for a thick beryllium target for 
proton energies 0.90 Mev to 1.5 Mev and for thin 
beryllium targets in the range 0.30 Mev to 3.0 


* Now at the National Research Council of Canada, 
Ottawa, Canada. . 

1H. R. Crane, L. A, Delsasso, W. A. Fowler, and C. C. 
Lauritsen, Phys. Rev. 47, 782 (1935). 

*L. R. Hafstad and M. A. Tuve, Phys. Rev. 48, 306 
(1935). 

*R. G. Herb, D. W. Kerst, and J. L. McKibben, Phys. 
Rev. 51, 691 (1937). 

4S. C. Curran, P. I. Dee, and V. Petrzilka, Proc. Roy. 
Soc. A169, 269 (1938). 
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Mev. Measurements of gamma-ray energies were 
made at each of two resonances. 


EXPERIMENTAL 


Thin beryllium targets were made by evapo- 
rating beryllium metal within an evacuated brass 
chamber. The beryllium was heated by placing it 
in a molybdenum or tantalum boat which was 
brought to white heat for a fraction of a minute 
by passing a current of forty to fifty amperes 
through the boat. Some difficulty in evaporating 
beryllium is encountered in this method if a 
V-shaped ribbon is used for a boat since the 
beryllium forms a low melting point alloy with 
molybdenum and tantalum so that the ribbon 
melts before the beryllium evaporates. The boats 
used were made by taking a ribbon of molybde- 
num or tantalum metal about 8 mm wide and 
5 cm long and rolling it into a cylinder about 3 
mm in diameter and 5 cm long. The ends were 
then pinched so that the cylinder forms a boat 
which encloses the beryllium completely except 
for a slot about a millimeter wide along the top. 


- Although the beryllium alloyed with the tantalum 


or molybdenum, two or three evaporations could 
be carried out with the same boat before the 
reaction penetrated the boat so as to make holes 
in it. The beryllium was evaporated onto a sheet 
of molybdenum or tantalum target-backing from 
a boat of the same material. 

The protons were accelerated by the Wisconsin 
concentric-electrode electrostatic generator’ and 
the intensity of the gamma-rays was measured by 


®R. G. Herb, C. M. Turner, C. M. Hudson, and R. E. 
Warren, Phys. Rev. 58, 579 (1940). 
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Fic. 1. Gamma-ray and neutron intensities from beryllium as a function of the energy of the bombarding protons. 


Curve E represents the neutron yield; the remaining curves represent gamma-ray yields. Above the proton ene 
the gamma-ray intensity scale is reduced by a factor of fifty. The bac 


1.9 Mev 
ound, taken on tantalum and molybdenum 


targets, is shown at the bottom of the figure by points corresponding to the upper curves. 


G-M counters with a target chamber and counter 
arrangement similar to that used by Plain, Herb, 
Hudson, and Warren.® 


PROTON ENERGY SCALE 


The proton energy scale for the thin target 
curves was calibrated against the 2.03-Mev 
Be(p, ) threshold by running a neutron yield 
curve simultaneously with the gamma-ray yield. 
This value of the threshold had been reported by 


Haxby, Shoupp, Stephens, and Wells’ and had 


*G. P. Plain, R. G. Herb, C. M. Hudson, and R. E. 
Warren, Phys. Rev. 57, 187 (1940). 

™R. O. Haxby, W. E. 7 , W. E. Stephens, and W. 
H. Wells, Phys. Rev. 58, 103 (1940). 


also been measured in this laboratory* by com- 
paring it to the 0.440-Mev Li(p, y) resonance. 
For the thick target curves the generator voltage 
was calibrated against the 0.862-Mev 
resonance which had previously been compared 
with the 0.440-Mev Li(p, y) resonance. 

An absolute measurement of the Li(, ) 
threshold has been made recently by Hanson and 
Benedict® which, if correct, would revise the 
proton energy scale. Based on this new determi- 
nation the Be(p, m) threshold would be at 2.06 
Mev instead of the accepted 2.03 Mev. In this 


* A. O. Hanson and D. L. Benedict, Phys. Rev. 65, 33 
(1944). 
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VIELO 


104-MEV PROTONS 


2.52-MEV PROTONS 


COINCIDENCE 


OL. i i i 1 i 1 L 


1 2 3 4 
ABSORBER THICKNESS (ma ot Al) 


Fic. 2. Absorption curves for secondary electrons caused 
by gamma-rays from llium for 1.04-Mev and 2.52-Mev 
bombarding protons. The curves are extrapolated back, 
0.30 mm, to correct for the absorption by the glass walls of 
the counters. The short horizontal lines indicate the thick- 
ness of aluminum which reduced the yield to one-half. 


revised scale all values of proton energies would 
be raised by 1.5 percent above the previously 
accepted values which have been used in this 


report. 
EXCITATION CURVES 


The gamma-ray yield as a function of the 
energy of the bombarding protons is shown in 
Fig. 1 with the yield in arbitrary units. The thick 
target curves A and B show resonances at ap- 
proximately 0.975, 1.06, and 1.16 Mev. Curves C 
and D show the gamma-ray yield from a thin 
beryllium target on molybdenum backing, the 
molecular beam being used for C and the atomic 
beam for D. Curve E shows the neutron yield 
obtained simultaneously with the gamma-ray 
yield D. Curves F and G show the gamma-ray 
yield from a thin beryllium target on tantalum 
backing, the molecular beam being used for F, 
_the atomic for G. The background yield was 
taken on a molybdenum target for curves C, D, 
and E and on a tantalum target for F and G, and 
is shown by the corresponding points at the 
bottom of Fig. 1. Above the proton energy 1.9 
Mev, the gamma-ray intensity scale is reduced 
by a factor of fifty. The thin target yields show 
resonance peaks at proton energies 0.975, 1.06, 
1.36, 2.52, and possibly at 0.86 and 1.13 Mev. 


The possibility that the hydrogen-filled G-M 
counters may have been counting neutrons was 
tested by placing four cm of paraffin between the 
target and the counter. The count was reduced 
by less than one percent showing that neutrons 
do not make any substantial contribution to the 
measured gamma-ray yield. 


ENERGY OF GAMMA-RAYS 


The energy of the gamma-rays was measured 
by obtaining the coincidence yield of two G-M 


- counters as a function of the thickness of alumi- 


num absorbers placed between the counters. 
Figure 2 shows the results of the absorption 
measurements of the secondary electrons caused 
by the gamma-rays. The abscissa scale gives the 
thickness of aluminum absorber. To correct for 


the absorption by the glass walls of the counters - 


the curves are extrapolated back, 0.30 mm, to the 
short vertical line. This correction figure had 
been determined by Plain, Herb, Hudson, and 
Warren.* The thickness of aluminum absorber 
which reduces the coincidence yield to half its 
viaue for no absorber is used as a measure of the 
gamma-ray energy and is shown by short hori- 
zontal lines in Fig. 2. The values are 3.3 mm Al 
for 1.04-Mev protons and 1.5 mm Al for 2.52- 
Mev protons. 

The counters had been calibrated previously® 
by using the 6.2-Mev gamma-rays from fluorine 
and the 17.5-Mev gamma-rays from lithium 
giving half-value thicknesses of 2.8 and 7.3 mm 
Al, respectively. Figure 3 shows the resulting 
calibration curve. If we assume that the radiation 
is monochromatic, the energy of the gamma-rays 


HALF-VALUE ABSORPTION THICKNESS 


GAMMA-RAY ENERGY - MEV 


Fic. 3. Calibration curve showing half-value absorption 
thickness as a function of gamma-ray energy. 
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GAMMA-RAYS FROM Be 


from beryllium is found to be 7.5 Mev for 1.04- 
Mev protons and about 3.0 Mev for 2.52-Mev 
protons. 

Placing 3-mm lead between the target and the 
counter did not reduce the intensity of the 7.5- 
Mev radiation appreciably, showing it has no 
large soft component. 


DISCUSSION 


Several resonance maxima appear in the excita- 
tion curve of Fig. 1 although the group at 1 Mev 
may not be completely resolved. The maxima at 
0.35 Mev and 0.67 Mev reported by Curran, Dee, 
and Petrzilka* were not observed. 

A summary of the results of the absorption 
measurements is given in Table I where the 
values in brackets are obtained if the new proton 
energy scale is used. The reaction energy given 
in the table assumes a radiative capture process, 


i.e., 
Be®+H'!—B!°+ Ap. (1) 


This reaction has 6.56 Mev of available energy 
plus the resonance energy of the protons. For 
1.04-Mev protons the total energy available is 
6.56 plus 9/10 of 1.04 or 7.50 Mev. This agree- 
ment between the gamma-ray energy and the 
total energy available shows that the reaction is 
resonance capture with a single transition to the 
ground state. 

The reaction at the 2.52-Mev resonance is more 
difficult to interpret. In a radiative capture reac- 
tion the available energy would be 6.56 plus 9/10 
of 2.52 or 8.83 Mev, whereas the energy of the 
radiation is about 3.0 Mev. The intensity of the 
radiation was roughly fifty times greater than the 
intensity at the 0.975-Mev resonance for the thin 
target G and about one hundred times greater for 
the thin target D. For the latter target the neutron 
yield, which also shows a resonance at 2.52 Mev, 
was estimated roughly to be of the same magni- 
tude as the gamma-ray yield. This suggests that 
the neutrons and the gamma-rays are produced 
in the same reaction. 

Some other possible reactions are given below. 
In calculating the reaction energies the following 
atomic mass values given by Haxby, Shoupp, 
Stephens, and Wells’ are used ; Be®= 9.01484 m.u., 
B*=9.01600 m.u., (m'—H')=0.000806 m.u. 
Other mass values are taken from Barkas.® 


*W. H. Barkas, Phys. Rev. 55, 691 (1939). 


A neutron producing reaction is 
(2) 


where the reaction energy Q2 is —1.83 Mev. 
However, B® is unstable and may return to Be® 
by emitting a positron or capturing a K electron 
or it may disintegrate into Be® plus H'. In the 
first case 


where the maximum positron energy E, is 0.06 
Mev. For K-electron capture, reaction (2) above 
is followed by 


with an energy evolution Q,’ of +1.08 Mev. In 
the third case 


B*—Be®+H'+Q2”, 
where Q,”’ is +0.10 Mev. If in place of the latter 


TABLE I. Gamma-ray energies. 


Reaction en- 


Proton Half-value Gamma 
energy (Mev) (mm of Al) energy (Mev) re Mev) 


3.3 7.5 


7.50 
(7.51) 
8.83 
(8.86) 


1.04 
(1.055) 

2.52 1.7 3.0 
(2.56) 


the non-capture disintegration 


takes place, then Q,”’ is —1.73 Mev. 
Theoretical considerations’® show that a par- 
ticle producing process is much more probable 
than a radiative capture process and since this 
resonance is above the Be(p,m) threshold it 
would be expected that the (p, ) reaction would 
be the most probable. The observed gamma- 
radiation would have to come from an excited 
residual nucleus in the Be(p,m) reaction but 
consideration of the energy evolution in the 
above three cases of reaction (2) shows that 
3.0-Mev gamma-rays are energetically impossible 
for 2.52-Mev incident protons. It seems certain, 
then, that the 3.0-Mev gamma-rays and the 
neutrons are not produced in the same reaction. 


10H. A. Bethe, Rev. Mod. Phys. 9, 69 (1937). 
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Two other processes have been observed by 
numerous investigators for low energy protons. 
They are 


(3) 


and 
Be®+H!—Li'+ Qu, (4) 


where Q3;= +0.44 Mev and Q4=+2.03 Mev if 
Be’ and Li® are left in the ground state. 

Among other factors the probability of a 
reaction depends on the energy evolution and on 
the height of the potential barrier to the particle 
emitted. The probability is greatest for a high 
energy evolution and for a low potential barrier. 
The potential barrier is zero for neutrons and 
higher for alpha-particles than it is for deuterons. 
It would be expected that reactions (2), (3), and 
(4) given above are almost equally probable since 
in each case the energy evolution compensates 
for the height of the potential barrier. Skaggs" 
has observed that for 0.262-Mev protons the 
a-particles and the deuterons are emitted in 
almost equal numbers. It is necessary, then, to 
examine whether 3.0-Mev gamma-radiation is 
energetically possible from an excited state of one 
of the product nuclei of reaction (3) or (4). 

In the reaction 


Li’+H'—Be®+Q, 


Delsasso, Fowler, and Lauritsen observed 17- 
Mev and 14-Mev gamma-radiation but no radia- 
tion between 2 and 10 Mev. This was considered 
experimental evidence that Be® has an excitation 
level at about 3.0 Mev but the transition from 
this state to the ground state by gamma-radia- 
tion is forbidden. Rather, the Be’, in the 3-Mev 
excited state, disintegrates into two alpha-par- 
ticles with the excitation energy going into the 
kinetic energy of the alpha-particles. This inter- 
pretation has been confirmed experimentally and 
theoretically more recently by Wheeler and 
several others who place the excitation level at 
2.8 Mev. Furthermore, if the value of the gamma- 
ray energy measured here is correct the radiation 


uL.S. s, Phys. Rev. 56, 24 (1939). 
#%L. A. Delssaso, W. A. Fowler, and C. C. Lauritsen, 


Phys. Rev. 51, 391 (1937). 
13 J. A. Wheeler, Phys. Rev. 59, 16 (1941) and 59, 27 


(1941). 


HUSHLEY 


cannot occur in reaction (3) because of energy 
considerations. 

Assuming that the a-particle has no stable 
excited state, it appears probable that the ob- 
served 3.0-Mev gamma-radiation is caused by an 
excited state of Li® in the Be(p, a) reaction. This 
reaction has 2.03 Mev of available energy plus 
the energy of the incident protons so that it is 
energetically possible to have the Li® excited up 
to about 4 Mev. This could be confirmed by 
determining the energy of the a-particles emitted 
in the reaction. One would expect a-particles of 
about 6.2-mm range from the reaction in which 
Li®™ is formed as well as 17.3-mm alphas corre- 
sponding to the formation of unexcited Li*®. 
These values are calculated from the energy- 
momentum relations given by Livingston and 
Bethe." 

An alternative explanation of the gamma-rays 
is to assume that radiative capture with cascade 
emission is as probable as a neutron producing 
reaction. That is, B!°, excited to 8.83 Mev, would 
return to the ground state in, say, three steps. 
The known excitation levels of B'® are at 0.55, 
2.15, 3.45, and 7.50 Mev so that the 3.0-Mev 
measured value of the radiation would represent 
the average effect of several different y-ray lines. 

The phenomenon at the 2.52-Mev resonance 
may then be interpreted as a resonance for the 
capture of a proton forming B'** which later 
disintegrates in one of several ways, such as 


hy, hv + hv3 
—Li* + He!—Li'+ He*+ Av 
—B*+n'. 


It is expected that the last two are the most 
probable of the three mentioned. 

It would have been desirable to repeat the 
excitation curve using thinner targets and to 
make more extensive absorption measurements, 
but the work was cut short before this could be 
done. 

The author is greatly indebted to Dr. A. O. 
Hanson for advice and help in the experimental 
work and to the Wisconsin Alumni Research 
Foundation for financial assistance. 


4M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 
9, 245 (1937). 
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TZy Electric Circuit Models of the Schrédinger Equation 
able GABRIEL KRON 
~~ General Electric Company, Schenectady, New York 
ree (May 25, 1944) 
This Equivalent circuits are developed to represent the Schrédinger amplitude equation for one, 
plus two, and three independent variables in orthogonal curvilinear coordinate systems. The 
it is networks allow the assumption of any arbitrary potential energy and may be solved, within any 
| up desired degree of accuracy, either by an a.c. network analyzer, or by numerical and analytical 
b circuit methods. It is shown that by varying the impressed frequency on a network of inductors 
y and capacitors (or by keeping the frequency constant and varying the capacitors), it is possible 
tted to find by measurements the eigenvalues, eigenfunctions, and the statistical mean of various 
s of operators belonging to the system represented. The electrical model may, of course, be replaced 
hich by an analogous mechanical model containing moving masses and springs. At first the network 
rre- for the one-dimensional wave equation for a single particle in Cartesian coordinates is developed 
Lis in detail, then the general case. A companion paper contains results of a study made on an a.c. 
‘ network analyzer of one-dimensional problems: a potential well, a double barrier, the harmonic 
rsy- oscillator, and the rigid rotator. The curves show good agreement, within the accuracy of the 
and instruments, with the known eigenvalues, eigenfunctions, and “tunnel” effects. 
rays 
ade THE ONE-DIMENSIONAL SCHRODINGER and capacitors and in each state the currents and 
an EQUATION voltages are sinusoidal in time. The state is 
. I’ Cartesian coordinates, the equation is changed by varying the frequency of the im- 
"4 pressed voltage. 
I ay The basic network concepts will be introduced 
| ve i ee E)y=0. in detail in connection with the first type of 
ae circuit (with additional comments on the second 
ine Three types of equivalent circuits may be type of interpretation). The third type of model, 
‘tie established. although it offers attractive analogies, will be 
noe 1. The circuit contains positive and negative only cursorily treated. 
resistors and in each state the currents and REPRESENTATION OF ENERGY OPERATORS 
voltages are constant in time. The state is , nee 
changed by varying the resistances, correspond- If the wave equation is multiplied by Ax and 
ing to a change in eigenvalue (energy level). the operator p= —1hd/dx is introduced, the equa- 
2. Although negative resistances are available x Bax , 
— for use with a network analyzer, in practice it is : 
more convenient to use a second type of circuit, 
the in which the positive and negative resistors are Tho 
1 to replaced by inductors and capacitors and the d.c. 
nts, currents and voltages are replaced by a.c. cur- arnt . | 
d be rents and voltages of fixed frequency. The use of 
the second type of interpretation is equivalent to _! 
multiplying the wave equation by 1=+/—1. The potential energy persion 
ntal In the diagrams to follow, unless otherwise 
arch stated, the inductors (whose reactance at the 
fixed frequency is denoted by Xx) may also be re 
" viewed as positive resistances of value X ;, and the | 
ys. 


capacitors (whose reactance is denoted by —X,) 
as negative resistances of value —X,. 
3. The third type of circuit contains inductors 
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The tote! energy operctor 


Fic. 1. Representation of energy operators by’admittances. 


‘ 
, 
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The Haomiltonion operator H* 


The operator H-E 
Fic. 2. Representation of the sum of operators by an 


interconnection. 
tion becomes 
— EAxy=0, 
where 


H=—+V=T+V. 
2m 


The three energy operators will be represented 
in the following manner: 1. The kinetic energy 
operator T= p/2m is represented (Fig. 1a) by a 
set of equal positive resistors (or inductive coils) 
in series extending from —« to +o. The 
impedance of each coil is (2m/h?)Ax. 2. The po- 
tential energy operator V is represented (Fig. 1b) 
by a set of isolated positive resistors (inductive 
coils). The admittance of each coil is VAx. These 
values vary from™point to point. 3. The total 
energy operator —E is represented (Fig. 1c) by a 
set of isolated equal negative resistors (capaci- 
tors). The admittance of each coil is the unknown 
quantity EAx. 

The sum of the two operators T and V forming 
the Hamiltonian operator H is represented by the 
interconnected network of Fig. 2a. That is, a 
summation of energy operators is represented 
electrically by an interconnection of the com- 
ponent networks. 

Finally, the sum of the Hamiltonian H and the 
total energy —£E operators is the resultant 
network of Fig. 2b extending from — © to +. 


OPERATION ON THE WAVE FUNCTION ¢ - 


The wave function y will be assumed to be 
represented (Fig. 3) by the differences of po- 
tential’! appearing between the junctions of the 

1It is possible to establish in one dimension a dual 
network in which y is represented by a current instead of 


a vol . However, in two and three dimensions the dual 
require ideal transformers. 


three types of coils and the ground connection 
(impedanceless wire). The function y varies 
along x, but is constant in time. An eigen- 
function y’ of H depends on the time as follows: 
exp [—1(E/h)t.] 

The differences of potentials appearing be- 
tween two junctions are Ay=(dy/dx)Ax+--- 
where higher order terms in the Taylor series 
development are neglected. 

The result of the operation ay or BAy (where a 
and £8 are operators) is represented by currents 
flowing in the respective component networks, as 
shown in Fig. 3. In particular: 1. The currents 
flowing in the capacitors are EAxj. 2. The cur- 
rents flowing in the vertical inductors are VAxy. 
3. The currents. flowing in the horizontal in- 
ductors are BAy = (h?/2m) (dp/dx). 4. The currents 
flowing out of the horizontal inductors at their 
junctions are A(BAy) = (h?/2m) (d*y/dx*) Ax, repre- 
senting (— TAxj). 

Since EAxy represents the currents flowing in 
the capacitors and HAxy=(T+ V)Axy those in 
the inductors, the equation HAxy = EAxy simply 
states Kirchhoff’s second law, that at each of the 
junctions of four coils the currents flowing into 
the positive resistors (inductors) are equal to 
the currents coming from the negative resistors 
(capacitors). 

Both the voltages and currents appear as 
standing waves. When the space distribution of 
y is such that it attenuates, the network may be 
terminated at any point by an equivalent 
impedance. When y does not attenuate the circuit 
may be terminated by a short circuit at any point 
of zero voltage, or by an open circuit at any point 
of zero current. These points are easily deter- 
mined on the a.c. network analyzer by trial. 

It should be noted that the equivalent circuit 
introduces an indeterminacy between quantities 
measured in the vertical coils and those in the 
horizontal coils. If, for instance, the vertical 
voltage y is known at a certain value of x, say xo, 
then the horizontal current (h?/2m) (dp/8x) is not 


Fic. 3. Representation of ¥ by voltages and ay by currents. 
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known at that particular xo, only at some inde- 
terminate value between x9 and x9+Ax/2 or 
xo —Ax/ 2. 


EIGENVALUES AND EIGENFUNCTIONS 


- Let it be assumed as an example that the V 
function is a potential well (Fig. 4a). Then in the 


' network the corresponding positive resistances 


assume either a constant or a zero value as shown 


_in Fig. 4b. In practice it is sufficient to extend the 


network to, say, twice the width of the potential 
well on either side. 

Let now a d.c. (or a.c.) generator be inserted 
anywhere in the network parallel with one of the 
negative resistances (inductors), as shown. If the 
values of all the negative resistances —EAx 
(capacitors) are simultaneously varied by the 
same amount, it will be found that the current 
(reactive current) in the generator varies and at 
some value of EAx becomes zero. 

It should be noted that while a current (reactive 
current) flows in the generator the circuit does 
not satisfy the differential equation, since at one 
point in the network (where the generator is con- 
nected) the currents do not add up to zero (as 
required by the equation) but to the generator 
current. That is, while a current flows in the 
generator, the voltages y do not represent solu- 
tions of the differential equation. Hence, only 
those network conditions are of interest in which 
the generator current is zero. 

Now a value E of the negative resistances, at 
which the generator current becomes zero, repre- 
sents a state at which the circuit is self-supporting 
and has a continuous existence of its own without 
the presence of the generator, as the negative 
resistances just supply the energy consumed by 
the positive resistances. (If the circuit contains 
inductors and.capacitors, the circuit is a resonant 
circuit and it oscillates at its basic frequency.) 
E is then an eigenvalue E,, while the voltage 
distribution across the capacitors (Fig. 4c) gives 
the corresponding eigenfunction yp. 

When the generator current is positive the 
circuit draws energy from the source, and when 
the current is negative the circuit pumps back 
energy into the source. At zero generator current 
the circuit neither gives nor takes energy, and 
theoretically the generator may be removed. All 
values of E at which the current crosses the axes 


and becomes zero are eigenvalues of the equation 
and the corresponding voltage distribution curves 
are eigenfunctions. When the energy level E 
overflows the well, the discrete spectrum of 
eigenvalues changes into a continuous spectrum 
and the generator current is zero at all greater 
values of E. 

When the energy E changes sign, the negative 
resistances become positive resistances and at no 
value of —E may the circuit be self-supporting 
(as it contains only positive resistances). That is, 
the equation has no negative eigenvalues. 


©) Eigen- function 


Fic. 4. Measurement of y. 


THE STATISTICAL MEAN OF OPERATORS 


To bring the measurements at the different 
energy levels to the same base, it is necessary to 
normalize the measured y values so that the new 
values of y satisfy the equation 


The measured y functions are normalized by 
plotting the square of y. If the area under the 
curve is 1/N?, all values of y are multiplied by NV. 
Then Ny is the normalized y. (Actually N may 
contain an arbitrary phase e*.) 
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Fic. 5. The variable-frequency model. 


The statistical mean of an operator a for a 
state y is defined as 


+0 


avya= f 


Since (aAx)y is a current (reactive current) 
flowing through an admittance, ¥*(aAx)y is the 
power (reactive power) in a single admittance. 
Hence the total power in a complete set of 
similar admittances, namely, 


+00 
Power => e*yve=> anh Ax, 


represents the statistical mean of the corre- 
sponding energy operator. That is: 1. The total 
power in all the vertical negative resistors is the 
average value of E. (That is, EZ itself, since y is an 
eigenfunction of H). 2. The total power in all the 
vertical resistors is the average value of the po- 
tential energy V. 3. The total power in all the 
horizontal resistors is the same as the total power 
in the vertical units, representing the average of 
the kinetic energy T= p?/2m. 

That is, the total power in all positive resistors 
is the same as the total power in all the negative 
resistors, or 


avyH = V=av,E=E. 
m 


THE THIRD MODEL 


Let the wave equation be divided by iw., where 
w.=+/w=(E/h)', and multiplied by Ax: 


1 dy VAx 
—— — —Ar'+ +hiadx 
2m Ax Ox? lide 
or 
Ax Ax 
—(H-E)y=0. 
We 


‘currents as in the first model. Their variation in 


— 


*fiax 


Fic. 6. Model of a free particle. 


In the present case: 

1. The kinetic energy operator T is represented 
(Fig. 5) by a set of equal inductors in series, 
whose inductance L; is (2m/h?)Ax. 2. The po- 
tential energy operator V is represented by a set 
of unequal coils in parallel, whose inductance L» 
is 1/VAx. 3. The total energy operator —E is 
represented by a set of equal capacitors in 
parallel whose capacitance is now Ax. (In 
the second model the capacitance was the 
unknown EAx.) 

The operand y is again represented by voltages 
and the result of the operation ay by the same 


time now is sinusoidal, with 2zf.=w.. 

Instead of varying the magnitude of the 
capacitors, now the frequency of the generator is 
varied, thereby varying the admittance of the 
capacitors, hw,.=E (and those of the inductors). 
Again when the generator current becomes zero 
the circuit is oscillatory and self-supporting and 
the network represents a stationary solution of 
the equation. The corresponding eigenvalue is 
E=hw=h(w,)?, rather than hw., because of the 
simultaneous variation of the reactance of the 
inductive coils. The eigenfunctions y of the model 
and of the equation are, however, identical. 

As the currents in.the horizontal inductors are 
(h?/2miw,)dy/dx, the results of an operation on 
by the momentum operator p=(h/i)0/dx are 
these currents divided by h/2mw,. Hence, in the 
third model the momentum operator p may be 
represented by a set of equal horizontal coils with 
inductance L = Ax/hw.. 

One slight disadvantage of this third model is 
that as the energy E changes sign, the reactance 
of the capacitor jw. cannot change signs. Since in 
most cases no eigenvalues exist in the negative 
energy range, this disadvantage is of little conse- 
quence. Of course, the second model with fixed 
frequency and variable capacitors works in all 
cases, since the capacitors simply become in- 
ductors when £ changes sign. 
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ELECTRIC CIRCUIT MODELS 


THE FREE PARTICLE IN ONE DIMENSION 


An interesting special case occurs when the 
potential V is zero everywhere. The one-di- 
mensional equivalent circuit of such a free 
particle is a conventional transmission line ex- 
tending to infinity in both directions (Fig. 6) in 
which the series inductance is 2mAx/h? and the 
shunt capacitor is hAx. 

It is well known that such a transmission line 
may maintain a standing wave at any frequency 
w. between zero and infinity drawing no current 
from the generator. That is, the positive energy 
values form a continuous spectrum. If the trans- 
mission line is considered as the second type of 
model with variable capacitors, then at negative 
energy values E the capacitors also become 
inductors and the line cannot maintain a standing 
wave. The corresponding free particle also has no 
eigenvalue at the negative energy levels. 


MODELS ALONG CURVILINEAR AXES 


The three-dimensional Schrédinger equation 
for a single particle is 


i? 
(V-E)y=0, 
2m 


where V? is the Laplacian operator in curvilinear 
coordinates. 

In order to establish a physical model for it, 
it is necessary to change it to a tensor density 
equation.2 The above equation in orthogonal 


*G. Kron, Proc. I. R. E. 32, 289-299 (May, 1944). 


Fic. 7. Model of the Schrédinger equation with three 
independent variables. 


curvilinear coordinates may be changed to a 
tensor density form by multiplying it by 
hyhoh3s= giving 
dy 
hy du! 


dy 
Ou? he du? 


oy 
V—E)y=0. 


If the equation is multiplied through by 
Au'Au*Au', it represents the surface integral of 
grad y around the six faces of a cube of space with 
volume hyh2hsAu'Au?Au*. The width Au* may be 
arbitrary and different in the three directions. 
The corresponding equivalent circuit is shown 
in Fig. 7. For a free particle (V=0) it represents 
a generalization of the conventional one-dimen- 
sional transmission line to three dimensions. 
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A.C. Network Analyzer Study of the Schrédinger Equation 


JANUARY 1 AND 15, 1945 


G. K. Cartert anp Gaprie, Kronft 
General Electric Company, Schenectady, New York 


(Received May 25, 1944) 


Utilizing the type of circuit described by Kron tests were carried out on the a.c. network 


analyzer. The tests reported include only one-dimensional circuits. Measurements were made of 
eigenvalues and eigenfunctions for the particular cases of the linear oscillator, the rectangular 
potential well, the double rectangular barrier, the single barrier, and the rigid rotator. These 
tests show the circuit representation to be valid and to be adaptable to solution on the network 


analyzer. 


ALTERNATING CURRENT NETWORK ANALYZER 


HE alternating current network analyzer! 
consists of a set of adjustable inductance, 
resistance, and capacitance units, each connected 
to a pair of flexible cords and plugs. Connections 
between units to form any desired network are 
made by inserting the plugs in adjacent jacks in 
a jack panel. As many units as desired can be 
connected to a common point. Alternating cur- 
rent electric power is supplied by a motor- 
generator set to individual generator units so that 
several different voltages, independently adjust- 
able in both phase and magnitude, can be in- 
_serted in different parts of the network. A 
centrally located set of instruments (voltmeter, 
ammeter, wattmeter, and varmeter) can be con- 
nected to any unit or circuit by a set of key 
switches. 


METHOD OF TEST 


The circuit used is that of Fig. 3 of the com- 
panion paper,’ repeated here as Fig. 1 for con- 
venience. The circuit was extended, not from 


4s 
Fic. 1. Equivalent circuit used. 


t Analytical Division, Central Station Engineering 
Divisions. 

tt Consulting Engineer. 

1H. P. Kuehni and R. G. Lorraine, Trans. A.I.E.E. 57, 


67 (1938). 
2G. Kron, Phys. Rev. 67, 39 (1945). 
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minus to plus infinity, but to the limit of the 
circuit elements available. 

Discrete energy levels are characterized in the 
circuit usually by a standing voltage wave which 
attenuates to zero outside of some special region. 


CURRENT 


GENERATOR 
$e¢s BB’ 


4 
+ 
| 
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Fic. 2. Generator a as a function of energy levels 


Fig. 3. 


For this type of wave, the circuit should be 
terminated in an equivalent impedance. Since no 
attempt was made to do this, such levels should 
show an error which increases with the voltage 
remaining at the end of the line. On the other 
hand, the continuous energy levels in general 
exhibit a standing wave which is a space sinusoid 
of constant amplitude extending to infinity in at 


least one direction. Hence, at the end of a finite 


length of line, such a wave must have a loop if the 
line is open-circuited, a node if it is short- 
circuited, or other space phase for other types of 
termination. Since only one length of line and one 
type of termination (open circuit) were used in 
these tests, only particular levels and only one 
wave at each level were actually measured. 

For the line elements, inductors and capacitors 
were used, the former being arbitrarily denoted 
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ANALYZER STUDY OF SCHRODINGER EQUATION 


vi 


Fic. 3. Eigenfunctions of the linear harmonic oscillator. 


as negative admittance units. Curves of generator 
reactive current versus energy level E (admit- 
tance yg) were measured for several generator 
locations along the line; from these the eigen- 
values (values of yg for zero reactive current) 
could be determined and the corresponding 
eigenfunctions (voltage distributions) measured. 
Since qualitative, rather than quantitative, 
results were the prime objective at this time, no 
special measures were taken to set admittance 
values closer than the nearest step on the 
analyzer units. ' 


LINEAR HARMONIC OSCILLATOR 


Taking Ax as unit distance,’ we represent the 
linear oscillator‘ by y,=h?/2m=1.10 and by 
yy = V=0.00704x*. The resulting generator re- 

* See Appendix for discussion of scale factors. 


_*V. Rojansky, Introductory ntum Mechanics (Pren- 
tice-Hall, Inc., New York, 1942). 


active current as a function of energy level yg=E 
is shown in Fig. 2. Ideally, the current curve 
shown in Fig. 2 should always slope upward 
toward the right (increasing values of 1/E), 
passing successively through zero, positive in- 
finity, and negative infinity ; the zero points then 
represent conditions in which the circuit should 
maintain its oscillations even without the ex- 
ternal driving force. Practically, the presence of 
resistance in the network forces the current to 
pass through additional zero points in going from 
large plus values to large minus values; but these 
additional zeros are readily recognized because 
the curve has the negative slope in passing 
through them, and they can therefore be ignored. 
The solid line is for one generator position; the 
dotted lines are check curves for other generator 
positions. Those eigenvalues at which the solid 
curve tries to reach zero but does not succeed 
correspond to eigenfunctions such that this 
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TABLE I. 


Measured 
eigenvalues 0.0921; 0.269; 0.442; 0.610; 0.775; 0.921 


Calculated 
eigenvalues 0.0880; 0.264; 0.440; 0.616; 0.792; 0.968 


TABLE II. 


0.0593; 0.216; 0.465; 0.800 


Measured eigenvalues 
0.0552; 0.220; 0.486; 0.832 


Calculated eigenvalues 


TABLE III. 
Measured | 4.14; 5.12; 6.10; 7.13; 8.12; 9.13 


generator position is near a node of the eigen- 
function ; it is not practical to supply the circuit 
losses from such a point. 

Figure 3 shows the measured eigenfunctions up 
to an energy level almost as great as the value of 
V at the end of the line (V=1.10). The measured 
eigenvalues compare with the calculated values 
as in Table I. The maximum discrepancy is only 
about 5 percent. At the lowest value, the dis- 
crepancy is largely caused by the approximate 
representation of the small values of V; while at 
the highest value, the line termination is proba- 
bly responsible. While the eigenfunctions have 
not been normalized for comparison, it can be 
seen that these too are in good qualitative agree- 
ment with the known forms for such functions.* 


POTENTIAL WELL 


The second case is that of a rectangular po- 
tential well. Taking again h?/2m=1.10 and Ax as 


G. K. CARTER AND G. KRON 


} E=0.216 


Fic. 4. An eigenfunction of the rectangular potential well 
in the discrete spectrum. 


* Reference 4, p. 24. 


unit distance, we represent the well by V=0 from 
—6 to +6 and by V)=1.10 outside this region. 
One of the eigenfunctions and eigenvalues in the 
discrete energy spectrum E< Vo, is shown in 
Fig. 4. Comparison of the measured with the 
calculated eigenvalues is given in Table II. The 
lowest figure in this case was not quite at zero 
current but represented the maximum inductance 
of the board units; the highest again is subject to 
some correction because of improper termination. 
Nevertheless, the agreement is very satisfactory. 
Again, though not normalized, the eigenfunctions 
show good agreement with the calculated shapes.* 

Figure 5 shows a measured eigenfunction in the 
continuous spectrum, E> Vo. This is, of course, 
only the particular function corresponding to a 
loop of the wave at each end of the existing line. 
For the same energy level, infinitely many other 
functions are theoretically possible; and a finite 


TRI 


E=1155 


Fic. 5. An eigenfunction of the rectangular potential well 
in the continuous spectrum. 


number of these, depending upon the number of 
circuit units per wave-length and the number of 
units available, could be obtained on the analyzer 
by simply shifting the nodes in space. Similar 
statements are true for any energy level E> Vo, 
of course. The figure is given here to show the 
performance of the circuit and the analyzer in the 
continuous spectrum region. 


DOUBLE BARRIER 


The double rectangular barrier was represented 
by Vo=1.10 between +5 and +10 and between 
—5 and —10, with V=0 elsewhere; Ax is again 
unit distance, and h?/2m = 1.10. Figure 6 gives six 


measured eigenfunctions corresponding to nodes 
at the ends of the line. The three curves on the 


*Reference 4, p. 156. 
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Fic. 6. Eigenfunctions for double rectangular potential barrier. 


left represent virtual binding.* The three on the 
right apparently represent a penetration phe- 
nomenon, related to the next section. The func- 
~ tion which is large on the outside of both barriers 
and small between was not obtained in this series 
of tests, because only one generator was used to 
supply resistance losses, and these are not readily 
transmitted through the low voltage region. The 
curves shown are again but individuals out of 
infinite families at every energy level. 


RECTANGULAR BARRIER 


The final test in this group was the single 
rectangular barrier represented by V=0 for x <0, 
Vo=1.10 between 0 and 3, and V,=0.40 for x>3; 
as before, Ax is unity and A?/2m=1.10. Figure 7 
shows some of the eigenfunctions measured. The 
three on the left represent penetration from the 
left, so to speak, at energy levels below Vi, be- 
tween V, and Vo, and above Vo; the three curves 


* Reference 4, p. 225. 


on the right, penetration from the right at levels 
above 


RIGID ROTATOR 


Legendre’s equation, which arises in the prob- 
lem of the rigid rotator,* can be written in the 
form 


—I(1+1) }@=0. 
~-101+1)| 


oa 

This is of the form of the Schrédinger equation, 
where h?/2m is replaced, by (1—x*), V by 
m?/(1—x*) and E by 1(/+-1); the range of x here 
is from —1 to +1. Tests on this equation were 
carried out only for m=0; one of the resulting 
measured @ functions with the corresponding 
value of / is shown in Fig. 8. Comparison with the 
theoretical values of /, in Table III, shows very 
good agreement in the range measured. The 
measured functions have not been normalized ; 


* Reference 4, p. 437. 
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Fic. 7. Eigenfunctions for rectangular barrier between potential levels: Left of barrier V=0; 


but calculated end points, indicated by the 
dotted end extensions, line up very well with the 
measured curves. 


SUMMARY 


The tests reported here demonstrate the appli- 
cability of the electric circuit equivalent for 
determining eigenvalues and eigenfunctions of 


one-dimensional equations of the Schrédinger . 


type. At the same time, they indicate that it is 
practical to use for this purpose existing a.c. 
network analyzers. Experience in these tests leads 
also to the conclusion. that the equivalent circuit 
together with a small acquaintance with trans- 
mission line theory goes far in aiding the visu- 
alization of the characteristics to be expected of 
the solutions. 

The significance of the equivalent-circuit 
method of studying characteristic values is that 
it is just as simple to study cases which are not 
readily solved analytically as it is to study the 
idealized cases given here. 


barrier height Vo=1.0; right of barrier Vi=0.40. 
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APPENDIX—SCALE FACTORS 


The equation represented in most of these 
tests is Schrédinger’s equation 


dy 
—-— —+[V(x)-E]y=0. 


The equation represented by the circuit of Fig. 1 
corresponds to 


d ty Je=0 


When the two equations represent the same 
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system, they give the proportions 
h?/2m V(x) E 


where K is the scale factor between the energy NL $2610 
level EZ and the admittance in the equivalent 
circuit. In Figs. 3-7, the unit of horizontal dis- 
tance has been taken as Ax; the physical distance Fic. 8. An eigenfunction of the rigid rotator. 
represented by this unit is determined by the 
equation 


h2/2m yy 
Ky, V(x) » 
In the case of the potential well, for example, 


(Ax)? = W/ = which is the case of Exercise 3.* Other proportions 

Vo for the well are obtained either by changing the 
: . number of Ax units in a or by changing the 
The half-width of the well is ratio of yy» to ya. 


(Ax)? 


a=6Ax * Reference 4, p. 155. 
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The Energy Spectrum of the Primary Cosmic Radiation 


KusaKa 
‘Smith College, Northampton, Massachusetts 


(Received October 11, 1944) 


ROM results of the experiments carried out 
in India, in Mexico, and the United States, 
Millikan, Neher, and Pickering! have put forward 
the hypothesis that the primary cosmic rays are 
electrons* whose energy arises from the complete 
transformation of the rest mass of atoms such 
as carbon, nitrogen, oxygen, and silicon which are 
abundant in interstellar space. Since the high 
altitude experiments of Schein, Jesse, and Wollan? 
indicate that the primary cosmic rays are not 
shower producing and hence are probably pro- 
tons, we have made a comparison with the ex- 
perimental results of Millikan and his colleagues 
and theoretical calculations based on the as- 
sumption that the distribution in number of the 
primaries is proportional to the inverse cube of 
the energy. The choice of the inverse cube is be- 
cause of the fact that the secondary particle, both 
the hard and the soft components, is known to 
have approximately this energy distribution. 
The main experimental results obtained by 
Millikan et al. are the variations with magnetic 
latitude of the energy brought in by cosmic rays 
arriving vertically, and the total energy brought 
in from all directions. It has been pointed out by 
these experimenters that the former is propor- 
tional to the area under the curve in which the 
rate of vertical coincidence is plotted against 
_altitude expressed in an equivalent height of 
water; and that the latter is proportional to the 
area under similar curves in which the ordinate 
is the rate of ion production in an ionization 
chamber. 
With our assumption of the energy spectrum 
and with the geomagnetic theory of cosmic radia- 


1 Millikan, Neher, and Pickering, Phys. Rev. 61, 397 
(1942); H. V. Neher and W. H. Pickering, Phys. Rev. 61, 
407 (1942); Millikan, Neher, and Pickering, Phys. Rev. 63, 
234 (1943). 

* Note added in proof: It has recently been pointed out 
by Warren [D. T. Warren, Phys. Rev. 66, 156 (1944) ] 
that the atom annihilation hypothesis is quite independent 
of the nature of the primary particles, so that they may 
be protons as well as electrons. 

* Schein, Jesse, and Wollan, Phys. Rev. 59, 615 (1941). 
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tion developed by Lemaitre and Vallarta,’ it is a 
simple matter to calculate the latitude variation 
of the vertical intensity which is inversely pro- 
portional to the minimum energy of arrival in the 
vertical direction, which has been done by these 
authors.‘ The calculation of the total intensity 
involves the integration of the product of the 
distribution in energy and the solid angle sub- 
tended by the main cone. The latter was obtained 
from the computations of Koenig® and the 
integration was carried out graphically. 

The results of the calculations and the experi- 
mental results of Millikan et al. are plotted in the 
accompanying figures. Since only the relative 
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VERTICAL INTENSITY 


va 


10 1s 20 
MAGNETIC LATITUDE 


Fic. 1. Vertical intensity as a function of magnetic 
latitude. The experimental any are for Bangalore (mag. 
ms nr Agra (mag. lat. 17° N), and Peshawar (mag. lat. 


values are obtainable in both the theoretical 
calculations and the observations, the theoretical 
values are normalized to give 100 at the equator; 
and the relative values of the experimental re- 
sults are adjusted to give best agreement with the 
theoretical curve. For the experimental points, 
the statistical error was estimated to be about 3 
percent and this is shown on the graphs. Figures 1 
and 2 show that, for the results on the vertical 


a os % Lemaitre and M. S. Vallarta, Phys. Rev. 50, 493 
‘In transforming the energy of the particles from stérmer 
units to electron volts, it is necessary to take definite values 
for the charge and mass of the primary particles. We have 
here assumed them to be protons, but it does not make any 
essential difference if we take them to be electrons. 
5H. P. Koenig, Phys. Rev. 58, 385 (1940). 
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SPECTRUM OF PRIMARY COSMIC RADIATION 
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TOTAL INTENSITY 
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Fic. 2. Total intensity as a function of magnetic lati- 
tude. en experimental points are for the same positions as 
in Fig. 1. 


and total intensity that were obtained in India, 
the agreement is within the experimental error. 
From Fig. 3 we see the values obtained in Mexico 
and in the United States. The points for Acapulco 
(magnetic latitude 25.6°), Valles (mag. lat. 31.1°), 


and Victoria (mag. lat. 32.8°) do not agree as - 


well, seeming to be just barely within the limits of 
statistical error. The points for Pasadena (mag. 
lat. 40.7°), St. George (mag. lat. 44.8°), and 
Pocatello (mag. lat. 51.0°) are completely off the 
theoretical curve. This divergence is probably 
caused by the fact that the spectrum of the 
primaries no longer follows the inverse cube law 
for energies of the order of 5 Bev, which is 
necessary for the vertical arrival of the primaries 
at Pasadena. 

Another point which has been put forward as 
evidence of the discontinuous nature of the pri- 
mary spectrum is that the observed vertical 
intensity at Acapulco, Mexico (mag. lat. 25.6°) is 
8.3 percent lower than that at Peshawar, India 
(mag. lat. 24.8°). On the atom annihilation 
hypothesis, the intensities at both places should 
be the same, while the continuous distribution 
hypothesis holds that the intensity at Acapulco 
should be higher. The observed sea level longi- 
tude effect gives a maximum intensity off of the 
South American coast and a minimum west of the 
Malay Peninsula with amplitude variation of 
about 4 percent. However, it is to be noted that 


VERTICAL INTENSITY 


10 20. + «30 
MAGNETIC LATITUDE 


Fic. 3. Vertical intensity as a function of magnetic 
latitude. The ee eg points are for Acapulco (mag. 
lat. 25.6° N), Valles (mag. lat. 31.1° N), Victoria (mag. lat. 
32.8° N), Pasadena (mag. lat. 40.7° N), St. George (mag. 
lat. 44.8° N), and Pocatello (mag. lat. 51.0° N). 


the theoretical calculations of Vallarta® give a 
maximum off of the west coast and a minimum 
in the Pacific Ocean. From the theoretical stand- 
point, there should be no difference in intensity 
at Acapulco and Peshawar because of the longi- 
tude since they are both just about midpoint 
between the longitudes at which the maximum 
and minimum intensities occur. It seems more 
reasonable to take the theoretical result rather 
than the observed sea level longitude effect since 
the discrepancy in phase possibly may be caused 
by local surface effects which do not influence the 
cosmic rays arriving at the top of the atmosphere. 
The fact that Acapulco is 0.8 degree farther 
north than Peshawar should make the intensity 
about. 2.0 percent higher at Acapulco; but this 
percentage is within the statistical error. 

Thus, generally, we may conclude that the 
agreement between the theoretical calculations 
and the experimental results is within the experi- 
mental error if the distribution is taken to be a 
curve which follows the inverse cube law for high 
energies and which flattens out on the low energy 
side at about 6 Bev. This interpretation, then, is 
offered as a possible alternate to the atom 
annihilation hypothesis derived from the experi- 
mental results of Millikan and his colleagues. 


*M. S. Vallarta, Phys. Rev. 47, 647 (1935). 
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PHYSICAL REVIEW 


Letters to the Editor 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 


Evidence for the Artificial Production of a New 
Neutral Radiation 


GERHART GROETZINGER, P. GERALD KRUGER, AND Luoyp SMITH 
The Nuclear Radiations Labor . Department of Physics, 
University of Ilinois, Urbana, Illinois 


January 3, 1945 


ROMPTED by indications, found by one of us, of a 
new non-ionizing radiation produced by a cyclotron, 
which manifested itself by the appearance of high energy 
electrons, we have performed a series of experiments to 
establish its existence definitely. 

Since gamma-rays and neutrons (which are produced in 
large numbers in a cyclotron) may also give rise to electrons 
of several million volts energy, it was necessary to devise 
methods of detection which could distinguish between these 
various radiations. 

Two Geiger-Mueller counters were set up at a distance 
of four meters from the target of the cyclotron and almost 
in direct line with the deuteron beam hitting the target. 
The axes of the counters were five cm apart, the line joining 
them being perpendicular to the direction of the deuteron 
beam. The counters had a wall thickness so that only 
electrons of energies higher than 2 Mev have an appreciable 
chance to trip both counters in coincidence. The cyclotron 
was completely enclosed by water tanks of 120-cm thick- 
ness. The tanks on the side facing the counters contained a 
saturated solution of boric acid. A lead house of 19-cm 
minimum thickness surrounded the counters on all sides. In 
addition, this cubical lead house was surrounded on the 
. sides by 0.5-mm sheets of cadmium and by boric acid solu- 
tion in water, at least 23 cm thick. The top and bottom of 
the house were shielded by Cd sheets and 16 cm of paraffin. 

With this arrangement we recorded the number of co- 
incidences and simultaneously the total number of 10.5- 
Mev deuterons striking the target of the cyclotron. The 
beam current was of the order of 30-50 microamperes and 
the number of coincidence counts about twice background. 
Table I shows the results for Mn and W targets. 

In order to establish how many of the coincidences were 
caused by fast and slow neutrons, the water cans and Cd 
shields between the counters and the target were removed. 
This should have increased the neutron intensity by a 
factor of 100 inside the lead house, but no measurable effect 
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on the counting rate was observed. Even the 1500 Ib. of 
boric acid which were dissolved in the 120-cm thick water 
tanks did not affect the counting rate in the case of the W 
target. This experiment, and further work to be described 
in a more complete report which will appear later, show 
that the observed counting rates are not the result of 
neutrons. 

The following experiments were performed to determine 
to what extent these coincidences could be caused by 
gamma-rays. Removing 9.5 cm of lead from the top and 
sides of the lead house increased the counting rate by a 
factor of 25 in the case of the W target, and by a factor of 
22 for the Mn target. Since for gamma-rays of maximum 
penetrating power in lead this factor should be 80, it is 
concluded that a radiation more penetrating than gamma- 
rays exists. In the case of the Mn target, an additional layer 
of lead, 9.5 cm thick, placed between the target and the 
counters reduced the counting rate only by a factor of 4. 
In this connection it should be mentioned that-.separate 
experiments showed that the counting rate was not affected 
by removing 9.5 cm of lead from the other three sides of the 
house, indicating that all of the coincidences are caused by 
radiation arriving from the direction of the cyclotron. 

Second, the counter telescope was turned through 90 
degrees, so that the line joining the axes of the counters was 


TABLE I. Coincidences per deuteron ( X10"). 


Lead absorber 
9.5 cm 19 cm 


Target Neutron 

material shields 28.5 cm 
Ww in 37.7+1.2 1.5+0.2 
Ww out 
Mn in 46.0+1.7 2.1403 0.53 +0.28 
Mn out 2.2+0.6 


1 Marcel Schein (private communication). 


-parallel to the deuteron beam. Under these conditions, the 


counting rate was unchanged for lead thicknesses of 19 cm 
and 28.5 cm, while it was increased by a factor of 2.5 for 
the 9.5-cm thickness. This increase in the latter case indi- 
cates that most of the radiation penetrating 9.5 cm of lead 
still consists of gamma-rays coming from the cyclotron; the 
isotropism of the electrons recorded in the former cases 
indicates that a different radiation is producing them. 
Other measurements indicated that the energies of these 
last-mentioned electrons are higher than 5 Mev. 

Similar experiments were performed with P, Sb, and 
NaCl targets and it was found that the yield does not 
change appreciably with target material. 

From these data it appears that the radiation here ob- 
served is more penetrating than neutrons in water or 
gamma-rays in lead. Therefore it seems likely that the 
radiation consists of neutral particles of low mass (low 
mass mesotrons). In this case the observed electrons would 
be decay electrons which are emitted after the mesotrons 
are slowed down or stopped. It is believed that there is 
evidence for charged mesotrons of low mass in cosmic 
radiation.! 
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_ Absorption of 2.8-Mev Gamma-Rays in Lead» 


GERHART GROETZINGER AND LLOYD SMITH 
Department of Physics, University of Illinois, Urbana, Illinois 
January 3, 1945 


N a recent paper! an absorption coefficient of 0.405 
for 2.85-Mev gamma-rays in lead was reported. This 
value, taken together with the result of a similar measure- 
ment in copper, was interpreted as indicating an “‘insuffi- 
ciency” in the Klein-Nishina formula for this energy, since 
the absorption coefficient predicted by the theory (which 
also takes into account absorption caused by the photo- 
electric effect and pair production) is 0.46 cm=2 On the 
other hand, complete agreement was found at 1.30 and 1.14 
Mev. The discrepancy of more than 10 percent at 2.85 Mev 
is the more striking since the theory has been confirmed at 
higher and lower energies by many workers. 

The reported low absorption coefficients were obtained 
by measuring the ionization produced by the gamma-rays 
of radioactive Na after passing through various thicknesses 
of lead and copper. Since this radiation is assumed to con- 
tain both 1.38- and 2.85-Mev gamma-rays,’ the absorption 
coefficient of the high energy component had to be deduced 
from the absorption of the mixture. 

To provide further information on this matter we have 
measured the absorption of Na gamma-rays in lead by an 
arrangement which, for two reasons, practically excluded 
the effect of the low energy component. First, we used a 
minimum absorber thickness of 9.5 cm, so that the ratio of 
the high energy to the low energy quanta is raised to 3:1 
according to the theory or even higher according to the 
above-mentioned paper;! second, we recorded the gamma- 
radiation with a twofold Geiger-Mueller counter arrange- 
ment of high resolving power, the counter walls being of 
such a thickness that only electrons of energies greater than 
2 Mev had an appreciable chance of penetrating both 
counters and being recorded. 

On account of a strong activity in the room in which the 
experiment was performed‘ the counters were surrounded 
on all sides by at least 19 cm of lead, except for an opening 
on the side facing the Na source in which lead was placed 
to make up a number of absorber thicknesses between 9.5 
and 19 cm. The source was placed 40 cm from the first 
counter, a position which allowed a suitable counting rate 
even at the greatest thickness. The effect of the imperfect 
geometry of this arrangement was overcome by pointing 
the counter telescope toward the source, causing it to be 
insensitive to radiation proceeding in any other direction 
than straight from the source. This directional effect occurs 
because most of the electrons which are produced by 2.8- 
Mev gamma-rays and which have sufficient energy to pene- 
trate both counters are Compton electrons emitted in a 
forward direction. This was proved by turning the counter 
telescope through 90°, which reduced the counting rate by 
a factor of 3. 

Using lead absorber thicknesses of 9.50, 12.13, 14.60, and 
19.0 cm we found an absorption coefficient of 0.467+0.009 
cm~, The absorption curve is of pure exponential form, 
indicating that the 1.98-Mev radiation was indeed elimi 
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nated by the method used. The result thus substantiates 
the present theory. 


5 Cork and R. W. Pidd, Phys. Rev. 66, 227 fyi). 
Heitler, The Quantum Theory of Radiation (Clarendon Press, 
Oxford, 1936). 

3 Actually ents by P. Gerald Kruger and W. E. 
that there are four gamma-ray lines between 2.6 and 2.9 
in this ion the absorption ient varies very slowly with energy, 
this fact is of minor importance for absorption measurements (private 
communication). 

4 The arrangement was primarily set up for a different purpose. 
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The Nuclear Isomerism of Gold 


MARCELLUsS L. WIEDENBECK 
Depariment of Physics, University of Notre Dame, Notre Dame, Indiana 
November 27, 1944 


TRONG activity has been observed when gold is 
irradiated with high energy x-rays. This activity has 
been observed when self-quenching counters, having gold 
cathodes, have been irradiated directly; and it can also be 
detected when a gold sheet is irradiated and quickly placed 
before a counter with a thin (1.5 mil) aluminum window. 

This activity is attributed to a metastable state of stable 
zsAu'*? produced by line absorption in the following 
reaction: 

The period, evaluated from a decay curve extending over 
one minute, has been found to be 7.50.5 second and the 
energy of the metastable level, which has been evaluated 
from the measurement of the energy of the conversion 
electron by absorption in aluminum foils, is 250 kev. 

Since the spin of gold in the ground state is 3/2, the spin 
of Au* is probably 9/2 or 11/2. 

It is of interest to note that gold, along with platinum!' 
and mercury,? form a triad similar to the one formed by 
silver, cadmium, and indium. 

A detailed study of the energy levels which combine with 
the metastable state is in progress and will be reported in 
the near future. 

Friedlander and Wu, P 


Kamen, and Ruben, Phys. Rev. 1937); 
Thornton, Phys. Rev. 52, 239 (1937). 


“Pool, Cork. and 


Gamma-Ray Absorption 
J. M. Corx 


University of Michigan, Ann Arbor, Michigan 
December 23, 1944 


N a recent paper! the linear absorption coefficient for 
certain gamma-rays in lead and copper was reported. 
Attention was given particularly to the 2.75-Mev radiation 
from sodium (24), showing an apparent discrepancy between 
observed and computed values of the absorption coeffi- 
cients in lead and copper. To examine more conclusively the 
behavior of the other radiations mentioned, namely, that 
from zinc (65) and from cobalt( 60), much stronger sources 
of these gamma-rays were studied. 

A sample of radio-cobalt of 32 millicuries equivalent 
intensity was kindly lent by Professor R. D. Evans of the 
Massachusetts Institute of Technology and a sample of 
zinc was bombarded as a probe in the cyclotron in this 
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TABLE I, Linear absorption coefficients in cm=!. 


Absorber 1.14 Mev 1.30 Mev 
Iron 0.426 0.400 
0. 0.440 
Lea 0.678 0.610 


laboratory for many hours. With these very strong samples 
it was possible to make observations through much thicker 
layers of the absorber than had been done before and hence 
more accurately determine the absorption coefficients. 
These values together with the absorption coefficients for 
iron are collected in Table I. 

The same values are shown graphically, along with the 
curves computed? by Heitler, in Fig. 1. There is some 
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Fic. 1. Absorption coefficients, observed and computed (Heitler). 


uncertainty in the energy of the cobalt radiation as it has 
been reported* to be complex, with energies 1.1 and 1.33 
Mev. The value of the absorption coefficient recorded is 
that after passing through 9 cm of lead, and hence 1.3 Mev 
is used. The assumption of a lower value for the energy 
would have increased the observed discrepancy. 

It is apparent that in every case the experimental points 
are below the computed Heitler curves, leading to a con- 
clusion similar to that expressed in the previous paper, 
namely, an inadequacy in the Klein-Nishina formula at 
high energies. 

M. Cork and R. W. Pidd, Phys. Rev. 66, 227 (1944). 

: Heitler, The Quantum Theory of Radiation (Clarendon Press, 

L. G. Elliott, Phys. Rev. 62, (1942); 


. Deutsch and L 
Mandeville and H. W. Fulbright, Phys. Rev. 64, 265 


A New Method for Determining Thresholds in 
y-n Processes 


MARCELLUS L. WIEDENBECK AND C. J. MARHOEFER 
Department of Physics, University of Notre Dame, Notre Dame, Indiana 
October 20, 1944 


-g N exact determination of the threshold energy of 
A photo-disintegration is of importance in nuclear 
physics and is especially useful in determining the relative 
masses of various nuclei. Up to this time, only two elements 
have been quantitatively studied in this connection, 
namely, beryllium and deuterium. However, it seems 
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probable that many more will be studied in the near future, 
and a discussion of the methods previously used and of the 
new method developed in this laboratory thus seems 
justified. 

The earliest method used! consisted in irradiating a 
sample of beryllium or deuterium with the y-rays from 
natural sources, ThC”, etc., the energy of which is known, 
and observing the energy of the proton liberated from the 
process, ,D?++—>,H!+on'. The threshold energy was then 
the energy of the y-ray minus twice the energy of the 
proton. The values obtained for the threshold determined 
in this manner vary from 2.16 to 2.26 Mev. More recent 
experiments indicate that the threshold for deuterium is 
2.18 Mev and that of beryllium is 1.63 Mev. 

With the introduction of various devices to produce high 
energy electrons, it became possible to vary easily the 
maximum energy of the x-rays and thus to study the cross 
section for disintegration as a function of the accelerating 
potential. In general, by this method the emitted neutrons 
were slowed down with paraffin and then detected by means 
of the activity produced in silver or rhodium. The activity 
obtained from such an arrangement is proportional to the 
total cross section for the process and when plotted against 
the accelerating voltage, all other quantities kept constant, 
gives a curve? of the form A =K(V—V:)? where K and p 
are experimental constants and V; is the threshold po- 
tential. Such a curve for deuterium is plotted in Fig. 1 
(broken curve). The value of the threshold can be obtained 
by extrapolating this curve to zero activity. 

The method described above has the disadvantage that 
the value chosen as the threshold depends only upon the 
few experimental points near the threshold where the 
activity is low. For a higher degree of accuracy, it is 
imperative that several points be obtained very close to the 
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Fic. 1. Neutron counting rate versus accelerating potential. Broken 
curve, } t.... from experimental points marked as crosses, obtained by 
in between deuterium ‘nen and detector. Solid 


ter, 
curve without the use of 


The present method makes use of the resonance absorp- 
tion of neutrons by rhodium or silver. It is well known, 
from theory* and from experiments on the excitation of 
nuclei by x-rays,** that the intensity of any isochromat in 
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the thick target continuous x-ray spectrum increases 
linearly with applied electron accelerating potential greater 
than the energy of the isochromat. Therefore, in a 7- 
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Fic. 2. Neutron counting, rate versus accelerating eye for 
disintegration of beryllium. The curve was obtained without the use 


of paraffin. 


n process, the number of neutrons of a given energy should 
also be linear with the applied voltage V, when 


A 
v= Va 


(V; being the threshold potential; A, the atomic weight of 
the nucleus being disintegrated; and V,, the energy of the 
neutron being considered). 

If no paraffin or other hydrogen-containing materials is 
present to slow down the emitted neutrons, the detector 
will be sensitive to essentially only one such neutron “‘line”’ 
of energy equal to the resonance energy of the detector and 
will not be affected by the faster neutrons. Thus the 
activity vs. accelerating potential curve should give a 
straight line intersecting the abscissa at the threshold 
potential. 

This conclusion was tested with both deuterium and 
beryllium. The detector used was an argon-ether filled 
counter with a rhodium cathode. Small samples of deuterium 
and beryllium were bombarded for two minutes by the 
x-rays produced by a beam current of 100 microamperes 
striking a thick gold target. The activity was taken as the 
number of counts above background obtained during the 
two minutes after the irradiation was stopped. 

The activity is plotted (solid curves in Fig. 1 and Fig. 2) 
from the threshold to 3.2 Mev as a function of the applied 
potential. It is seen that in both cases a straight line is ob- 
tained, which when extrapolated to zero activity, gives the 
threshold for the process, namely, 1.630+0.006 for beryl- 
lium and 2.185+0.006 for deuterium. 

Thus, the threshold can be determined from a linear 
curve drawn through many points separated by very con- 
siderable distances. By this method the thresholds can be 
determined with a high degree of accuracy. 

N, Feather, Nature 13, $63 (1935): Proc. Roy. Soc. 151, 479 (1935). 
Guth: Phys. Rev, 99. 325 (1941) 
iedenbeck, Bull. Am. Phys. Soc. 17, Nos. 


*B. Waldman and M. Wi 
4, 5 (1942). 
*M. Wiedenbeck, Bull. Am. Phys. Soc. 19, Nos. 3, 5 (1944). 
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Showers of Penetrating Particles 


O. SALA AND G. WATAGHIN 
Department of Physics, University of Sio Paulo, Sio Paulo, Brazil 
November 24, 1944 


OMPARATIVE studies of showers of penetrating 
particles at various altitudes with different materials 
were begun in July and August, 1944, in Campos de Jordao 
(Brasil) at an altitude of 1750 m and latitude 23°, and in 
Sao Paulo (altitude 750 m and latitude 23° 5’). 

The multivibrator circuit,! elaborated by M. D. de Souza 
Santos, was adopted in connection with counters of fast 
type (alcohol-argon mixture). The resolving time of the 
coincidence set used was ~6X 10 sec. Tlie rate of chance 
4-fold coincidences in all experiments was negligible. The 
efficiency was tested before and after the experiments and 
was >97 percent. The arrangements XV and XVI are 
indicated in Fig. 1. They are of the type used in previous 
experiments in order to observe showers of penetrating 
particles.* 

Fourfold coincidences were observed between counters 
fully surrounded by lead of thicknesses not smaller than 
10 cm Pb and separated also by 10 cm of lead. 

In the experiment XVI the arrangement of counters and 
lead was the same as the experiment XV, and only an 
absorber of water was added. A total amount of about 750 


40+ 


40 


xv 


liters of water was used forming a layer of 80-cm thickness. 
The preliminary results are given in Table I. 
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Although the number of observed 4-fold coincidences is 
very reduced it seems to us that the results can be con- 
sidered as an indication of an increase of 23 percent in the 
frequency of penetrating showers, because of the additional 
layer of water. Indeed the observed increase is more than 3 
times the standard error. Obviously this layer functions as 


TaBLe I, Number of 4-fold coincidences in arrangements XV and XVI. 


4-fold time freq. 
arrang. coinc. min. 10° min.“ 
S. Paulo XV 73 36300 2.0+0.3 
C. de Jordio XV 155 21830 7.120.5 
C. de Jordio XVI 290 32440 8.9+0.5 


an absorber and as a source of secondary radiation. Our 
observations seem to indicate that groups of particles 


of water of only 80 cm. 
mn Recently V. Regener* and Schein, Iona, and Tabin‘ ob- 
_ served production of particles in paraffin. Janossy and 
Rochester observed production of groups of penetrating 
secondaries in lead. Evidence on production of groups of 
penetrating particles in heavy materials was obtained also 
in cloud-chamber photographs.*® 
The comparison of frequencies in the arrangement XV 
at two altitudes (1750-750) gives an increase by a factor 
-~3.5. Thus the shower producing radiation is rapidly 
absorbed in a layer of 1000 m of air (at an atmospheric 
pressure ~66 cm Hg, the absorbing mass being ~100 
g/cm?*). This finding and the fact that shower producing 
radiation generates showers of penetrating particles in layer 
of water ~80 g/cm? seems to indicate that this radiation 
could not be responsible for showers of penetrating particles 
under clay in a depth of 50 m water equivalent’ (observed 
in S. Paulo). It seems possible that there exist several types 
of showers of penetrating particles produced by different 
kinds of rays. Further studies are in progress. 
1M. D. de Souza Santos, An. Acad. Bras. de Ciéncias 12, 183 (1940). 


2 Wataghin, de Souza Santos, and Pompeia, Phys. Rev. 57, 61 (1940); 
57, ge (1940); 59, 902 (1941). 

V. H. Regener, Phys. Rev. 64, 250 (1943). 

‘Scheie. Iona, and Tabin, Phys. Rev. 64, 253 (1943). 

5 L. Janossy and G. D. Rochester, Nature 150, 633 (1942). 

* E. O. Wollan, Symposium sobre Raios Cosmicos, Acad. Bras. de 
Ciéncias, 123 (1941). D. Hughes, ibid. 67 (1941). W. E. Hazen, Phys. 
Rev. 65, 67 (1944). 
we sobre Raios Cosmicos, Acad. Bras. de Ciéncias, 155 
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penetrating more than 30 cm of Pb are produced in a layer - 
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The Two-Body Problem in Birkhoff’s and 
Einstein’s Theories 
CARLETON W. BERENDA 
Physics Department, Massachusetts State College, Amherst, Massachusetts 
November 17, 1944 

N a recent article by the late G. Birkhoff and others,! 

a new general theory of relativity (with a flat space- 
time) was further advanced.? The authors contend, in sup- 
port of their theory, that it provides a solution of the two- 
body problem while Einstein’s theory does not. However, 
Einstein, Infeld, Hoffmann, and Robertson in 1939 and in 
1940? submitted a solution of the two-body problem which, 
to the approximation desired, predicts a rate of advance of 
periastron per revolution of 


P.= 


Relativistic units of mass (cm) are used. 
Birkhoff’s theory predicts, to his approximation, a value 
of 


| Qn 
a(1—e*)° 


This gives 


mim: 
P. 3 (mi+m2)” 


or, for m:=mz (the optimum test case): 
P,—P.=P./12. 


Hence, the maximum difference obtainable between the 
two theories certainly falls within the limits of experimental 
error, and may arise from differences in the approximation 
methods employed by the two theories. If any preference 
can be given to Birkhoff's theory, it must be on the ground 
of simplicity of its general mathematical (geometrical) 
structure. All other predictions are the same as in Einstein's 
theory. The special assumption, in Birkhoff’s theory, of a 
“disturbance velocity of matter equal to c”’ for subatomic 
particles, seems a reasonable postulate in the light of the 
equations of retarded potentials in Maxwell's electrody- 
namics. The assumption of flat (homaloidal) space-time in 
a gravitational field is as acceptable per se as it is in N. 
Rosen's‘ or in A. N. Whitehead’s® theory of relativity. 

1 Barajas, Birkhoff, Graef, and Vallarta, Phys. Rev. 66, 138-143 
Nat. Acad. Sci. 29, 231 (1943). 


3 Einstein, Infeld, Hoffmann, and Robertson, Ann. Math. [1] 39, 65, 


101; [2] 41, 455. 
+N. Rosen. Phys. Rev. 57, 147-155 (1940). 
SA. N. Whitehead, Principle of Relativity. 


fc 
h 
al 
di 
R 
sp 
ar 
ele 
M 
fer 
da 
ap 
sic 
G. 
to 
fur 


- 


TY 


ue 


PHYSICAL REVIEW 


VOLUME 67, NUMBERS 1 AND 2 


JANUARY 1 AND 15, 1945 


Proceedings of the American Physical Society 


MEETING AT CHICAGO, DECEMBER 1 AND 2, 1944 


HE 263rd meeting of the American Physical 
Society was held in the Museum of Science 
and Industry at Chicago on Friday and Satur- 
day, the first and second of December, 1944. The 
registration came to one hundred and twenty-six. 
President Dempster assumed the functions of the 
Local Committee; the gratitude of the Society is 
due him for the arrangements and to the Museum 
of Science and Industry for the excellent facilities 
which they put at our disposal without charge. 
The meeting was distinguished by a symposium 
of fourteen invited papers on spectroscopy. The 
Secretary acknowledges with thanks the aid of 
R. S. Mulliken, H. H. Nielsen, and O. Struve in 
seeking and finding speakers able in these difficult 
times to illustrate so many of the facets of 
spectroscopy, and expresses the thanks of the 
Society to the speakers themselves. Their names 
and the titles of their papers are appended to 
these minutes. At one of the sessions opportunity 
was provided for the planning of the proposed 
Division of Spectroscopy, and many and valuable 
were the opinions and the suggestions which were 
received. 

Twenty-six contributed papers were submitted 
for the meetings; the abstracts are printed 
hereinafter. 

The dinner of the Society was held at the 
Windermere West Hotel on Friday evening with 
an attendance of almost one hundred. The after- 
dinner speakers were P. Klopsteg, J. Franck, 
R.S. Mulliken, and H. H. Nielsen, Mr. Klopsteg 
speaking of the scientific institutions of Australia 
and the other three on the topic of the Division of 
Spectroscopy. 

The Council met on Friday afternoon, and 
elected to Fellowship two candidates and to 
Membership the unprecedented number of no 
fewer than one hundred and eighty-two candi- 
dates; their names follow. The Council further 
appointed an Organizing Committee for the Divi- 
sion of Spectroscopy composed of K. K. Darrow, 
G. R. Harrison, R. S. Mulliken, and three more 
to be selected by those three. The Council 
further authorized, in response to numerous peti- 


tions, the formation of a Division of Chemical 
Physics in the Society, and appointed as its 
Organizing Committee the following: O. Beeck, 
F. G. Brickwedde, K. K. Darrow, J. E. Mayer, 
R. S. Mulliken, and O. K. Rice. 

The Society has lost through death the Fellows 
and Members here listed: John M. Adams, H. A. 
Clark, Milo A. Durand, Sir Ralph Fowler, and 
F. W. Stallmann. 

Advanced from Membership to Fellowship: M. L. 
Huggins. 

Elected to Fellowship: W. T. Szymanowski. 

Elected to Membership: A. R. Anderson, D. L. 
Anderson, D. E. Andrews, Jr., John Angus, D. A. 
Baldwin, T. F. Ball, N. H. Barbre, J. S. Barlow, 
Stanley Bashkin, J. R. Beatty, Richard Bellman, 
A. L. Bennett, Theodore Berlin, Marietta Blau, 
S. D. Bloom, J. C. Boonshaft, I. B. Born, John 
Bowman, L. L. Boyarsky, I. R. Brenholdt, W. G. 
Brock, C. J. Burbank, E. A. Burrill, F.X. Byrnes, 
Leon Camp, R. R. Carhart, C. W. Christensen, 
R. J. Christensen, H. H. Q. Chun, W. H. 
Clohessy, B. L. Cohen, R. E. Corby, N. H. Coy, 
W. S. Cramer, Lawrence Cranberg, J. A. Craw- 
ford, L. L. Cruise, O. R. Cruzon, J. R. Curry, H. 
W. Curtis, Jack Dainty, E. B. Dale, S. L. Dart, 
J. R. Davis, E. S. Dayhoff, W. B. Denny, P. M. 
Doty, W. W. Drake, L. F. Drummeter, G. E. 
Dewell, B. G. Eaton, C. J. Engle, L. C. Erich, 
W. F. Fairbank, Abraham Fineman, M. D. Fiske, 
E. R. Fitzgerald, Harold Fleisher, E. P. Fricke, 
M. M. Freundlich, J. A. Galvin, W. C. Geer, M. 
I. Goldberg, Oscar Goldman, J. E. Goldman, 
P. C. Goldmark, L. S. Goodman, Kathryn 
Greene, Jack Greenfield, R. M. Hadley, W. C. 
Hall, E. A. Hamacher, W. A. Hargreaves, S. P. 
Harris, C. W. Harrison, J. R. Holmes, Leroy 
Hughbanks, H. W. Hunter, E. C. Y. Inn, B. A. 
Jacobsohn, A. T. Jaques, D. G. Jelatis, R. R. 
Johns, H. C. Johnson, Christian Johnston, J. F. 
Kane, Rockwell Kent III, M. J. Klein, R. W. 
Koza, Arnold Kramish, G. L. Krieger, F. W. 
Kuether, W. A. Ladd, J. R. Laugham, Martha 
H. Lawson, Boris Leaf, Irvin Levin, Harold 
Levine, J. W. Lewis, W. F. Libby, P. J. Loatman, 
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M. J. Lun, W. B. Lurie, W. F. Main, E. A. Marx, 
J. E. Mayer, W. S. McAfee, F. T. McClure, W. J, 
McGonnagle, E. S. McKee, K. A. McMurtrio, 
Edward Melkonian, R. P. Miller, William Miller, 
D.J. Montgomery, J. R. Morgan, F.S. Mortimer, 
Nathan Most, H. R. Molton, R. E. Mueser, H. 
M. Muncheryan, J. D. Orndoff, Eleanor Oshry, 
Alfred Owyang, Murray Peshkin, R. P. Peterson, 
R. L. Platzman, R. V. Pound, W. W. Pratt, 
William Primak, Wilfrid Rall, D. H. Ransom, R. 
R. Rau, W. C. Redman, R. B. Reddy, G. W. 
Renner, R. C. Reed, F. G. Rest, D. P. Riley, H. 
L. Ritter, F. L. Roth, Paul Rudnick, Stanley 
Ruthberg, D. S. Sakon, A. L. Schawlow, R. E. 


Schofield, A. H. Sharbaugh, M. J. Sheehy, W. A. 
Sisson, V. D. Snyder, K. H. Sommermeyer, 
Anthony Sperduto, R. E. Stephens, D. W. 
Stewart, J. L. Stewart, W. G. Stroud, Mrs. 
George Sutton, Preston Taulbee, E. A. Taylor, 
R. W. Treharne, P. S. Turner, Douglas Venable, 
G. C. Wallick, H. F. Webster, R. A. Weiss, A. M. 
Weitzenhoffer, G. F. Welch, R. M. Whaley, D. 
H. Wilkinson, F. H. Willis, J. M. Willson, W. F. 
Witzig, Victor Wouk, Dorothy Wrinch, K. K. 
Wyckoff, Sidney Visner, T. F. Young, L. W. 
Zabel, M. G. Zabetakis, Henry Zatzkis, Paul 
Zilezer, Kurt Zuber. 
Kart K. Darrow, Secretary 


SYMPOSIUM ON SPECTROSCOPY 


Address of Welcome by K. K. Darrow. 


Forbidden Lines in the Laboratory: the Past and the Future. S. Mrozowsk1, University of Chicago. 

Partial Selection Rule for Sensitized Fluorescence. J. G. WiINANS, University of Wisconsin. 

Some Features of the Spectroscopy of the Atmosphere. O. R. WuLr, University of Chicago. 

Double-Bond Absorption in the Vacuum Ultraviolet. J. R. PLatt, Northwestern University. 

The Near Ultraviolet Absorption of Pyridine Vapor. H. Sponer, Duke University. 

Spectroscopic Properties of Chlorophyl and Porphyrines. E. Rasinowitscu, Massachusetts Institute of Technology. 
Measurements of Line Width in the Infra-Red Spectra of CO2, N.O, and NH;. A Ape Anp E. F. BARKER, University 


of Michigan. 


The Infra-Red Absorption Spectrum of Deutero-Methyl Chloride. A. H. Nie_sen, University of Tennessee, AND H. H. 


NIELSEN, Ohio State University. 


Degenerate Modes of Vibration and Perturbations in Polyatomic Molecules. W. H. SHAFFER, Ohio State University. 
Absorption Spectra of the Planets. N. Boprovnixorr, Perkins Observatory. 

Progress in the Interpretation of Stellar Spectra. O. Srruve, Yerkes Observatory. 

The Lines of Hydrogen and Oxygen in Stellar Atmospheres. S. CHANDRASEKHAR, Yerkes Observatory. 

Physical Problems Related to the Hydrogen and Helium Lines in Stellar Spectra. M. K. KroGpaut, Yerkes Observatory. 
The Absorption Laws for Gases in the Infra-Red. J. R. N1ELSEN, University of Oklahoma, V. THORNTON AND E. BROCK- 


DALE, Phillips Petroleam Company. 


ABSTRACTS OF CONTRIBUTED PAPERS 


El. A Simple Counting System for Alpha-Ray Spectra. 
S. ROSENBLUM AND W. Y. CHANG, Princeton University.— 
When several tungsten wires are stretched parallel in front 
of a brass plate, it is found that each wire answers alpha- 
particles very well but is not affected by strong beta- or 
gamma-rays. Hence each wire behaves as an alpha-ray 
counter. The voltage of the plate relative to each wire is 
about negative 3000 v, when the distance of each wire from 
the plate is 1.5 mm, the operation taking place in air of 
atmospheric pressure. The conditions of reliable operation 
have been carefully examined and well defined. If a ca- 
pacity of about 500 cm is connected across each counter, 
the power output of even a single power tube is high enough 
to operate a mechanical counter. Removal of the capacity 
needs an amplifier of the multivibrator type to actuate a 
Cenco high voltage impulse counter, but reduces the 
resolving time of each counter considerably. The interfer- 
ence between the different counters can be eliminated by 
erecting aluminium walls between any two adjacent wires. 


When the system of counters is made airtight and each wire 
is properly coupled to an amplifier of the above type, they 
can be used in a deflection chamber under the action of the 
magnetic field to determine the alpha-ray spectra. The 
counter voltage should be well stabilized. The distribution 
form of the Po alpha-particles determined by these counters 
agrees well with that obtained from the photographic line 
(see the other abstract by W. Y. C.). 


E2. A Study of the Alpha-Ray Spectra by the Cyclotron 
Magnet. W. Y. CHANG, Princeton University.—An alpha- 
ray magnet spectrograph has been constructed, consisting 
of the Princeton cyclotron magnet and a Plexiglas de- 
flection chamber, in which the alpha-particles can be bent 
into a semicircle of 75-cm maximum diameter. A balance 
has been devised to measure the magnetic field, which is 
sensitive to about 1 in 10,000. A non-linear deflection coil 
has been used to indicate the variation of the field and can 
also be employed to control directly the magnetic field by 
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allowing a light spot reflected from the mirror of the coil to 
fall on to a twin photo-cell, which then regulates the cur- 
rent of the generator exciting field. Three methods of de- 
tection have been used according to the different strengths 
of the sources. They are, namely, the ordinary photographic 
method, the counting method, and the method of photo- 
graphic tracks. The counters, specially developed for this 
purpose, consist of eight tungsten wires stretched in front 
of a brass plate and are described in another abstract. By 
using the track method a source of one microcurie or 
smaller can be detected. The forms of the energy distri- 
bution of ‘the Po alpha-particles have been determined, 
respectively, by these three methods and agree fairly well 
with one another. The half-width is smaller than 4 mm, i.e., 
about 0.01 Mev. Hence according to the dimensions of the 
chamber, two lines separated by 0.01 Mev can be resolved. 
The intensities a few mm away from the maximum are 
practically zero, being quite different from the results of 
Roy Ringo.* 
* Roy Ringo, Phys. Rev. 58, 942 (1940). 


E3. The Nuclear Excitation of Silver. Marcettus L. 
WIEDENBECK, University of Notre Dame.—The excitation 
of silver' has been studied in the region from the threshold 
for the process to 3.3 Mev. The thick target x-ray* excitation 
curve was obtained by the direct bombardment of an argon 
(90 percent)-ether (10 percent) counter with a silver 
cathode. This curve shows sharp increases in slope when the 
maximum energy of the x-rays reaches a new energy level. 
These levels indicate that the threshold for the process 
Ag+y—Ag* is 1.18+0.03 Mev and higher activation levels 
were also found at 1.59, 1.95, 2.32, 2.76, and 3.13 Mev, re- 
spectively. The decay period of the metastable state was 
obtained from a curve extending over 6 half-lives and was 
found to be 40.4+0.2 sec. The excitation was also obtained 
by direct electron bombardment of silver in the process 
Ag+evr>Ag* 


1 Alvarez, Helmholtz, and Nelson, Phys. Rev. 57, 600 (1940). 
2. R. Feldmeier and G. B. Collins, Phys. Rev. 59, 937 (1941). 


E4. A “Fast” Gamma-Ray Transition. RosaLyn S. 
YaLow M. GoLpHaBER, University of Illinois —We 
have continued the investigation of the “fast"’ gamma-ray 
transition, first detected by O'Neal and Scharff-Goldhaber' 
through the presence of a “‘tellurium-like” satellite among 
the Sb x-rays emitted in the decay of Te". This decay leads 
to two excited states of Sb", of 0.23-Mev and 0.61-Mev 
energy. These energies, determined by absorption of the 
gamma-rays in several elements, are in good agreement 
with the energies found by Hales and Jordan? for the in- 
ternal conversion electrons from Te". By studying the 
coincidences between the x-ray satellite and the internal 
conversion electrons we could show that the transition from 
the 0.23-Mev level to the ground state is the “fast” transi- 
tion, taking place before the K-electron captured by Te™ 
is replaced. Its internal conversion coefficient is surprisingly 
large, >0.5, though one K-electron is missing! The mean 
lifetime of this state is <2X10- sec. and the corre- 


sponding half-width >30 ev. The transition from the 0.61- 
Mev level to the ground state is “‘slow,”’ and its half-width 
<1 ev. 


1R, D. O'Neal and G. Scharff-Goldhaber, Phys. 
2 E. B. Hales and E. B. Jordan, Phys. Rev. 202 ( 


. 62, 83 (1942). 
1943). 


ES. Attempt to Detect Nuclear Resonance Absorption 
of Gamma-Rays. GERTRUDE SCHARFF-GOLDHABER, M. 
GOLDHABER, AND RosaLyn S. YALow, University of Illinois. 
—Since past attempts to discover nuclear resonance ab- 
sorption with ThC” gamma-rays were unsuccessful, it 
seemed of interest to search for resonance absorption with 
a gamma-ray which appeared to be particularly suited for 
such an experiment: the 0.23-Mev gamma-ray' emitted 
from an anomalously wide level of Sb"*', resulting from the 
decay of Te". However, we were unable to detect a nuclear 
absorption of this gamma-ray in Sb; if present it is <3 
percent of the atomic absorption and at least a hundred 
times smaller than expected. A search was therefore made 
for a low lying excited state of Sb", but no gamma-rays 
above 10 kev were found. If the 0.23-Mev transition leads 
to the ground state of Sb™, a new factor must be responsible 
for a shift of the gamma-ray line from the resonance posi- 
tion. As this “‘fast’’ gamma-ray is emitted with one K-elec- 
tron missing, while it would havg to be absorbed in a 
neutral Sb" atom, it may be that the interaction of the 
nucleus with the “hole in the K-shell” leads to a consider- 
able shift of the nuclear level from its ‘‘normal’’ position. 


1 See preceding abstract by R. S. Yalow and M. Goldhaber. 


E6. Radiation Widths of Highly Excited Nuclei. M. 
GoLDHABER, University of Illinois —Recent evidence for a 
wide low energy level in a medium weight nucleus (Sb™) 
appears to be difficult to reconcile—within the framework 
of the compound nucleus theory—with the sharp states 
found for resonance capture of slow neutrons. In that 
theory the existence of sharp resonance states, ~9 Mev 
above the ground state, is ascribed to the “slowness” of 
gamma-ray transitions, as observed for some heavy 
radioactive nuclei. For medium weight nuclei it may be 
necessary to treat the sharp “initial capture state” formally 
as a particular state of the combination Z” +n rather than 
as one of the many states of the highly excited compound 
nucleus Z**1, If such a view were taken many experimental 
results would have to be reinterpreted. In particular, the 
empirical rule that sharp resonance states are predomi- 
nantly found for nuclei of odd mass number, may take on 
new significance. For an odd M nucleus the magnetic 
interaction with a neutron is comparable in magnitude to 
the energies found for neutron resonance. For even M nuclei 
this interaction is usually absent. This correlation suggests 
that the influence of the magnetic forces on the position and 
width of the “‘initial capture state’’ may be important. 


E7. Effect of Nuclear Electric Quadrupole Moment on 
the Energy Levels of a Diatomic Molecule in a Magnetic 
Field. BERNARD T. FELD AND E. Lams, Pupin 
Physics Laboratories, Columbia University —In molecular 
beam experiments' performed on diatomic molecules of zero 
electronic spin, effects were observed which seemed to indi- 
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cate an electrical interaction between the nucleus and the 
rest of the molecule. It was suggested that these effects are 
caused by an interaction between a nuclear electrical 
quadrupole moment and the electric field of the molecule. 
We have calculated the energy level system and expected 
molecular beam resonance spectrum for heteronuclear 
molecules in which one nucleus has quadrupole moment and 


_ the other zero spin. The calculations were carried out, by ~ 


the use of perturbation theory, for nuclear spins 1 through 
9/2 in weak and strong magnetic fields. It is shown that the 
observed effects could be owing to a nuclear electric 
» quadrupole moment. Analysis of unpublished data on reso- 
nances in LiF give for Li’ the value qQ=1X10~* erg (Q is 
the nuclear electric quadrupole moment; g, a quantity, 
depending on the distribution of charge in the molecule, 
which must be independently calculated). This may be 
compared with 0.3 10~* erg for D in HD and D; 


1 Kusch, —, and Rabi, Phys. Rev. = _ (1939). 
2 Kellogg, Rabi, Ramsey, and Z Rev. 57, 677 (1940). 


E8. Neutron-Proton Scattering and the Meson Theory 
of Nuclear Forces. J. M. Jaucn, Palmer Physical Labora- 


tory, Princeton University—The ratio R=o(x)/o(x/2) of 


the scattering cross section for 14-Mev neutrons on protons 
has been calculated assuming for the interaction the ex- 
pression which results from the charge symmetrical meson 
theory with two kinds of mesons of masses « and y(y>«). 
The previously determined values for the coupling constant 
and the mass ratio u/«' are used. The value of the ratio R 
is essentially determined by the sign of the interaction in P 
states. The sign is such as to give repulsive forces for the 
charge symmetrical theory. Numerical calculations have 
been carried out to an accuracy for R of ~1% for all states 
of angular momentum =2. The effect of the higher states 
can be estimated with an approximation method and is 
found to be negligible. The result is R= 1.11. This value is 
at variance with the recently determined? ratio of 0.52+0.03 
for the same energy. This discrepancy supported by other 
similar calculations* is a strong argument against “"~ 
a symmetrical meson theory. 
. M. Jauch and Ning Hu, Phys. Rev. 65, 289 (1944). 
Amaldi et al., Naturwiss. 30, 582 (1942). 
°L Hulthén, Arkiv f. Mat. Astronom. o. Fysik 29, No, 33 (1943); 


B. Ferreti, Nuovo Cimento [XXI] 1, 25 (1943); J. Schwinger and W. 
Rarita, Phys. Rev. 59, 556 (1941). 


E9. The Influence of the Recoil of Heavy Particles on 
the Nuclear Potential Energy. J. Lerre Lopes, Princeton 
University—An attempt to improve the value of the 
quadrupole moment of the deuteron in Schwinger’s mixed 
meson theory of nuclear forces! was made by taking into 
account the recoil of the nucleons during the emission and 
absorption of mesons. The interaction energy operator of 
the system proton-neutron was calculated by perturbation 
theory and represented by an integral operator in mo- 
mentum space. In the pseudoscalar and vector symmetrical 
theories the tensor force has a singularity of the type 1/r* 
in configuration space; in Schwinger’s mixture it has a 
singularity of the form 1/r*. This singularity still remains if 


relativistic terms are considered, even if recoil is neglected. 
An evaluation of the quadrupole moment of the deuteron 
is therefore only possible with a cut-off procedure and this 
would make unimportant the relativistic and recoil contri- 
butions as well as the mixed theory itself. 


1J. M. Jauch and Ning Hu, Phys. Rev. 65, 289 (1944). 


E10. The Adaptation of the Cauchois Spectrograph to 
Artificial Radioactive Sources. J. E. Epwarps, M. L. Poot, 
AND F. C. BLaKE, Ohio State University.—To identify x-rays 
emitted by an artificial radioactive source it is usually 
necessary to distinguish x-ray wave-lengths characteristic 
of two neighboring elements. The feeble x-rays from such 
sources can be photographed with a curved crystal spectro- 
graph. Quartz and mica crystals have been used in the 
transmission (Cauchois) type curved crystal spectrograph 
to photograph characteristic Ka and K@ x-ray lines from 
radioactive sources. Quartz crystals adapted to this type 
of spectrograph by DuMond and Watson were found to be 
most efficient for the longer wave-length range of the 
instrument from about 0.7A to 1.9A. Mica crystals proved 
to be more efficient for shorter wave-lengths from about 
0.36A to 1A. Spectrographs with curvatures of 8” and 15” 
were used with both mica and quartz crystals. Broad 
sources resulting from chemical separations may be used in 
this type spectrograph. Line sharpness is decreased by the 
broad source but the Ka lines of two neighboring elements 
are easily identified. The characteristic silver x-rays associ- 
ated with the decay of 6.7-hour Cd'®.!0 were photographed 
with a mica spectrograph. The nickel x-rays associated 
with the decay of Cu™ have been identified by use of a 


quartz spectrograph. 


Ell. Radioactive Zirconium and Columbium. M. L. 
Poot AND J. E. Epwarps, Ohio State University.—Ex- 
tended observations on the decay of Zr® and Cb® show a 
genetical relationship between these two isotopes. Zr™ is 
characterized by a half-life of 67.8 days, a gamma-ray of 
0.85 Mev, a strong B-ray of 0.29 Mev, and a weak one of 0.8 
Mev. The ratio of the intensities is approximately 8. Cb” 
is characterized by a half-life of 38.7 days, a gamma-ray of 
0.78 Mev, and soft electrons of 0.140 Mev. The gamma- 
radiation from Zr® rises in intensity with time after a 
chemical separation of the Cb®; the 8-radiation, however, 
decreases. Cb x-rays have been observed by absorption and 
also by photography in a curved crystal spectrograph. The 
ratio of the number of x-rays to that of y-rays is approxi- 
mately 4. A new reaction Y* (a, n) Cb” gives Cb” very 
strongly. The half-life is 10.1 days. A gamma-ray of 1.1 Mev 
is emitted. X-rays have been established by absorption 
measurements and are about as intense as the y-rays. By 
the spectrograph these x-rays are shown to be those of Zr. 


E12. Progress in the Determination of the Number of 
Artificially Produced Radioactive Atoms. J. D. KuRBATOV 
AND M. H. Kurpatov, Ohio State University.—A determi- 
nation of the number of atoms obtained by some means of 
activation requires the development of a method for — 
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quantities less than 10~* gram. A possible solution of this 
problem has been tested using the adsorption phenomena 
of divalent ions on different hydrous oxides. The method 
consists in obtaining adsorption isotherms for known minute 
quantities of divalent ions under conditions such that the 
exponent in the activated adsorption equation is consider- 
ably below unity; then an unknown initial quantity can be 
found from the adsorption ratio. The adsorption constants 
have been evaluated for radium since a technique is known 
for its accurate determination in quantities of less than 10-* 
gram. The investigations have been extended to obtain the 
adsorption isotherms of minute quantities of barium and 
strontium using the radioactive tracers Ba™ and Sr®*, 
However, it has been found experimentally that with 
conventional quantities of adsorbent, that is more than 10 
mg, the exponent in the adsorption equation became unity 
for quantities of Ra, Ba, and Sr below 10~* gram and the 
fraction of ions adsorbed was constant. Thus, to evaluate 
minute quantities of Ra, Ba, and Sr it appears necessary 
to use the adsorbents, ferric hydrous oxide or manganic 
oxide (IV), in quantities less than 10~ g. 


Fl. The Hardy Coefficient for Instrument Oil. S. 
BLOOMENTHAL, Automatic Electric Company.—The method 
of inclination, ascribed to W. B. Hardy, appears to be’ the 
most straightforward of many ways for comparing lubri- 
cants. The coefficient of static friction u for glass on glass 
was measured by the inclination method on an apparatus 
which functions automatically. A number of instrument 
oils, suggested as possible substitutes for genuine porpoise 
jaw oil, which is unobtainable at the present time, were 
tried. Addition of only 1 percent glyceryl monoricinoleate 
to a pure lubricating fluid caused a decrease in u from 0.26 
to 0.13 because of boundary layer formation. The disclosure 
on last August 15 by G. E. Barker of the composition of a 
non-spreading synthetic instrument oil, made available to 
the public by the Elgin National Watch Company, is note- 
worthy in this connection. This oil contains a stable aralkyl 
ether with an aliphatic chain of at least five carbon atoms 
and a benzene nucleus connected to the ether oxygen atom. 
u for M56A Elgin Watch Oil was found to be 0.17 when 100 
grams total load is used. This compares with 4 =0.26 for a 
light oil derived entirely from petroleum, 


F2. Equations of State and Phase Changes in Films on 
Liquids and Solids. Wmt1am D. HARKINS AND GEORGE 
Jura, University of Chicago and Universal Oil Products 
Company.—Oil films on water seem not to be stable unless 
they are monomolecular but adsorbed films on solids may 
be monomolecular but more generally are multimolecular 
in the sense that they may attain thicknesses equivalent to 
from five to ten or even more molecular layers. Monolayers 
on water exhibit three condensed and three non-condensed 
phases, or six in all. The relation between film pressure (x) 
and molecular area (¢) is linear (r=b—ae, equivalent to 
log ~/po= B—A//v*) for any condensed phase, but for liquid 
intermediate or expanded films or for gaseous films the 
equations are more complicated. At sufficiently low temper- 
atures films adsorbed on solids are often condensed at the 


higher vapor pressures, but gaseous at very low pressures. 
At higher corresponding temperatures either liquid ex- 
panded or liquid intermediate phases may occur. Phase 
changes which are first order in three dimensions usually 
become second order in two dimensions, though in the 
latter the condensation of a vapor may be either first or 
second order. 


F3. The Effect of the Subphase upon the Properties of 
Physically Adsorbed Films. GEorGE JURA AND WILLIAM D. 
Hargis, University of Chicago and Universal Oil Products 
Company.—Almost all theoretical approaches to the prob- 
lem of the physical adsorption of gases on the surfaces of 
solids have assumed a negligible interaction between the 
solid and the film. Recent careful experimental work shows 
that this is incorrect. A few effects which show this are 
listed. (1) When the effective cross-sectional area of nitro- 
gen on 100 solids is calculated it is found that three maxima 
are obtained at 14.0, 15.2, and 16.2A?*. It is found that in 
some cases at least there is a relation between the area 
exhibited and the composition or the crystal structure of 
the solid. (2) The film of propyl alcohol on barium sulphate 
is monomolecular, but polymolecular on anatase (titanium 
dioxide) and quartz. The thickness of the film is related to 
the lattice constants of the solids. (3) The transition from 
a gaseous to an intermediate film is first order at 1000A? per 
molecule for one solid while for another of the same crystal 
structure and particle size, the transition is found to be 
second order at 50A* per molecule. 


F4. Molecular Weights of Proteins as Determined by 
Light Scattering. W. HELLER AND H. B. KLEVENS, 
Chemistry Department, University of Chicago——Molecular 
weights can be determined by light scattering by using 
either the theoretical principle of Mie\—light scattering 
is considered as an effect of material units—or the principle 
of Einstein*—light scattering is considered as an effect of 
fluctuations in concentration. Debye’s recent equation 
based on the latter principle could be verified on suspen- 
sions of polymers and on protein solutions by using a 
spectrophotometric-interferometric method (18-cm layer). 
The symmetry of radiation was tested by means of 
Tyndall spectra.‘ Aqueous solutions of egg albumin gave, 
for example, the following data (pw: 4.65-4.68): 


47,000+ 1500 
64,000 + 2000 
81,000+ 3000. 


Immediately after preparation (4°C) 
24 hours later (4°C) 
72 hours later (4°C) 


The ultracentrifugal and osmotic pressure data (literature) 
are 40,500-44,000 and 42,000-45,100, respectively. Similar 
agreement between light scattering and ultracentrifugal 
data was found for serum-albumin. It is apparent that the 
method is extremely sensitive towards aggregation (de- 
naturation) of proteins. Symmetrically radiating polymer 
suspensions gave an agreement within 10 percent with 
ultramicroscopic data of H. Oppenheimer. 

1 Mie, Ann. d. Physik 25, 377 (1908). 

2A. Einstein, Ann. d. Physik 33, 1275 (1910). 

Debye, J. A 

‘W.H and 


. Phys. 15, 338 sage 
. Vassy, Phys. Rev. 63, 65 (1943). 
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FS. Diode Modulation. G. H. Fetr ann A. D. Barey, 
Department of Electrical Engineering, University of Illinois. 
—It is shown analytically that a perfect diode in series with 
a pure resistance load is capable of producing frequency 
components required in amplitude modulation either when 
the effect.of the diode is considered as a series of pulses 
analysed by means of Fourier series, or when the diode is 
considered as a synchronous commutator. The end result 
is identical in either case. The experimental verification of 
the mathematical result is given for a pure resistance load, 
and it is shown that if account is taken of the resistance of 
the diode the agreement is good. It is also shown that the 
maximum modulated output takes place when the load 
resistance is about equal to the diode resistance. The 
problem of the parallel resonant circuit load is discussed, 
and it is shown that the diode under these conditions be- 
haves like a variable effective resistance, the magnitude of 
which is controlled by the magnitude and sign of the 
modulating signal in series with the carrier. When the Q of 
the load is high, complete modulation is obtained with about 
half of the signal voltage required with a resistance load. 


F6. The Multiple Production of Mesotrons in Paraffin 
at High Altitudes. W. G. Srroup AND MARCEL’ SCHEIN, 
University of Chicago.—By modifications of a technique 
used in previous balloon flights,| more information re- 
garding the production of mesotrons in the stratosphere has 
been obtained. The number of mesotron showers, involving 
three or more particles, under three and six cm of paraffin, 
have been measured and compared with the curve for the 
mesotron intensity as a function of atmospheric pressure. 
The results show that at altitudes between 13 and 16.5 km 
10 percent of the recorded penetrating particles are ac- 
companied by mesotron showers of at least three particles 
below the paraffin. The angular divergence of these showers 
is between 9 and 20 degrees. Showers of penetrating par- 
ticles were also detected below 6 cm of paraffin plus an 
additional 8 cm of lead. It was found that in the above- 
mentioned altitude range more than 20 percent of the 
penetrating particles are accompanied by mesotron showers 
below the paraffin plus lead. Many of these showers were in 
coincidence with those below the paraffin only which 
demonstrates that a large fraction of the particles generated 
in the paraffin are able to penetrate 8 cm of lead. 


1M. Schein, W. P. Jesse, and E. O. Wollan, Phys. Rev. 59, 615 
fees M. Schein, M. Iona, Jr., and J. Tabin, Phys. Rev. 64, 253 


F7. Investigation of Bursts Observed in Two Thin- 
Walled Ionization Chambers. Konrap L. KINGSHILL (in- 
troduced by Marcel Schein) anp LLoyp G. Lewis, Uni- 
versity of Chicago.—Coincidences in time between bursts of 
ionization occurring in each of two thin-walled ionization 
chambers were measured at Chicago. The apparatus used 
here was the same as that used at Echo Lake, Colorado 
(elevation 3100 meters). The cumulative size-frequency 
distribution curve for all bursts recorded in each chamber is 
essentially in agreement with the data obtained by previous 
investigators. The altitude comparison was made with the 


data obtained at Echo Lake. The ratio of burst rates found 
at Echo Lake and at Chicago for bursts caused by 50 or 
more cosmic-ray particles was approximately 10:1. This is 
in general agreement with Geiger-Mueller counter coinci- 
dence measurements. For bursts caused by more than 150 
particles the corresponding ratio is 200:1. Since burst 
coincidence rates measured in two ionization chambers 
showed a similar altitude effect when the chambers were 
separated by a distance of 1 meter, it can be assumed that 
both the large single bursts and the coincidence bursts are 
caused by the same type of atmospheric shower. The rate of 
occurrence of burst coincidences was found to be very small 
in comparison to the rate for bursts occuring in a single 
chamber. This demonstrates that most of the atmospheric 
showers must have a very small lateral width of their region 
of high particle density. Provided that the cascade theory 
of showers is valid for electron energies higher than 10" ev, 
it can be concluded that sea-level showers exhibiting such 
narrow high density regions (of 1-meter extension) cannot 
originate from primary electrons. 


F8. Investigation of Extensive Atmospheric Showers in 
the Stratosphere. P. AuGeR, A. RoGOZINSKI, AND MARCEL 
SCHEIN, University of Chicago ——To measure the frequency 
of extensive atmospheric showers at high altitudes a 
counter apparatus of large linear dimensions (10 meter) 
was constructed and sent to the stratosphere. The essential 
part of the apparatus consisted of six large Geiger-Mueller 
counters (area 125 cm?) arranged in a horizontal plane. The 
maximum separation between counters was 7 meters. 
Another set of 4 small counters (area 15 cm?*), separated by 


’ $0-cm distance, was placed in the middle between the large 


counters. Six different combinations of threefold coinci- 
dences were registered simultaneously. The coincidence 
circuits had a resolving time of 2 10~* second and there- 
fore chance coincidences could be neglected. The results 
show that from sea level to an altitude of 15 km the shower 
rate increases by a factor of 4 for showers of 2-m average 
spread and by a factor of 30 for showers of 50-cm average 
spread. It was also found that in the high altitude no 
coincidences occurred between the set of small counters and 
any of the large ones which was frequently the case at sea 
level. This demonstrates that at 15-km altitude the air 
showers consist of a very narrow and dense bundle of 
particles. According to the cascade theory these showers 
must be in the initial stages of their development. 


G1. A Method for Obtaining Depolarization Factors of 
Raman Lines. Forrest F. CLEVELAND, Jilinois Institute of 
Technology.—In many cases it is possible to determine the 
structure of a molecule by interpreting the Raman and 
infra-red spectra with the aid of the group theory selection 
rules for two or more assumed structures. To do this, one 
must know how many lines are depolarized (depolarization 
factor = p= 6/7) and how many lines are polarized (p <6/7). 
This requires very precise values of those depolarization 
factors whose values are near 6/7 in order to decide whether 
these lines are polarized or depolarized. Consequently, the 
method previously used has been revised with a view 
toward more precise results. A 90° split-field, Polaroid disk 
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is placed in contact with the window of the Raman tube, 
with the dividing line horizontal, and is sharply focused on 
the slit of the spectrograph. A single exposure is made. 
Calibration marks corresponding to known relative in- 
tensities are placed upon each plate by use of a ground-glass 
screen and a small incandescent lamp. One-minute ex- 

ures are used for the calibration marks and correction 
for failure of the reciprocity law is made by use of informa- 
tion obtained in preliminary experiments. Intensities are 
obtained by use of a Gaertner microdensitometer. Results 
thus obtained for CCI, will be compared with data obtained 
by Reitz. 


G2. Raman Spectrum of Dimethyldiacetylene. ARNOLD 
G. MEISTER AND Forrest F,. CLEVELAND, IJilinois In- 
stitute of Technology —Raman frequencies, relative inten- 
sities, and depolarization factors were obtained for 
dimethyldiacetylene. The relative intensities and the 
depolarization factors were obtained with the use of a 
Gaertner microdensitometer. A tentative assignment was 
made of all the observed Raman frequencies assuming that 


the molecule has the symmetry Dy. The chain frequencies 
were established by comparison with Wu and Shen's! 
analysis of the vibrational spectrum of diacetylene. 
Crawford’s® analysis of the dimethylacetylene Raman spec- 
trum was used in determining the methyl group frequencies. 
Three cases of Fermi resonance were observed, but there 
were also two cases where only a single line was observed in 
spite of the fact that the frequency and symmetry condi- 
tions for resonance were satisfied. The carbon-carbon triple- 
bond frequency which appears at 2183 cm for diacetylene 
is shifted considerably, appearing at 2264 cm for 


dimethyldiacetylene. 
1T. ¥. Wu, Vibrational Spectra and Structure of Polyatomic Molecules 


(National University of Peking, Kunming, China, 1939), p. 259. 
2B. L. Crawford, Jr., J. Chem. Phys. 7, 555 (1939). 


G3. A Difficulty in Analysis of Some ClO, Bands. 
J. B. Coon, Department of Physics, Indiana University.— 
The rotational structure of the violet absorption spectrum 
of ClO; cannot be explained in terms of energy levels of a 
rigid, approximately symmetrical top. If these levels are 
used, an additional term linear in K is needed. The spec- 
trum can, however, be explained using the levels* 


F= +1) 


Replacing k by K gives levels of a rigid, approximately 
symmetrical top. Here the symbol A refers to the quasi- 
symmetry axis of the approximately prolate top model. 
L, K, and k refer, respectively, to the total angular mo- 
mentum excluding spin, component of this along the axis, 
and the component of the angular momentum of molecular 
rotation along the axis. Though & is not integral, it is as- 
sumed that any change in k is integral. y= F’— F” (AL=AK 
=Ak=0) was fitted to the observed Qg branches. Within a 
given band it was found that 2A4A is constant. Using 
constants of the spectrum, k was calculated and k—K was 


found to be a positive constant within a given band, 1.68 
for (000-000) band. 


* H. Sponer and E. Teller, Rev. Mod. Phys. 13, 75 (1941). 


G4. New Evidence for the Magnetic Current. Fe.ix 
EHRENHAFT, New York, New York.—In a homogeneous 
electric field through which flows a constant electric cur- 
rent, a single magnetic pole circulates. If the single mag- 
netic pole also bears an electric charge, this movement is a 
spiral (electrodynamics). Observation discloses that in a 
constant vertical homogeneous magnetic field, single test 
bodies move in a spiral path in gases as well as in liquids 
with velocities that are measurable in the dark field of the 
microscope. These velocities are of the order of 10°10" 
smaller than the velocity of light, showing that they bear 
electric as well as magnetic charges (magnetodynamics). If 
they bear a magnetic charge alone, their movement is 
polar.* New experiments will be demonstrated in gases 
where the permanent magnet is used. Covering the vertical 
magnetodes of the electromagnet with glass caps, the 
opposing faces of which are fitted with platinum surfaces, 
one obtains in the space between the caps vertical homo- 
geneous fields where the electric and magnetic force lines 
are parallel to one another. When the two fields are in 
action simultaneously a circulation of the liquid can be 
observed with the microscope and a twisting of the stream 
of electrically evolved gas bubbles rising from the lower 
platinum surface occurs. The direction of twist as well as 
the direction of circulation follows the rule of Ampere. New 
facts about magnetolytic processes will be given. 


* F. Ehrenhaft, Nature 154, 426 (1944). 


GS. Magnetic Saturation of Microscopic Particles. 
BROTHER GABRIEL KANE, Manhattan College.-—Microscopic 
particles of nickel have been observed to move in the direc- 
tion of the lines of force in a vertical homogeneous magnetic 
field when illuminated by a strong beam of light.' The 
velocity of the particles is a function of the field strength. 
The direction of motion reverses with a reversal of the 
magnetic field. For low field intensities, the velocity is 
directly proportional to the field strength. For higher 
intensities, the velocity either approaches a limiting value 
or reaches a maximum with increasing field and then 
decreases. Increase of field strength may even reverse the 
sense of direction of the motion. All of these motions can be 
accounted for if it is assumed that the force acting on a 
particle is gH, where g is the magnetic charge and H the 
effective field strength. The value of H is determined from 
the equation H=H’—NI used for the macroscopic state. 
H’ is the applied field, N the demagnetizing factor, and J 
the intensity of magnetization. N is a function of the field 
strength.? 

Ann. de physique 13, 151 (1940); Phys. Rev. 57, 562 


uses) et: I 230, 381 (1940) 
2 PL Phys. Rev. 31, 165 (1910). 


G6. The Rotation of Electrolytes under the Influence of 
a Magnetic Field. Cuartes B. Reynoips, Federal Com- 
munications Commission, New York, New York.—The rota- 
tion of solutions of CuSO,, H,SO,, and FeCl; under the 
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influence of a magnetic field has been reported.! The author 
has repeated this experiment using many additional 
electrolytes of ferro-, para-, and diamagnetic materials, 
organic and inorganic in nature. The experiments performed 
_ to determine whether or not the motions could be explained 
by the Lorentz force, F=q(vXH/c), or in terms of some 
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other well-known theory, will be described. It will be shown 
that the rotations are not explained by the Lorentz force, 
concentration gradient, or pure electrochemical action. A 
satisfactory explanation can be shown if the existence of a 
magnetic current is assumed. 

1F, Ehrenhaft, Bull. Am, Phys. Soc. 19, 3, 10 (June 23, 1944). 


Proceedings of the American Physical Society 


MEETING AT PASADENA, CALIFORNIA, DECEMBER 16, 1944 


HE 264th meeting of the American Physical 
Society was held at the California Institute 

of Technology, Pasadena, California, on Satur- 
day, December 16, 1944. The attendance varied 
from 90 to 125, the largest at any Pacific Coast 
meeting since the start of the war. Thus the trend 
of increasing numbers, started a year ago, con- 
tinues in evidence. The audience was, however, 
almost exclusively from Southern California. 
Nearly 50 members and guests had luncheon 
together in Dabney Hall, the shift from the 
Athenaeum having been forced by war conditions. 
The morning session opened at 9:30 A.M. and 
adjourned at 12:30 p.m. It was devoted to three 


invited papers, viz: 


The Conceptual Content of the Quantum 
Theory, S. Epstein, California Institute 
of Technology. 

The Logical Foundations of Quantum Me- 
chanics, HANS REICHENBACH, University of 
California at Los Angeles. 

The 1944 Values of the Physical Constants 
of the Quantum Theory, RaymMonp T. Brrcg, 
University of California. 

The afternoon session opened at 2 P.M. and 
adjourned at 4:30 p.m. Abstracts follow of the 
twelve contributed papers presented at that 
session. Number 10 was read by title. 


R. T. BirGE 
Local Secretary for the Pacific Coast 


ABSTRACTS OF CONTRIBUTED PAPERS 


1. A Possible Cause of the Variability with Time of the 
Incoming Cosmic Rays at High Latitudes. Ropert A. 
MILLIKAN, H. Victor NEHER, AND WILLIAM H. PIcKER- 
ING, California Institute of Technology.—(1) We present 
new evidence that a new band of rays which we interpret 
as helium annihilation rays comes in vertically at about 
the predicted latitude. (2) We present a discussion of the 
possible composite character of the so-called silicon an- 
nihilation band and of the so-called oxygen annihilation 
band. (3) We bring forward an explanation of a possible 
cause of our observed variability in the incoming cosmic 
rays at high latitudes. (4) We make a new and more ac- 
curate determination of the value of the field sensitive 
and the non-field sensitive components of the incoming 


cosmic rays. 


2. Some Recent Investigations on Hydrodynamic Sta- 
bility. C. C. Lin, Guggenheim Laboratory, California Insti- 
tute of Technology (Introduced by R. A. Millikan).—Recent 
investigations of the stability of two-dimensional parallel 


flows bring out that the physical mechanism underlying 
hydrodynamic stability in an incompressible fluid can be 
understood in terms of vorticity transfer.’ In particular, 
the condition of vanishing curvature of the velocity profile 
in Rayleigh’s criterion is identified with an extremum of 
vorticity distribution, and the physical picture describing 
the associated instability has been worked out in some 
detail. Conditions similar to Rayleigh's criterion have now 
been obtained for curved and axially symmetrical flows, 
and these can also be interpreted in a similar manner. 
Investigations of the stability in a compressible fluid are 
being carried out with the collaboration of L. Lees. 


' Proc. Nat. Acad. Sci. (October, 1944). 


3. Pressure Flow of a Turbulent Fluid between Two 
Parallel Infinite Planes. P. Y. Cuou, Guggenheim Labora- 
tory, California Institute of Technology (Introduced by R. 
A. Millikan).—This investigation is based upon a general 
theory developed by the author: ‘On velocity correlations 
and the solutions of the equations of turbulent fluctuation.” 
Two solutions of the problem are given. The first is the 
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consequence of the equations of mean motion and of double 
correlation, while, in the second solution, all the equations 
of mean motion and of double and triple correlations are 
utilized. It is shown that the first solution holds good only 
in the central portion of the channel, while the second 
covers the entire channel except the immediate vicinity of 
the channel wall. The behavior of the mean squares of 
turbulent fluctuations is discussed qualitatively. The pres- 
ent calculation also indicates that the velocity defect 
(U.— U)/U., which is independent of the Reynolds number 
of the mean flow, according to von Karman,? may also be 
independent of the magnitudes of the turbulent velocity 
fluctuation, after the flow has reached its steady turbulent 
state. 


vom of pl Com. 
1 +P. 


4. Frequency Elimination in Spirotron Systems for Ac- 
celerating Ions and Electrons. L. E. Dopp, University of 
California at Los Angeles—The spirotron design changes 
the speed of an electron or ion by components, with rela- 
tively low frequency of potential-reversal. It combines the 
cyclotron and magnetic-electron-microscope principles. 
Further modification eliminates frequency by using con- 
stant potential. This is getting back somewhat to former 
designs before the advent of the cyclotron. However, a 
difference is that here smaller diameter magnetic fields are 
permitted. Multiple electrode pairs, oppositely arranged on 
the surface of a cone, are replaced by a single pair of long 
straight electrodes, of the same general character as before, 
parallel to the magnetic field. The particle moves on the 
surfaces of a sequence of right-cylindrical helixes. These, 
with successive changes of radii, are joined together by a 
single straight line—the axis of the one electrode pair— 
parallel to the axes of all cylinders, but located in the 
surface of each, this common line and the cylinder axes 
being confined to the same plane. The particle travels a 
whole revolution between accelerations, instead of a half- 
revolution as in the cyclotron. Its path during this time is 
on one of the cylinders. Constant potential avoids fre- 
quency difficulties from relativistic change of mass at higher 
velocities. 


5. Minimum Conditions Necessary to Achieve the Ve- 
locity of Escape. H. S. Seirert, California Institute of 
Technology.—The maximum velocity which may be 
achieved by a rocket directed vertically upward depends 
upon five basic parameters: (1) The velocity (relative to 
the rocket) with which mass is ejected; (2) the fraction of 
the total mass which is finally ejected, or ‘‘loading factor;” 
(3) the initial acceleration; (4) the dependence of accelera- 
tion upon time (flight law); (5) the dependence of frictional 
resistance upon time and velocity. The flight laws of con- 
stant acceleration and constant thrust have been con- 
sidered. Although constant thrust produces higher ter- 
minal velocities than constant acceleration, the terminal 
accelerations produced in an efficient rocket by constant 
thrust are greater than a human passenger could tolerate. 
Assuming an initial acceleration of g, an energy loss due 


to friction equal to the kinetic energy of escape, constant 
thrust, and a step rocket consisting of three parts of mass 
ratios 16:4:1, each with a loading factor of 0.9, the attain- 
ment of escape velocity would require an ejection velocity 
over 3000 meters per second. 


6. The Relativistic Increase in Ionization of Cosmic-Ray 
Electrons. WAyNE E. Hazen, University of California.— 
The preliminary experiments of Corson and Brode' have 
been continued with improved techniques and a stronger 
magnetic field (4500 gauss). It is believed that uncertainties 
in condensation efficiency and curvature no longer con- 
tribute significantly to the errors. The overlapping of drop 
images and its variation with track diffuseness probably 
limit the accuracy of both absolute and comparative ioniza- 
tions to 10 percent. The minimum ionization, which was 
obtained from 35 tracks of P® beta-rays with Hp>4500, 
is 40 ion pairs per cm in air at N.T.P. The average ioniza- 
tion of 15 cosmic-ray electrons with Hp of 1-8 10° is 53. 
These figures for ionization are for the ionization produced 
as a result of all collisions in which the primary electron 
loses less than ~ 600 ev. Theory predicts an increase of 40 
percent, which agrees with the experimental results within 
the estimated uncertainty. 


1D. R. Corson and R. B. Brode, Phys. Rev. 53, 773 (1938). 


7. The Diffusion Problem and the Theory of Connected 
Events. DonaLp A. DarLtnG, California Institute of Tech- 
nology (Introduced by R. A. Millikan).—An extension of 
the one-dimensional random walk to the case where the 
particles possess inertia leads to a statistical problem in 
the theory of non-independent, or Markoff, processes. The 
modification of the “homogeneous random process” in this 
case is deduced, and the meaning to be attached to random 
events extended to the continuum is studied. The analogue 
to the ordinary diffusion equation becomes an integro- 
differential equation 

aptt, d*plt, x) 
at 


where p(t, x) is the probability that the particle (possessing 
inertia) has a displacement x at a time #, and ®(t, y) is a 
solution to the ordinary diffusion equation. Solutions of 
the equation are deduced by considering the moments of 
the probability distribution. Some of the equations used in 
turbulence and Brownian movement problems are special 
cases of the integro-differential equation, with a suitable 
interpretation of the dependence function ®. 


8. Magnesium Ion Source for High Intensity Mass 
Spectrograph. T. H. P1, California Institute of Technology 
(Introduced by R. A. Millikan).—Magnesium ions were 
produced by bombarding magnesium vapor by electrons 
emitted from oxide-coated cathodes. The ions were then 
focused into a thin sheet by means of an electrostatic lens. 
This ion sheet was sent into the magnetic lens, which 
focused the sheet into a small spot where the different 
isotopes were collected. The ion current obtained at the 
collector was 120 microamperes with good resolution. 


10wn 
orce, 
m. A 
ofa 
tum 
itute 
Me- 
of 
ants 
2GE, 
and 
the 
that 
St 
lying 
n be 
ular, 
rofile 
m of 
ibing | 
some 
now 
lows, 
nner. 
1 are 
Two 
bora- 
y R. 
neral 
tions 
the 


66 | AMERICAN PHYSICAL SOCIETY 


9. The Movement of Bed Particles in a Turbulent per unit atomic weight; H/Z = 1837.17 +0.24, ratio of mass 
Stream. Hans ALBERT EINSTEIN, California Institute of of hydrogen atom to that of electron. 


Technology (Introduced by R. A. Millikan).—The current 
theory of suspensions in turbulent streams uses an equi- 
librium condition between the continuous settling of the 
particles in the surrounding water and the vertical turbulent 
exchange of water masses to establish the law for the ver- 
tical distribution of sediment concentration. At the bed 
itself this argument breaks down because here the exchange 
of liquid approaches zero. However, an exchange of par- 
ticles between the bed and the moving layers exists, and 
like all such acts of exchange, is governed by the general 
rules of probability. The probability of a particle moving 
out of its position in the bed can be related directly to the 
probability of the adjoining liquid having a certain velocity. 
If one introduces the measured probabilities for instan- 
taneous velocities near the bed, the law for the movement 
of particles out of the bed can be predicted up to some 
general constants. The law obtained describes very closely 
the sediment movement in laboratory flumes and natural 
rivers. 


10. The Rydberg Constants and the Value of e/m. CHao- 
WanG HsveEu, California Institute of Technology.—The 
wave-lengths of the H-alpha and D-alpha lines were meas- 
ured with a Fabry-Perot interferometer directly against 
the primary standard of wave-length, 6438.4696IA. The 
results of the measurements are: 6562.85337+0.00022I1A 
for H-alpha, and 6561.06811+0.00005IA for D-alpha. 
Following William’s assignment,' they are the respective 
wave-lengths of the strong component-group (145). From 
the above values, the Rydberg constant for hydrogen 
Ru, for deuterium Rp, for helium Rye, and for infinite 
mass R,,, and also the values of e/m and E, the atomic 
weight of electron, were calculated. They are as follows: 


Ra=109,677.5795 +0.0037 cm=, 
Rp = 109,707.4226+0.0008 cm=!, 
Rue= 109,722.2694 +0.0024 cm-, 
= 109,737.3111+0.0041 
e/m= (1.75883+0.00022) X 10’e.m.u./g, 
E= (5.48740+0.00069) x 


For convenience, the result for E is also expressed in the 
following ways: 1/E = 1822.36+0.23, number of electrons 


1R. C. Williams, Phys. Rev. 54, 559 (1938). 


11. The Turbulent Jet. P. Y. Coou, Guggenheim Labora- 
tory, California Institute of Technology (Introduced by R. 
A. Millikan).—The first solution of the jet problem based 
upon the equations of turbulent fluctuation was given by 
C. C. Lin,’ who derived the so-called constant ‘‘exchange 
coefficient” from the equations of mean motion and of 
double correlations by using the values of the triple correla- 
tions at the center of the jet. The agreement between theory 
and experiment of the velocity and temperature distribu- 
tions for the plane and axial jets is satisfactory. This con- 
stant exchange coefficient has been assiimed recently by 
Gértler, based upon physical arguments.” Since 1941 the 
general theory of turbulence has been improved and the 


’ jet problem has again been attacked by S. S. Huang, and 


the solution agrees with Lin’s result for the plane and axial 
jets and with Tollmien’s for the half-jet under the mo- 
mentum transport theory. According to the latest version 
of the general theory, the constant exchange coefficient 
still follows as a first approximation from the equations of 
mean motion and of double correlations, substantiating 
Lin and Huang’s results. It is also shown that these equa- 
tions provide a basis for the more accurate solutions to. 


the problem. 


1 Science Report of National 
Ht f. angew. Math. u. Mech. 
*S. S. Huang, master’s thesis, National 


12. Radio Direction Finding at a Wave-Length of 1.8 
Meters. W. H. PICKERING, J. Davip, AND W. F. Hornyak, . 
California Institute of Technology.—The accuracy attain- 
able with an Adcock type antenna has been investigated. 
Particular attention has been directed to bearing errors due 
to polarization of the transmitted wave and to nearby ob- 
jects behaving as parasitic radiators. The antenna con- 
sisted of two vertical quarter-wave rods spaced three- 
quarters of a wave-length. The effects of various types of 
grounding plane have also been investigated. With favor- 
able conditions an accuracy of 0.2 degree is attainable, but 
polarization errors may be as large as 10 degrees. For best 
determination of the null direction, the receiver output is 
recorded on a paper strip synchronized with the rotation 


of the antenna, 


1941. 
Hua University, 1943, 
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